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program  for  the  prediction  of  Solid  Propellant  Rocket  Motor  Performance.  The  com- 
puter program  described  herein  will  be  referred  to  as  the  SPP  program,  and  describes 
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the  results  obtained  to  date. 
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Volume  III  of  this  report  is  a Program  User's  Manual  which  describes  the 
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to  Interpret  the  output.  A sample  case  Is  also  Included. 
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A flexible,  modular,  fully  automated,  solid  rocket  motor  performance  prediction 
program  has  been  developed.  The  program,  which  has  been  given  the  acronym 
SPP  is  based  on  six  pre-existing  computer  codes.  These  codes  have  been  Inte- 
grated and  modified,  as  required.  To  supplement  the  theory,  where  necessary, 
and  to  increase  the  flexibility  of  the  pregram,  a number  of  existing  and  newly 
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computed  as  a series  of  Independent  efficiencies.  Hie  program  currently  treats 
the  following  losses:  two-dtmensianal/two-phase  (coupled),  nozzle  erosion, 
kinetics,  boundary  layer,  combustion  efficiency,  submergence.  The  program 
predicts  average  delivered  perfonnar.ee,  as  well  as  mass  flow,  pressure,  thrust, 
impulse,  end  specific  impulse  as  functions  of  time  and  trajectory. 

In  order  to  assess  the  validity  of  the  SPP  program,  calculated  results  were  com- 
pared to  firing  data  for  four  different  types  of  motors.  While  conclusive  state- 
ments regarding  the  accuracy  and  range  of  validity  of  the  SPP  program  cannot  be 
made  until  additional  verification  efforts  are  conducted,  the  results  of  these  four 
test  cases  were  encouraging.  These  calculations  also  served  to  demonstrate 
the  desirability  of  eliminating  some  of  the  present  limitations  of  the  program. 
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PREFACE 


"his  Is  Volume  I of  a three-port  report  which  describes  a computer 
program  for  the  prediction  of  Solid  Propellant  Rocket  Motor  Performance.  The 
computer  program  described  herein  will  be  referred  to  as  the  SPP  program. 

Volume  I of  this  report  describes  the  engineering  analysis  which  was 
used  In  developing  this  computer  program  and  the  results  obtained  to  date. 

Volume  II  of  this  report  Is  a programming  document  of  the  computer 
program  which  was  developed  under  this  contract.  It  includes  a subroutlne-by- 
subroutlr.e  description  of  all  of  the  elements  of  the  SPP  program. 

Volume  III  of  this  report  Is  a Program  User’s  Manual  which  describes 
the  input  necessary  to  execute  the  SPP  computer  program  and  the  Information 
required  to  Interpret  the  output.  A sample  case  is  also  Included. 
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ABSTRACT 

A flexible,  modular,  fully  automated,  solid  rocket  motor  performance  pre- 
diction program  has  been  developed.  The  program,  which  has  been  given  the 
acronym  SPP  is  based  on  six  pre-existing  computer  codes.  These  codes  have 
been  Integrated  and  modified,  as  required.  To  supplement  the  theory,  where 
necessary,  and  to  Increase  the  flexibility  of  the  program,  a number  of  existing 
and  newly  developed  seml-emplrlcal  correlations  were  Incorporated  Into  the  pro- 
gram. The  program  has  a general  three-dimensional  grain  design  capability, 
coupled  to  a cne-dlmenslonal  ballistics  analysis.  The  deviations  from  ideal 
performance  are  computed  as  a series  of  Independent  efficiencies.  The  program 
currently  treats  the  following  losses:  two-dlraenslonal/two-phase  (coupled), 
nozzle  erosion,  kinetics,  boundary  layer,  combustion  efficiency,  submergence. 
The  program  predicts  average  delivered  performance,  as  well  as  mass  flow,  pres- 
sure, thrust.  Impulse,  end  specific  Impulse  as  functions  of  time  and  trajectory. 

The  theoretical  models  and  empirical  correlations  upon  which  the  program 
Is  based  are  described,  together  with  the  assumptions  that  are  employed.  The 
known  limitations  of  the  analysts  are  presented  to  aid  the  user  In  determining 
the  range  of  applicability  of  the  program. 

In  order  to  assess  the  validity  of  the  SPP  program,  calculated  results  were 
compared  to  firing  data  for  four  different  types  of  motors.  While  conclusive  state- 
ments regarding  the  accuracy  and  range  of  validity  of  the  SPP  program  cannot  be 
made  until  additional  verification  efforts  are  conducted,  the  results  of  these  four 
test  cases  were  encouraging.  Thtse  calculations  also  served  to  demonstrate 
the  desirability  of  eliminating  some  of  the  present  limitations  of  the  program. 
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Nozzle  area 
Propellant  bum  area 
Port  area 

Characteristic  exhaust  velocity 
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Mass  fraction  of  l^1  species 

Specific  heat  of  gas 

Specific  heat  of  propellant 
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Particle  diameter 

Thrust 

Gravitational  constant 
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Heat  of  fusion 
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Specific  Impulse 

Ratio  of  propellant  bum  area  to  nozzle  throat  area 
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Molecular  weight,  also  Mach  No, 

burning  rate  pressure  exponent 

Pressure 

Rate  of  change  of  pressure 
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Gas  constant 

Universal  gas  constant 
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Time 

Temperature 
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1 . INTRODUCTION 


The  physical  phenomena  which  govern  the  overall  efficiency  of  a solid  rocket 
motor  are  quite  numerous  and  complex.  For  this  reason  the  prediction  of  solid 
rocket  motor  delivered  performance  has,  historically,  been  empirically  oriented. 
Total  reliance  on  correlations  of  firing  data  Is  undesirable,  however,  as  their 
ranges  of  applicability  cannot  usually  be  clearly  defined,  and  a priori  performance 
predictions  for  new  motor  designs  or  propellants  cannot  be  made  with  confidence. 

In  order  to  achieve  the  desired  generality,  performance  predictions  should,  where 
possible,  be  based  on  physically  realistic  models  of  the  controlling  phenomena. 

Over  the  past  ten,  to  fifteen,  years  Individual  facets  of  the  performance 
prediction  problem  In  solid  rocket  motors  have  been  the  subject  of  considerable 
study.  Reference  i contains  a survey  of  the  state-of-the-art  up  to  1971.  Some 
of  the  more  recent  developments  In  this  field  (Including  parts  of  the  present  Inves- 
tigation) are  discussed  In  Reference  2.  Terse  description  of  rocket  industry 
analytical  performance  prediction  capabilities  as  of  July  1974  may  be  found  in 
Reference  3.  Some  recent  analytical  and  experimental  work  is  described  In  Ref- 
erences 68  and  69. 

When  this  effort  was  Initiated,  analytical  models  for  most  of  the  physical 
processes  governing  motor  performance  had  been  formulated,  and  computer  programs 
based  on  these  models  were  In  existence.  Unfortunately,  these  computer  programs 
were  developed  Independently;  each  program  addressing  a particular  aspect  of  the 
overall  oroblem.  Few  organizations  possessed  a large  enough  operational  subset 
of  these  programs  to  make  analytical  performance  predictions.  This  situation  pro- 
vided the  primary  motivation  for  the  present  effort. 

It  was  believed  that  the  existing  technology  was  capable  of  yielding  adequate 
performance  predictions,  but  It  was  felt  that  a systematic  and  efficient  procedure 
for  Implementing  this  technology  was  required  In  order  to  foster  Its  use  on  a in- 
dustry-wide basis.  Thus,  the  primary  goal  of  this  program  was  to  develop  an 
automated,  computerized,  methodology,  based  on  existing  methods,  which  would 
be  capable  of  predicting  delivered  specific  Impulse  to  within  +2%  and  delivered 
thrust  and  total  impulse  to  within  +5%.  The  method  of  approach  utilized  to  achieve 
this  goal  is  outlined  in  Section  Since  we  can  expect  continuing  advances  in 
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the  state-of-the-art,  the  Solid  Rocket  Performance  Program  (SPP)  was  modularly 
structured  so  as  to  be  able  to  easily  accommodate  them. 

The  theoretical  techniques  employed  In  the  SPP  program  are  succintly  de- 
scribed In  Section  3.  This  report  Is  not  meant  to  be  a detailed  dissertation  on 
solid  rocket  performance  prediction  methods.  It's  purpose  Is  to  Indicate  the  methods 
selected  to  meet  the  objectives  of  this  study,  the  reasons  for  their  selection,  their 
primary  assumptions;  and  limitations.  Detailed  descriptions  of  the  theoretical 
models  and  solution  techniques  may  be  found  In  the  pre-existing  source  material 
referenced  throughout  the  text.  However,  In  cases  where  existing  methods  were 
significantly  modified,  or  extended,  details  of  the  changes  are  provided  herein. 

Pre-existing  and  newly  developed  performance  loss  correlations  are  presented 
In  Section  4.  Certain  of  these  correlations,  e.g.  combustion  efficiency,  are  re- 
quired since  theoretical  models  for  all  of  the  pertinent  phenomena  have  yet  to  be 
developed;  the  others  have  been  Incorporated  Into  the  program  for  the  sake  of  In- 
creased flexibility. 

No  performance  prediction  methodology  would  be  complete  without  considering 
such  effects  as  bum  rate,  nozzle  erosion  and  particle  size.  The  manner  In  which 
these,  and  other  Items,  were  addressed  Is  described  In  Section  5. 

The  limited  comparisons  of  analytical  predictions  with  test  data  that  have 
been  carried  out  to  date  are  described  In  Section  6.  Conclusions  drawn  from  this 
study  are  presented  In  Section  7. 1.  Section  7.2  contains  a list  of  the  more  signi- 
ficant limitations  of  the  SPP  program  and  our  recommendations  for  future  efforts 
directed  tcvards  alleviating  them.  This  Is  followed  by  Section  8,  which  contains 
the  references  cited  In  the  text. 

As  noted  In  the  Preface  this  report  Is  Volume  1 of  a, three  volume  set.  Since 
this  volume  and  Volume  3 (User's  Manual)  will  often  be  used  by  different  parties 
some  of  the  descriptive  material  contained  herein  has  also  been  Incorporated  Into 
the  User's  Manual. 
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In  accordence  with  die  last  two  of  the  listed  program  characteristics  the 
Solid  Rocket  Performance  Program  (SPP),  was  modularly  structured;  each  module 
(with  the  exceotion  of  the  control  module)  having  been  based  on  a pre-existing 
computer  program.  The  following  programs  form  the  foundation  of  the  SPP  program: 

• One  Dimensional  Equilibrium  Program  (ODE)^ 

• One  Dimensional  Kinetics  Program  {ODKJ^ 

• Two  Dimensional  Two  Phase  Program  (TD2P)® 

(7) 

• Turbulent  Boundary  Layer  Program  (TEL)  ' 

(8) 

• Three  Dimensional  Grain  Design  Program 

• Motor  Ballistics  Program  ® 

Each  of  the  first  four  of  these  programs  Is  a separate  module  In  Jie  TPP  pro- 
gram. while  the  last  two  programs  have  been  combined  In  a single  .nad'iic.  These 
programs  have  been  combined  with  a control  module  which  permits  th  .an  to  be  run 
together  automatically.  The  six  modules  that  comprise  the  SPP  progra-i  are  sche- 
matically shown  In  Figure  2-1.  The  primary  function,  or  functions,  of  each 
module  Is  also  Indicated.  All  of  the  programs  which  were  adopted  for  ^se  In  the 
SPP  program  have  been  modified  to  some  degree  to  better  suit  the  present  purpose. 
These  modifications  are  discussed  In  Section  3,  as  part  of  the  module  by  module 
description  of  the  program. 

The  SPP  program  has  also  been  designed  to  provide  the  user  of  the  code 
flexibility  In  selecting  which  of  the  modules  are  to  be  exercised  In  any  one  com- 
puter run.  The  control  module  allows  all  jf  the  modules  and  programs  to  be  ex- 
ercised sequentially  In  a single  run.  Often,  however,  this  Is  not  the  most  efficient 
approach.  Therefore,  the  program  also  allows  the  performance  related  modules  to 
be  exercised  Individually,  or  In  any  "aer  determined  combination.  A complete  per- 
formance prediction,  utilizing  all  of  the  program  modules  can  be  carried  out  in 
parts  as  efficiently  (from  a computer  time  standpoint)  as  If  only  a single  run  were 
made.  This  Is  made  possible  by  the  nature  of  die  linkage  structure  built  Into  the 
program.  All  of  the  data  whit*  needs  to  be  transferred  from  one  module  to  another 
is  written  out  on  an  external  logical  unit  (In  addition  to  being  automatically  trans- 
ferred Internally).  At  the  completion  of  a run,  this  linkage  data  can  then  be  punch- 
ed, or  stored  In  some  other  permanent  manner.  This  data  Is  then  available  for  use 
in  subsequent  runs,  thereby  eilmi-v  >>e  need  for  redundantly  exercising  a module 


In  accordence  with  die  last  two  of  the  listed  program  characteristics  the 
Solid  Rocket  Performance  Program  (SPP),  was  modularly  structured;  each  module 
(with  the  exceotion  of  the  control  module)  having  been  based  on  a pre-existing 
computer  program.  The  following  programs  form  the  foundation  of  the  SPP  program: 

• One  Dimensional  Equilibrium  Program  (ODE)^ 

• One  Dimensional  Kinetics  Program  {ODKJ^ 

• Two  Dimensional  Two  Phase  Program  (TD2P)® 

(7) 

• Turbulent  Boundary  Layer  Program  (TEL)  ' 

(8) 

• Three  Dimensional  Grain  Design  Program 

• Motor  Ballistics  Program  ® 

Each  of  the  first  four  of  these  programs  Is  a separate  module  In  Jie  TPP  pro- 
gram. while  the  last  two  programs  have  been  combined  In  a single  .nad'iic.  These 
programs  have  been  combined  with  a control  module  which  permits  th  .an  to  be  run 
together  automatically.  The  six  modules  that  comprise  the  SPP  progra-i  are  sche- 
matically shown  In  Figure  2-1.  The  primary  function,  or  functions,  of  each 
module  Is  also  Indicated.  All  of  the  programs  which  were  adopted  for  ^se  In  the 
SPP  program  have  been  modified  to  some  degree  to  better  suit  the  present  purpose. 
These  modifications  are  discussed  In  Section  3,  as  part  of  the  module  by  module 
description  of  the  program. 

The  SPP  program  has  also  been  designed  to  provide  the  user  of  the  code 
flexibility  In  selecting  which  of  the  modules  are  to  be  exercised  In  any  one  com- 
puter run.  The  control  module  allows  all  jf  the  modules  and  programs  to  be  ex- 
ercised sequentially  In  a single  run.  Often,  however,  this  Is  not  the  most  efficient 
approach.  Therefore,  the  program  also  allows  the  performance  related  modules  to 
be  exercised  Individually,  or  In  any  "aer  determined  combination.  A complete  per- 
formance prediction,  utilizing  all  of  the  program  modules  can  be  carried  out  in 
parts  as  efficiently  (from  a computer  time  standpoint)  as  If  only  a single  run  were 
made.  This  Is  made  possible  by  the  nature  of  die  linkage  structure  built  Into  the 
program.  All  of  the  data  whit*  needs  to  be  transferred  from  one  module  to  another 
is  written  out  on  an  external  logical  unit  (In  addition  to  being  automatically  trans- 
ferred Internally).  At  the  completion  of  a run,  this  linkage  data  can  then  be  punch- 
ed, or  stored  In  some  other  permanent  manner.  This  data  Is  then  available  for  use 
in  subsequent  runs,  thereby  eilmi-v  >>e  need  for  redundantly  exercising  a module 
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This  type  of  linkage  has  two  major  advantages.  The  first  is  an  obvious 
savings  In  computer  run  time  for  parametric  studies  in  which  the  output  from  one  or 
more  modules  is  unchanged  from  case  to  case.  The  second  is  the  ability  to  use 
other  analyses  Of  appropriate)  by  'fooling'  the  SFP  code  into  accepting  them  as 
an  integral  part  of  its  own  computational  procedure.  This  allows  the  user  to  sub- 
stitute die  output  from  any  other  code  or  analysis  into  the  SPP  prefonnance  predic- 
tion procedure,  provided  It  is  m Je  to  conform  to  die  proper  format. 

In  conclusion,  the  SPP  program  may  be  described  as  a modular,  flexible, 
code  which  combines  theory  and  empiricism  to  provide  accurate  predictions  of 
solid  rocket  motor  delivered  performance.  The  analytical  techniques  and  correla- 
tions currently  utilized  in  die  SPP  are  summarized  in  Table  2-1.  In  many  cases  the 
program  allows  the  user  (an  option)  to  directly  input  a performance  related  quantity. 
In  these  cases  (indicated  In  die  table)  the  pertinent  analytical  and/or  empirical 
calculations  are  bypassed,  in  favor  of  the  input  quantity. 

2.2  Definition  of  tosses 

The  basic  method  of  approach  used  in  the  development  of  the  SPP  computer 
program  is  to  calculate  the  maximum  or  ideal  ‘theoretical’  performance  of  a rocket 
motor  and  to  subtract  from  that  maximum  performance  the  losses  or  deviations  from 
the  "theoretical"  performance  which  are  known  to  occur  in  a real  rocket  motor. 

In  Implementing  such  an  approach  It  13  often  assumed  that  the  various  losses  can 
be  treated  independently.  As  noted  In  Section  2.3,  all  of  the  losses  should  not  be 
treated  independently.  There  are  some  Interactions  strong  enough  to  warrant 
consideration.  With  these  exceptions  in  mind,  the  concept  of  treating  the  losses 
Independently  appears  to  be  valid,  considering  the  success  of  this  approach,  and 
others  of  a similar  nature. 

A standardized  performance  nomenclature  does  not  exist.  Therefore,  for  the 
sake  of  clarity,  the  definitions  used  herein  are  listed  below. 

• Performance  - unless  otherwise  specified,  performance 

will  be  in  terms  of  specific  Impulse  I , 
delivered  to  a vacuum.  sp 

•Theoretical  I - the  maximum  possible  which  can  be  de- 

v ltvered  at  the  initial  area  ratio  by  a propel- 

lant at  a given  chamber  pressure  and  total 
enthalpy.  Referred  to  as  I, 

spth 


i ' 


) 
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SUMMARY  OF  SPP  PROGRAM  CAPABILITIES 


• Loss 

• Two  phase  flow  loss 

• Two  dimensional  or 
divergence  loss 

•Finite  Rate  Kinetics 
loss 

•Boundary  Layer  Loss 
•Erosion  Loss 
•Submergence  Loss 
•Combustion  Efficiency 

•Delivered  Performance 


the  decrement  from  the  theoretical  l£p  which 
can  be  attributed  to  a physical  phenomena 
not  Included  In  the  calculation  of  I..  . 

spth 

the  decrement  In  performance  due  to  finite 
velocity  and  temperature  differences  be- 
tween the  gas  and  condensed  phase. 

the  decrement  In  performance  which  Is  due 
to  the  momentum  of  the  rocket  exhaust  not 
being  totally  aligned  with  the  axis  of  the 
motor. 

the  decrement  in  performance  which  is  due 
to  Incomplete  transfer  of  latent  heat  to  sen- 
sible heat  caused  by  the  finite  time  required 
for  gas  phase  chemical  reactions  to  occur. 

the  decrement  In  performance  due  to  viscous 
forces  adjacent  to  the  nozzle  wall  and  heat 
transfer  to  the  nozzle  wall. 

the  decrement  In  performance  associated 
with  erosion  Induced  changes  In  the  nozzle 
exit  area  ratio 

the  decrement  In  performance  that  can  only 
be  attributed  to  the  fact  that  a given  nozzle 
is  submerged. 

as  used  here,  the  term  combustion  efficiency 
will  refer  only  to  the  degradation  of  I due 
to  a departure  of  chamber  total  temperature 
from  the  theoretically  calculated  total  temper- 
ature. 

the  calculated  performance  of  the  rocket 
motor  which  Includes  all  of  the  losses  con- 
sidered herein. 


2.3  Calculation  of  Delivered  Performance 

Within  fte  framework  already  discussed  (l.e.,  separation  of  Individual 

losses),  there  are  three  basic  methods  which  can  be  used  to  calculate  delivered 

I . The  first  method  Is  to  calculate  discrete  decrements  In  performance,  aI 
sp  sp 

and  to  subtract  these  losses  from  the  theoretical  specific  Impulse,  I . That 
is  Ptn 


spD  spth 


- Tj  Al, 


sp. 


loss 


(2-1/ 
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The  second  method  is  to  calculate  a series  of  efficiencies,  n's,  to  be  applied 
multipllcatively,  l.e. 


^Pp  “ ^ " »>10S3  • (2_2) 

Both  of  the  above  methods  are  the  same  through  first  order,  thus.  If  the  losses 
are  small  they  yield  the  same  result. 

The  third  method,  selected  by  the  J ANN  AT  Committee  on  Performance  Stan- 
dardization for  liquid  rocket  engines.  Is  a combination  of  the  first  two  methods. 
That  Is,  all  of  the  losses  except  the  boundary  layer  loss  are  treated  as  efficiencies 
while  the  boundary  layer  loss  Is  treated  as  a decrement.  Hence, 

IspD  “ Ispth  ” 1,1055 ' 4IapBL  <2'3) 

This  latter  method  of  approach  was  selected  for  Incorporation  In  the  SPP  code  for 
two  reasons.  The  first  reason  being  that  there  was  no  compelling  reason  to  de- 
viate from  what  has  been  established  as  a “standardized"  method,  and  secondly, 
that  any  other  combination  of  the  above  methods  yields  substantially  the  same  end 
result. 

Unless  a grain  design  and  ballistics  solution  Is  obtained,  the  only  other 
performance  related  quantity  that  Is  calculated  Is  the  specific  Impulse  delivered 

to  average  ambient  pressure,  I , defined  as  follows: 

spD 

a 

i =:  -L. 

spD  spD  P 
a c 

where: 

P a average  ambient  pressure  (pst) 

Q 

Pc  = average  chamber  pressure  (psl) 
r = average  nozzle  expansion  ratio 
CJ,  = theoretical  C*  (ft/sec) 

CD  = nozzle  discharge  coefficient  (from  (TD2P) 


^C‘th 

-egr- 


(2-4) 
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If  a grain  design  and  ballistics  solution  is  obtained  delivered  thrust  (vacuum 

and  ambient),  delivered  specific  Impulse  (vacuum  and  ambient)  and  mass  flow  are 

calculated  as  functions  of  time.  The  I values  are  calculated  by  applying  equa- 

sp 

ting  (2-3)  and  (2-4)  at  successive  time  Increments,  Delivered  thrust  values  are 
then  computed  as  follows*, 

FD(t)  = m(t)  I (t)  (2-5) 

u spD 

FD  (t)=  FD(t)  - Pa(t)  Ae  (2-6) 

a 

The  total  mass  flow,  m^,  total  impulse  1^,  time  averaged  specific  Impulse,  , 
and  time  averaged  mass  flow  rate,  m are  also  calculated  in  the  indicated  manneP, 


mT  = 

Jlbum  m(t)  dt 
0 

(2-7) 

Ij  = 

fHjum  F_  (t)  dt 
o a 

(2-8) 

T*pd  = 

VmT 

(2— S) 

5 

m'/tbum 

(2-10) 

where  tjJum  denotes  the  total  motor  bum  time. 

2.4  loss  Interactions  and  Reference  Performance  Values 

The  concept  of  defining  a theoretical  I and  various  individual  losses  Is 

sp 

a convenient  artifice  that  allows  us  to  predict  delivered  preformance.  In  addition 
to  being  artificial,  this  concept  of  treating  individual  losses  is  somewhat  sub- 
jective. The  number  of  “losses"  that  are  defined,  and  what  they  are  called, 
can,  and  does,  vary  from  one  performance  prediction  methodology  to  another. 

There  is  a constraint  on  this  subjectivism,  however,  which  unfortunately  is  some- 
times violated.  One  way  of  expressing  this  constraint  is  as  follows:  any  number 
of  individual  losses  may  be  defined,  as  long  as  each  physical  source  of  non-idea! 
performance  is  accounted  for  once,  and  only  once. 

•If  a grain  design  and  ballistics  solution  is  obtained,  but  a two  phase  flow  solution 
is  not,  then  thrust  is  calculated  somewhat  differently.  See  Section  3.3.2  for  the 
related  details. 
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2.4.1  Loss  Interactions 

In  defining  a set  of  losses  one  normally  tries  to  arrange  for  each  loss  to 
be  primarily  determined  by  a single  physical  phenomenon  or  flow  characteristic. 

Since  the  phenomena  controlling  motor  performance  are,  in  reality,  a nonlinear 
mutually  interacting  set  a one  to  one  correspondence  between  phenomena  and 
losses  can  be  achieved  only  in  a relative  sense.  Each  of  the  "real",  non-ideal, 
processes  Influences  each  of  the  losses,  but  to  varying  degrees.  In  order  to  satisfy 
the  aforementioned  constraint  all  of  the  significant  Influences  of  each  phenomena 
should  be  accounted  for,  but  only  once.  The  losses  listed  In  Section  2.2  were 
defined  such  that  each  Is  primarily  the  result  of  a single  physical  phenomenon, 
or  effect.  For  example,  the  kinetics  loss  Is  chiefly  a result  of  finite  rate  kinetic 
effects.  Table  2-2  (from  Ref.  10)  represents  an  attempt  to  quantify  the  degree 
to  which  a given  phenomenon  Influences  the  losses  that  are  not  directly  related 
to  It.  Most,  but  not  all  of  the  losses  and  phenomena  are  listed  in  this  table. 

The  table  should  be  read  down  and  to  the  left.  The  number  ihat  appears  In  each 
square  (except  the  diagonals,  which  represent  the  Interaction  of  a loss  with  Its 
chief  determinant)  Indicates  the  relative  degree  of  Interaction  between  the  Indi- 
cated loss  and  phenomenon.  The  key  to  Interpreting  this  number  is  given  In  the 
legend  below  the  table.  The  Importance  of  the  various  Interactions  Is  thus  char- 
acterized In  terms  of  their  potential  for  Influencing  predicted  delivered  performance. 

In  developing  the  present  methodology  It  was  decided  to  directly  Incorporate 
only  those  Interactions  which  were  felt  to  be  of  primary  Importance.  It  can  be  seen 
from  T.-ble  2-2  that  there  are  two  such  Interactions.  These  are:  the  effect  of  non 
one-dlmenslonal  flow  on  the  boundary  layer  loss;  the  effect  of  multiphase 
flow  on  the  divergence  loss.  3 

Two  dimensional  effects  have  an  Important  Impact  on  the  nozzle  boundary  i 

layer  because  of  the  sensitivity  of  the  boundary  layer  to  free  stream  conditions.  j 

The  wall  streamline  In  a typical  nozzle  differs  substantially  from  the  one  dimen- 
sional analysts  streamline  for  the  nozzle,  particularly  In  the  nozzle  throat  region 
where  the  heat  transfer  rates  are  the  highest. 

The  effect  of  multiphase  flow  on  divergence  loss  Is  Important  because  of  die 
manner  In  which  heat  Is  transfered  from  the  particles  to  the  gas.  Because  of  their 


2-9 


Yfr-m  rtfiiLiHuj 


inertia,  particle  concentrations  will  be  greater  near  the  nozzle  axis  than  the  nozzle 
wall.  Consequently,  there  will  be  greater  heat  addition  to  the  axis  streamline 
than  to  the  wall  streamline.  This  Is  demonstrated  by  the  fact  that  two  dimensional 
two  phase  flow  calculations  for  a perfect  gas  show  a greater  maximum  velocity 
for  the  axis  streamline  than  for  the  wall  streamline.  Thus,  two  phase  flow  effects 
cause  a reduction  In  the  percentage  of  misaligned  thrust  (lower  divergence  loss). 

It  should  be  pointed  out  that  the  estimates  given  In  Table  2-2  are  based  on 
typical  motor  configurations  and  operating  conditions.  The  Importance  of  the  various 
Interactions  should  be  reexamined  In  extreme  cases.  For  Instance,  In  a motor 
operating  at  very  low  pressure,  with  a metallzed  propellant,  the  coupling  between 
Incomplete  combustion  and  the  kinetics  loss  could  become  significant. 

2.4,2  Reference  Performance  Values 

Some  of  the  Individual  loss  calculations  directly  yield  a performance  decre- 
ment or  an  efficiency,  e.g. , the  boundary  layer  loss.  The  other  performance  loss 

calculations  yield  I values,  e.g.  kinetics  loss,  2-D,  2-phase  loss.  It  Is  not 
sp 

desirable  to  directly  utilize  any  of  these  I values  In  the  calculation  of  delivered 

performance.  Rather,  these  calculated  I values  should  be  related  to  a reference 

sp 

value  of  I ; thereby  producing  either  an  efficiency,  or  a performance  decrement.* 

Properly  defined  reference  values  are  one  of  the  keys  to  successful  implemen- 
tation of  a performance  prediction  methodology  based  on  separation  of  losses. 

One  must  make  sure  that  the  reference  computation  Includes  all  of  the  same  physical 
phenoir  :oa  and  assumptions  as  the  related  loss  calculation,  except  for  the  phen- 
omenon, or  phenomena,  responsible  for  the  loss,  as  defined.  Ideally,  this  reference 
calculation  should  be  made  by  the  same  computer  code  which  Is  selected  to  make 
the  loss  calculation,  since  this  would  minimize  deviations  due  to  numerical  tech- 
niques, truncation  and  roundoff  errors.  However,  such  an  approach  Is  not  always 
practicable.  For  example,  both  the  equations  for  the  two  phase  flow  and  chemical 
kinetics  losses  become  "stiff”  In  the  equilibrium  limit  and  result  In  unacceptably 
long  and  difficult  calculatlonal  procedures.  Even  11  sophisticated  numerical  tech- 


*ln  the  present  methodology  the  use  of  efficiencies  is  preferred.  It  Is  felt  that 
efficiencies  better  compensate  for  the  omission  of  things  like  real  gas  effects. 
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% ' nlques  are  used  to  calculate  the  "near"  equilibrium  case,  extrapolation  to  the  full 

| j equilibrium  solution  introduces  errors  which  can  negate  the  advantage  gained  by 

using  the  same  computer  code. 

In  developing  the  reference  performance  calculations  used  in  the  present 
methodology  the  following  criteria  were  used: 

The  conceptual  differences  between  the  reference  and  loss  calculations 
should  Include  only  those  mechanisms  which  cause  the  loss. 

Computational  efficiency  should  be  considered,  without  sacrificing  accuracy 
and  consistency. 

The  modifications  to  existing  computer  codes  should  be  Kept  to  a minimum 
to  minimize  program  development  costs. 

The  adopted  procedures  are  described  in  Section  3 . 
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ANALYTICAL  METHODS 


I 


The  *heoretical  methods  employed  In  the  SPP  program  were  all  adapted  from 
existing  computer  programs.  Thus,  the  basic  theories  ard  computational  tech- 
niques have  all  been  previously  documented.  As  a result  of  the  number  of  analyses 
contained  In  the  present  methodology  no  attempt  was  made  to  combine  these  pre- 
existing documents  Into  a single  entity.  The  purpose  of  this  section  Is  to  indicate 
the  methods  selected,  the  reasons  for  their  selection,  their  primary  assumptions 
and  limitations.  Wherever  applicable,  details  of  current  modifications  and/or 
extensions  to  the  pre-existing  analyses  are  Included.  This  section  Is  structured 
along  the  lines  of  the  program,  i.e.  the  methods  are  discussed  In  relation  to  the 
module  In  which  they  are  utilized.  The  Interfaces  between  the  modules  are  also 
described.  While  It  does  not  include  analytical  methods,  in  the  usual  sense  of 
the  word,  the  Master  Control  Module  Is  described  first,  since  It  Is  a key  element 
in  the  overall  program  structure . 
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Master  Control  Module 


The  master  control  module  (MC)  controls  the  execution  of  the  SPP  computer 
code.  The  MC  module  selects,  via  user  Input,  which  of  the  five  basic  calculation 
modules  aie  to  be  executed,  and  whether  calculated,  input,  or  empirical,  losses 
are  to  be  .^elected  for  the  calculation  of  delivered  specific  Impulse,  thrust,  and 
total  Impulse.  This  module  also  controls  all  of  the  internal  communication  bet- 
ween modules. 


In  addition  to  its  control  functions  the  MC  module  stores  selected  results 
from  each  of  t*-e  performance  loss  modules  and  then  uses  them  to  calculate  delivered 
performance.  Table  3-1  shows  the  variables  used  to  calculate  delivered  perfor- 
mance, and  Indicates  the  module  In  which  each  Is  calculated. 

1 The  delivered  performance  calculation  does  not  use  all  of  the  variables  11st- 

- ed  In  Table  3-1  In  their  given  form. 
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Where  necessary  the  MC  mocule  combines,  ratios,  averages,  etc.,  these 
variables,  in  order  to  generate  the  desired  quantities.  The  empirical  performance 
loss  estimates  are  also  computed  In  the  MC  module. 
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Table  3-1.  Variables  Used  In  Delivered  Performance  Calculation 


\ 1 " ■■  - 

Module 

Variable 

Theoretical  Performance 

Spth  spREQ  1 

Grain  Design  and  Ballistics 

Pc,  m,  r*,  L*,  nc*  as  functions  of  t. 

Kinetics 

Ispkin  at  '»• 

TWo  Dimensional  TWo  Phase  Flow 

I , as  functions  of  c and  C_.  i 

spTD2P  *****  u ‘ 

Turbulent  Boundary  Layer 

£1 

sPtbl 

The  equations  used  to  calculate  delivered  performance  have  been  given  In 
Section  2.3.  The  equation  fcr  vacuum  delivered  specific  Impulse,  previously  written 
in  symbolic  form  (equation  2-3),  is  written  here  In  the  manner  used  In  the  program. 

'spD  = !spth  ’taN  ' t,TD2P  * ’’C.  * 'ISUB  * '!spm  (3_1) 


The  theoretical  specific  Impulse,  I 


Is  losed  on  the  Initial  (t=0)  nozzle 


exit  area  ratio. 


The  loss  due  to  finite  rate  chemical  kinetics  Is  treated  as  an  efficiency 
factor,  The  kinetic  efficiency  is  calculated  as  the  ratio  of  an  I calcula- 

ted by  the  CDK  module,  to  an  I calculated  by  a special  option  In  the  ODE  module. 

sp 

This  special  ODE  option,  referred  to  as  the  "Restricted  Equilibrium"  option  causes 
the  ODE  program  to  compute  an  I based  on  the  more  restrictive  physical  assump- 

Sp 

tlons  employed  In  the  ODK  calculation.  This  procedure,  which  Is  further  described 
In  Section  3.4,  Is  r.ecessary  if  certain  losses  are  not  to  be  counted  twice. 

If  a motor  experiences  throat  erosion  during  the  course  of  a firing  the  nozzle 
exit  area  ratio,  and  hence,  the  specific  impulse,  varies  with  time.  The  performance 
decrement  associated  with  this  phenomena  Is  usually  called  the  erosion  loss. 

While  this  loss  can  be  estimated  based  on  one-dlmenslonal  calculations.  It  was 
felt  that  a two-dimensional  estimate  of  this  effect  would  be  superior.  Thus,  for 
eroding  nozzles,  ls  calculated  as  follows: 


3-2 


’»TD2P“’%32P 
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The  quantity  t^q2p  ls  time  averaged  two  dimensional  two  phase  efficiency. 
The  TD2P  module  calculates  an  efficiency  as  a function  of  the  area  ratio  t . Area 
ratio  js  a function  of  time  Is  established  from  the  Input,  or  internally  calculated, 
nozzle  erosion  rate.  This  allows  nxD2 P to  *■**  calculated  as 


As  defined  In  equation  (3-1)  Vjjy2p  Includes  the  'erosion*  loss.  In  addition 

to  the  coupled  two-phase  and  divergence  losses.  The  erosion  loss,  as  defined 

herein.  Is  equal  to  the  ratio  of  the  2-D,  2-phase  I at  the  time  averaged  expan- 
— ®P 

slon  ratio,  1 (t),  to  the  specific  Impulse  evaluated  at  the  Initial  expansion 

“PTD?P 

ratio,  I ( (,).  (Future  modifications  to  the  program  should  break  the  erosion 

spTD2P  1 

less  out  as  a separate  efficiency.) 


The  combustion  efficiency,  jjqj.  Is  obtained  from  the  average  empirical  c* 
efficiency  as  follows: 


Vet  = C'D  Vc* 


(3-4) 


tjce  = 1.0,  If  the  previous  relation  yields  a value 
greater  than  unity 

The  c*  t.  clency  Is  a correlation  cf  experimentally  measured  c*  values  (see  Sec- 
tion 4.1)  and,  therefore.  Includes  all  of  die  physical  effects  which  determine  the 
nozzle  mass  flow  rate.  One  of  these  effects,  the  combined  two  dimensional,  two 
phase,  effect  on  moss  flow,  is  also  Included  In  the  calculation  Gf  Vye2P‘  T*’e 
definition  of  jjC£  In  equation  (3-4)  prevents  the  mass  flow  effect  related  to  the  nozzle 
discharge  coefficient,  C^,  from  being  counted  twice. 

It  is  well  known  that  the  discharge  coefficient  for  single  phase,  two  dimen- 
sional, flow  Is  less  than  unity.  However,  It  is  also  well  established  that  In  non- 
equilibrium  two  phase  flew  more  mass  flows  through  a nozzle  at  a given  chamber 
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pressure  than  would  If  the  particles  were  In  velocity  and  thermal  equilibrium.  By 
Itself  this  effect  would  produce  CQ  values  greater  Qian  unity.  For  most  metallzed 
(olid  propellant  motors  the  net  result  of  these  two  opposing  effects  is  a discharge 
coefficient  greater  than  unity.  Since  ijc*  is  empirically  obtained,  and  subject  to 
scatter,  the  second  part  cf  equation  (3-4)  Is  added  to  prevent  meaningless  values 
of  tiggfroa  being  calculated  (i.e.  t?cE>U- 

If  the  combustion  efficiency  were  100%  {complete  combustion)  the  c*  efficiency 
of  the  motor,  rjc*.  would  he  given  by 

nc.  = &-S) 

Thus,  in  such  cases,  tjc*  would  usually  be  less  than  unity  for  metallzed  propellants, 
and  would  always  be  greater  than  unity  for  non-metalized  propellants.* 

The  submergence  efficiency,  r)sus,  is  based  on  an  empirical  correlation  and 
is  a function  of  the  throat  radius,  mass  fraction  of  condensed  phase,  length  of  sub- 
mergence to  length  of  the  internal  motor,  and  chamber  pressure.  This  correlation 
is  discussed  in  more  detail  in  Section  4.2.6. 

The  decrement  in  performance  due  to  the  formation  of  a turbulent  boundary 

layer  in  the  motor  nozzle  is  transmitted  directly  to  the  MC  module  and  is  reported 

as  a function  of  area  ratio.  The  decrement  in  I Is  taken  to  be  the  value  corres- 

sp 

ponding  to  the  initial  nozzle  exit  area  ratio.  This  decrement  is  assumed  to  be  in- 
dependent of  erosion  induced  area  ratio  changes,  and  throat  roughness. 


Delivered  performance  (I  only)  is  calculated  somewhat  differently  when 
5P 

based  completely  on  the  performance  loss  correlations  presented  in  Section  4. 
Instead  of  equation  {3»1),  the  following  relation  Is  employed, 

= ’’KIN  ^2D  ^2t  Vcz  ’’SUB  " 4lsp&  (3~6 


•Most  solid  rocket  motors  have  normalized  radii  of  curvature  w2  and,  in  the  ab- 
sence of  particles  Cq  tends  to  be  .99  f Cp  < 1 , or,  in  other  words,  with  perfect 
combustion  ric*  would  be  bounced  as  follows,  1 <i?c*S  1.01.  This  small  effect 
tends  to  be  masked  by  measurement  and  data  reduction  errors,  and  many  times 
test  engineers  are  reluctant  to  report  values  cf  Jjc*  greater  than  unity  due  to  tho 
mistaken  assumption  that  such  values  are  physically  Impossible  to  achieve. 
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In  the  «i wyHfUtd  performance  prediction  methodology  the  two  dimensional, 
tr2D,  and  two  phase,  tj2?,  efficiencies  are  treated  separately.  The  combustion 
efficiency , T!q£,  Is  simply  equated  to  the  empirical  c*  efficiency,  rtc*!  end  the 
so-called  'erosion  loss*  Is  not  directly  accounted  for.  If  an  estimate  of  the  ero- 
sion Ion  Is  desired,  and  a TD2P  solution  has  not  been  obtained,  one  can  compute 


the  ratio  I. 


(f{t"0)},  where  I {«  (t=0)}  Is  the  specific  Impulse  at 


the  Initial  nozzle  expansion  ratio,  usually  denoted  simply  by  I 


In  order  to 


compare  the  delivered  specific  Impulses  obtained  theoretically  (equation  (3-1)), 
and  empirically  (equation  (3-6))  an  estimate  of  the  erosion  loss  must  be  addend  erf 

to  (3-6),  or  equivalently,  I must  be  replaced  by  1 (7). 

£pth  spth 

the  Master  Control  Module  also  calculates  the  other  performance  related 

quantities  of  Interest  I , Fn  ft),  S~  etc. , using  equations  (2-4)  to  (2-10). 
spD  °a 

a 

3.2  Theoretical  Performance  Module 


For  a specified  propellant  combination  the  maximum  specific  impulse,  Isp  , 
which  can  possible  be  delivered  can  be  predicted  using  any  of  a number  of  01 
standard  computer  codes  developed  for  this  purpose.  A survey  conducted  by  the 
JANNAF  Performance  Standardization  Working  Group  found  sixteen  such  computer 
programs,  each  of  which  represented  an  Independent  development.  All  of  these 
programs  are.  In  principle,  capable  of  calculating  the  quantities  of  Interest,  i.e. 
shifting  and  frozen  I and  Cp,  exhaust  species  etc.  However,  many  of  these 
programs.  In  practice,  cannot  adequately  handle  the  solid  phases  required  In 
solid  rocket  motor  applications.  Of  the  remaining  applicable  programs,  two  stand 
out  by  virtue  of  their  excellent  characteristics.  These  programs  are: 

The  ODE  Computer  Program  developed  by  NASAARC  and 
described  in  NASA  SP  273,  Reference  4. 

The  Theoretical  igp  Program  developed  by  the  AFRPL  and 
described  In  Reference?  10  and  11. 

with  either  computer  program  area  ratio  expansion  calculations 
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yielding  theoretical  performance  can  be  carried  out  for  any  rocket  engine  of 
current  or  projected  use  provided  basic  thermodynamic  data  exists  for  the  pro- 
pellant. Each  computer  program  currently  contains  a vast  library  of  thermodynamic 
data.  A comparison  of  these  two  computer  programs  as  applied  to  the  proposed 
effort  is  presented  in  Hie  paragraphs  below. 


m 


the  OOE  Program 

This  computer  pregram  has  been  under  continuous  development  at  NASA/IRC 
since  the  mid  19S0's.  In  1968  it  was  adapted  by  the  ICRPG  (currently  JANNAF) 
Performance  Standardization  Working  Group  as  the  standard  computer  program  for 
theoretical  performance  computation.  A newer  extensively  revised  version  of  ODE 
is  now  available.  This  version  is  documented  In  complete  detail  in  Reference  4. 

Older  versions  of  ODE  were  unable  to  successfully  compute  equilibrium 
compositions  fer  aluminized  propellants.  This  difficulty  lead  directly  to  the  de- 
velopment of  certain  of  the  sixteen  chemical  equilibrium  codes  previously  men- 
tioned. The  more  recent  versions  . the  ODE  program,  however,  are  all  capable 
of  handling  multiple  solid  phase  products.  For  example,  the  SP273  version  of 
ODE  can  successfully  compute  compositions  with  large  amounts  of  solid  phase 
carbon;  a system  with  a history  of  difficulties. 

The  ODE  program  uses  a free  energy  minimisation  technique,  along  with  the 
following  assumptions,  to  calculate  the  vacuum  I and  related  theoretical  per- 
formance parameters. 

One  dimensional  flow. 

Chemical  equilibrium  between  Species  and  phases. 

Thermal  and  dynamic  (velocity)  equilibrium 

between  condensed  and  gaseous  phases. 

Adiabatic,  Isertroplc  flow. 

The  ODE  program  requires  curve  fit  data  for  enthalpy  vs.  temperature.  Seven 
coefficients  {Sj,  a2«  2j,  a^,  a5,  ag,  a^)  are  required  such  that: 


CP 

pT 


al  + ®2T  + ®3P  + aiP  + aST* 
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aj  + a2T/2  + a3T»/3  + a^/A  + agTVS  + ag/T 


9j  InT  + a2T  + agT3/2  + a4T3/3  + asT*/4  + ay 


Two  set*  of  coefficients  are  required  per  species  such  that: 

300 °K sis  1000°K  for  set  1 

1000 °K  s T s 6000°K  for  set  2 

The  fits  must  be  continuous  across  the  Juncture  at  1000°K.  NASA  SP273 
contains  coefficients  for  421  chemical  species.  A NASA/LRC  program  exists  for 
the  purpose  of  calculating  thermodynamic  data  and  outputting  the  above  curve 
fit  data.  This  program,  called  PAC,  Is  described  In  Reference  23.  A new  all 
FORTRAN  IV  version  of  PAC  Is  now  available  and  has  been  found  to  be  highly  sat- 
isfactory. 

The  curve  fit  data  Is  highly  compact  and  computationally  efficient  but  suffers 
from  two  major  draw  backs,  which  are: 

• certain  data  cannot  be  accurately  fit  with  only  seven  coefficients, 
especially  over  the  range  300°K  to  1000°K. 

• curve  fit  data  Is  awkward  to  produce  since  It  requires  use  of  the 
PAC  program  followed  by  careful  comparison  of  the  resulting  fit 
to  the  original  data. 

NASA  SP273  also  contains  elemental  composition  and  enthalpy  data  (l.e., 
reactant  cards)  for  62  commonly  used  propellants. 

Options  available  with  the  ODE  program  are: 

Point  Options: 

• pressure  and  temperature 

• enthalpy  and  pressure 

• entropy  and  pressure 

• density  or  volume  and  enthalpy 


Problem  Options: 

• rocket  option  (l.e.  chamber,  throat,  assigned  area  ratios,  and/or 
pressures) 
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AfRTI.  Theoretical  I Program 

This  program  also  uses  a free  energy  minimization  technique  to  calculate 

theoretical  performance,  and  employs  the  same  assumptions  listed  for  the  ODE 

program.  The  Aran,  program  contains  several  options  not  found  In  ODE,  e.g. 

Air  Augmented  I , Delta  Velocity,  Weight  Optimization;  however,  these  options 
sp 

are  not  of  Interest  for  the  present  purposes. 

The  Aran,  program  requires  curve  fit  data  for  specific  heat  vs.  temperature. 
Ten  coefficients  (Zj,  z^,  Zj,  z^,  zg,  zg,  z?,  zg,  Zg,  ZjQ)  are  required  such  that 

CP  " 21  + z2T_1  + Zj1'2  + Z4T"3  + 25T'4 

for  298°KsTsl200°K 
and 

Cp  - *6  + *7T  + *8T'*  «9T>  + *10T* 

for  1200°K  sir  6000°K. 

A curve  fit  program  and  an  extensive  library  of  fit  thermodynamic  data  Is 
maintained  at  AfRPL. 

Either  one  of  these  two  programs  would  be  very  satisfactory  for  the  purpose 
of  computing  theoretical  performance . Computer  run  times  for  both  programs  are 
short  and  are  In  no  way  a limiting  factor.  Computer  core  storage  requirements  are 
also  modest  (approximately  32K  words)  fci  both  programs.  Both  programs  feature 
excellent  input  methods.  The  ODE  program  was  chosen  for  incorporation  into 
the  performance  prediction  methodology.  The  deciding  factors  were  as  follows: 

1.  The  AFRPL  program  lacked  the  extensive  documentation  provided  with 
the  ODE  program. 

2.  The  ODE  program  has  been  adopted  as  a standard  by  JANNAF  and  Is  widely 
distributed. 

3.  The  ODE  program  was  already  Integrated  Into  the  Or.e  Dimensional  Kinetics 
Program  (ODlQ , which  was  adopted  for  the  kinetics  loss  calculation.  Re- 
placing ODE  with  another  program  would  have  required  a significant  effort 
to  link  the  new  program  to  ODK 
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3,2.2  Modifications  TO  The  ODE  Program 


For  the  present  purposes  the  ODE  program  was  modified  to  calculate  the 
transport  properties  required  by  the  Grain  Design  and  Ballistics  Module,  2-D,  2 
Phase  Module  and  Boundary  Layer  Module.  In  addition,  the  program  has  been 
modified  to  calculate  the  reference  condition  for  the  loss  In  performance  due  to 
finite  rate  chemical  kinetics.  This  calculation  Is  referred  to  as  the  "restricted 
equilibrium’  solution  and  Is  discussed  In  more  detail  In  section  3.4.  ODE  has 
also  been  modified  to  perform  an  equilibrium  calculation  at  the  exact  area  ratio 
corresponding  to  solidification.  To  reduce  the  length  of  the  program  the  subroutines 
not  required  for  the  performance  calculation  were  deleted.  Thus,  only  the  so- 
called  rocket  problem  option  of  the  ODE  program  can  be  exercised  In  the  present 
context. 

The  transport  properties  calculated  by  the  modified  ODE  program  are  the 
gas  phase  viscosity  u , thermal  conductivity,  y , and  Prandtl  number,  Pr.  The 
viscosity  and  thermal  conductivity  of  the  individual  gaseous  species  are  calculated 
from  formulas  given  In  Reference  20: 

4.15822xl0'8,  £ 7t 


(.45  + 1.32 


The  required  Leonard  Jones  parameters,  Or^.f^)  are  Internally  stored  In  the  computer 
program  for  206  species.  (See  Subroutine  MUK,  Volume  II , for  a list  of  the  species.) 
The  vlccoslty  of  the  mixture  Is  calculated  from  Wilke's  seml-emplrical  formula^21) 


« =S  f 
1=1 1 


“1 (1  % *11 

1/1 
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where  N Is  the  number  of  species,  X^,  the  mole  fraction  of  species  1,  and  Is 
defined  by: 

, f M.  l'1/2f  U.  1/2  M,  I/4I2 


, M w 1/2  M.  1/4  2 

*tf  = p72  1 + ^_  1 + ((i7)  (M j-5  (3"10) 

The  thermal  conductivity,  *,  Is  based  on  the  equation  given  by  Mason  and  Saxena 
which  Is  a slight  modification  of  Eucken*s  relation. 


N N x.  _1 

= S »,  1 + 1.0652  -1- 

1=1  1 ]=1  xl 

1^1 


The  Prandtl  number  Is  simply  given  by 


A subroutine,  0MEGA,  which  performs  a least  squares  fit  for  viscosity  In  the  form 
(T/T*)whas  also  been  added  to  ODE.  This  routine  was  necessary  because 
the  TDZP  and  TBL  programs  require  viscosity  to  be  In  this  form.  The  curve  fits 
treat  the  chamber,  throat  and  exit  points;  the  throat  values  are  treated  exactly. 

3.3  Grain  Design  and  Ballistics  Module 

3.3.1  Grain  Design  Calculations 

3.3. 1.1  Approach 

It  vms  required  that  the  SPP  program  be  capable  of  calculating  solid  rocket 
motor  delivered  performance  (l.e. , pressure  and  thrust  versus  time  as  well  as 
specific  Impulse  and  total  Impulse)  for  a wide  range  of  motor  designs  and  grain 
designs.  At  a minimum, end-burning,  circular  port,  star  port  and  spherical  geo- 
metries, and  at  least  two  other  qeometrles  representative  of  modem  production 
motors  (e.g.,  slotted  tube  and  flnocyl),  were  to  be  considered. 

The  problem  with  constructing  separate  specialized  subroutines  for  each 
geometry  Is  that  It  Is  both  cumbersome  and  limited.  A multiplicity  «f  Individual 
geometry  subroutines  would  add  greatly  to  program  length.  The  resultant  program 
would  also  be  limited,  since  the  list  of  potential  grain  geometries  is  exhaustive. 

In  addition,  the  geometry  calculations  do  not  stand  alone;  they  have  to  be  Integra- 
ted with  the  ballistics  calculations.  v separate,  specialized  geometry  cal- 
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culatlon  would  have  to  be  Individually  Integrated  with  the  ballistics  calculations. 
Clearly,  a general  method,  that  accomodates  a wide  range  of  geometries.  Is  a 
more  desirable  means  of  fulfilling  the  objectives  of  this  program. 

Several  grain  geometry  computer  programs  of  more  or  less  general  capability 
are  available.  Some  are  limited  to  the  geometry  problem,  some  Incorporate  com- 
bustion and  gasdynamtcs  to  provide  an  Integrated  Interior  ballistics  capability. 
These  candidate  programs  were  reviewed  In  order  to  select  the  most  appropriate. 
Selection  criteria  were  principally  geometric  capability  and  adaptability  to  general 
Interior  ballistics  analysis. 

3. 3. 1.2  Review  of  Existing  Computer  Programs 

3. 3. 1.2.1  Aerojet  Basic  Grain  Design  and  Interior  Ballistics  (564) 

This  program  solves  the  grain  geometry  and  Internal  ballistics  problems 
In  a sequential  manner.  The  geometric  solution  Is  two-dimensional,  and  the  In- 
| temal  ballistics  uses  a simple  burning  rate  law  and  cne-dlmenslonal  gasdynamtcs. 

Two-dimensional  grain  design  can  be  very  useful  for  symmetrical  motors 
with  sufficiently  large  length/diaraeter  ratio  that  end  effects  become  negligible. 

- However,  at  lower  length/diameter,  or  with  end  surfaces  burning,  or  with  partial 

bum  surface  restrictions,  even  a nominal  two-dimensional  geometry  would  re- 
quire a three-dimensional  treatment.  Consequently,  this  program  was  not  selected. 

3.3.1.  '..2  Lockheed  General  Grain  Design  Program  (470) 

This  grain  geometry  program  Is  basically  two-dimensional,  but  can  accom- 
modate cylindrical  grain  designs  In  three  dimensions.  An  Important  feature  of  this 
program  Is  that  It  includes  logic  to  account  for  nonuniform  burning  for  ballistics 
adaptability.  However,  It  was  not  selected  because  it  is  not  fully  capable  of  three- 
dimensional  treatment. 

3. 3. 1.2. 3 Thiokol  Generalized  Grain  Design  Program  (AGDA) 

The  AGDA  program  is  a three-dimensional  grain  geometry  program,  but  Is 
not  capable  of  handling  Internal  ballistics.  The  program  uses  the  Elkton-Moore 
method  which  Is  a vectorial  approach  to  the  geometric  problem.  It  Is  capable  of 
handling  a broad  range  of  grain  oo^igns,  the  complexity  manifesting  Itself  In  the 
effort  required  for  Input. 
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This  program  was  not  selected  because  It  does  not  Incorporate  combustion 
and  gasdynamlcs  analysis,  and  does  not  contain  logic  for  coupling  such  analysis 
to  the  geometry  calculations. 

3. 3. 1.2. 4 Boeing  Internal  Ballistics  Computer  Program  (D2-125286) 

The  Boeing  program  Is  a three-dimensional  grain  design  program  Incorpor- 
ating one-dlmenslonal  gasdynamlcs  and  erosive  burning.  However,  it  was  created 
for  a special  purpose  and  Is  therefore  deficient  In  some  respects  and  excessive  In 
others.  It  Is  deficient  In  that  It  Is  restricted  to  certain  types  of  grain  geometries. 
It  Is  excessive  In  die  detail  given  to  transients  and  high  acceleration  effects. 
Consequently,  the  program  as  presently  constructed  would  not  be  best  suited  for 
the  present  purpose. 


3. 3.1. 2. S Hercules  Grain  Design  and  Internal  Ballistics  Evaluation  Program  (64101) 

This  program  consists  of  two  parts,  the  Basic  Grain  Design  Program  and 
the  Internal  Ballistics  Evaluation  Subprogram. 

The  Basic  Grain  Design  Program  examines  a grain  In  three  dimensions  by 
analytical  simulation  of  drafting  techniques  used  In  developing  a grain  design. 

It  solves  for  geometric  parameters  used  In  Internal  ballistics  calculations,  and 
also  for  mass  properties  used  In  flight  simulation  calculations.  It  was  developed 
In  order  to  handle  the  most  complex  geometries  used  In  modem  production  motors, 
and  therefore  has  the  desired  general  capability. 


The  Internal  Ballistics  Evaluation  Subprogram  uses  simple  one-dlmenslonal 
gas-dynamics  and  erosive  burning.  It  uses  the  geometry  data  computed  for  each 
burning  station  In  the  Grain  Design  Program,  and  feeds  back  the  ballistics  data 
at  that  location.  It  therefore  accommodates  nonuniform  burning,  and  contains 
logic  for  coupling  grain  design  with  Interior  ballistics.  Ignition  and  taUoff  tran- 
sients are  computed,  but  more  for  appearance  than  accuracy. 

This  program  was  selected  because  of  Its  geometric  generality  and  bal- 
listics adaptability.  TWo  conectlble  deficiencies  were  detected:  the  method 
of  computing  the  radial  slots  In  segmented  grains,  and  the  enormous  size  of  the 
program  (largely  due  to  subroutines  not  required  for  the  present  purpose). 
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3. 3. 1.3  The  Selected  Program 


The  Hercules  64101  program  Is  described  In  detail  in  Reference  8.  'With 
all  of  its  options,  it  consists  of  six  links  and  sixty  subroutines.  Only  a fraction 
of  these  subroutines,  fifteen  in  all,  were  found  to  be  necessary  for  the  requirements 
of  the  present  program.  Most  of  these  are  contained  within  LINK  2 of  the  original 
program,  which  deals  with  the  geometry  calculations.  LINK  4 (a  pseudo  ballistics 
option)  and  LINK  S (plotting  routines)  were  eliminated  since  they  could  not  be  ex- 
ercised when  LINK  3 (the  ballistic-  - ubroutlnes)  is  used.  The  mass  properties 
subroutines  within  LINK  2 were  also  eliminated.  LINK  3 was  modified  to  incor- 
porate the  LPC  ballistics  analysis  described  in  Section  3.3.2,  and  LINKS  0 and 
1 (main  driver  and  common  regions  between  LINKS  2 and  3)  were  modified  in  ac- 
cordance with  the  other  insertions  and  deletions.  These  modifications  resulted 
in  an  80  percent  reduction  in  the  si2e  of  the  program.  Information  regarding  the 
operation  of  this  modified  program  is  presented  in  the  Users'  Manual  portion  of 
this  report.  Volume  3. 

Basically,  the  geometry  calculations  compute  volumes,  and  changes  in 
volumes.  Other  geometric  quantities  are  derivable  from  these  volumes.  Thus, 
for  example,  the  burr,  area  is  equal  to  the  rate  of  change  of  volume  divided  by 
the  rate  of  change  of  web: 


m 

= °PV 

(3-13) 

dm 

■ar 

= „ A.  *£. 
°p"b  dt 

(3-14) 

_ dV/dt 
dw/dt  ~ 

4MJ.I 

fiw/6 1 'mean  value 

(3-15) 

m = propellant  mass 

m = propellant  mass  flow  rate 

9„  = propellant  density 
P 

Ajj  = propellant  bum  area 
r = propellant  bum  rate 
w = propellant  web 
t = time 

V = propellant  volume 
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The  changes  In  web  over  a time  increment  are  calculated  from  the  ballistics 
subroutine,  and  the  corresponding  changes  in  volume  are  calculated  from  the 
geometry  subroutines.  Tha  computations  Iterate  on  flow  rate  until  a solution  for 
pressure  converges  Bum  area  is  output  as  a matter  of  information,  but  is  not 
used  as  such  In  the  computations.  Volumes  form  the  basis  of  the  computations  In 
order  to  best  satisfy  the  constraint  that  the  weight  expended  must  equal  the  weight 
loaded. 

Inputs  to  the  geometry  calculations  initially  consist  of  the  outer  dimen- 
sions of  the  propellant  grain,  which  is  taken  to  be  completely  filled  with  propellant. 
Next,  a series  of  inputs  subtracts  all  initial  void  volumes  representing  the  motor 
cavity.  Such  volumes  may  be  represented  by  combinations  of  cylinders,  cones, 
spheres,  and  prisms  located  by  spatial  coordinates;  these  volumes  may  Intersect 
each  other,  but  must  not  Intrude  into  propellant.  Supplementary  parameters  are 
available  to  define  an  Inhibited  surface,  a comer  round,  a segmented  grain  slot 
face,  and  the  symmetry  of  the  grain.  The  grain  Is  then  divided  Into  an  X-Y  com- 
putational mesh.  Web  fractions  are  also  defined.  When  the  web  burns  past  a 
specified  fraction,  information  may  be  output  or  Internally  transmitted  between  sub- 
routines. The  finer  the  mesh,  the  greater  the  accuracy,  the  smoother  the  output 
and  the  longer  the  computer  time.  The  method  by  which  the  changing  volumes  are 
calculated  for  each  mesh  segment  Is  presented  in  Ref . 8 and  will  not  be  repeated 
here. 

3. 3. 1.4  Program  Modifications 

As  far  as  the  geometry  calculations  per  se  are  concerned  the  only  modifi- 
cation required  was  In  the  method  of  computing  the  volumes  In  the  vicinity  of 
opposed-face  radial  slots  separating  segmented  grains.  In  the  course  of  using 
the  Lockheed  156-5  motor  (a  multiple  segmented  motor)  as  an  early  test  case, 
there  were  pressure  surges  in  the  ballistics  output.  These  were  traced  to  exces- 
sive volumes  being  expended  when  the  slots  crossed  the  reference  bum  Increments. 
In  this  type  of  situation,  the  calculations  were  expending  triangular  wedges  of 
propellant  (formed  by  the  end  face  and  the  port  to  a certain  length)  rather  tnan 
orthogonal  webs.  The  error  was  magnified  by  the  multiplicity  of  the  segments. 

This  deficiency  was  corrected  by  providing  Input  logic  to  key  the  existence  of  a 
slot,  and  a single  equation  for  the  annular  geometry  of  the  slot  to  be  used  at 
that  location. 
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A related  problem  was  encountered  with  end-burning  faces  In  general. 
Depending  upon  the  selection  of  the  mesh  size  and  bum  Intervals,  It  Is  possible 
to  have  the  same  volume  expended  counted  several  times.  This  can  be  avoided 
by  maintaining  a certain  relationship  between  these  Inputs,  as  discussed  In  the 
User's  Manual.  No  program  modification  was  required. 

The  aforementioned  80%  size  reduction  from  the  original  Hercules  program 
was  successfully  accomplished.  Care  had  to  be  exercised  to  assure  that  deletions 
would  not  destroy  functionality.  The  retained  subroutines  are  listed  In  the  Pro- 
gram Description,  Volume  2 of  this  report. 

The  revised  ballistics  analysis  Is  discussed  subsequently. 


3.3.2  Ballistics  Calculations 


3.3.2. 1 Approach 

The  Interior  ballistics  analysis  solves  the  conservation  equations  applied 
to  the  rocket  motor  cavity  to  predict  pressure  and  mass  flow  as  a function  of  time. 

In  its  most  rigorous  form,  a ballistics  calculation  solves  the  time  dependent 
equations  governing  conservation  of  mass,  momentum  and  energy.  The  simplest 
type  of  ballistics  calculations,  widely  used  In  rocket  motor  design,  employ  a 
steady-state  mass  balance.  The  sophistication  required  Is  determined  bv  the 
purpose  of  the  analysis  and  the  accuracy  required. 

A comprehensive  model  of  the  energy  exchange  processes  within  a rocket 

motor  was  prepared  by  LPC  for  the  prediction  of  duty  cycles  of  controllable  solid 
(13) 

rocket  motors  . The  sophistication  was  warranted  by  the  emphasis  on  transient 
phenomena:  predictions  of  Ignition  delay  times,  ignition  pressure  spikes,  oxtlngulsh- 
ment-relgnltion  phenomena  and  transient  performance  Integrals  were  required.  The 
program  would  be  recommended  for  such  purposes.  For  present  purposes  however, 
which  emphasize  steady-state  performance  prediction.  Its  size,  complexity  and  run 
time  would  not  be  justified. 

Computer  programs  such  as  Hercules  64101,  Aerojet  564  and  Lockheed 
241  use  a simplified  one-dlmenslonal  treatment.  Although  the  Hercules 
64101  program  was  selected  for  Us  geometric  capabilities.  It  would  not  be  re- 
commended for  the  computation  cf  Interior  ballistics.  It  employs  a three-stage 
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treatment  with  separate  analyses  for  Ignition,  steady-state  and  tailoff.  The  method 
of  interfacing  these  stages  is  ill-defined.  The  transient  stages  do  not  employ  an 
interna'  flow  analysis,  and  the  steady-state  stage  uses  a simplified  one-dimensional 
flow  solution  with  time  dependent  terms  omitted.  A better  approach,  to  avoid  ar- 
bitrary division  of  the  motor  operation  into  stages,  would  be  to  employ  one  set 
of  equations.  The  Seeing  program  takes  this  approach,  but  like  Reference  13  was 
developed  fer  a special  purpose. 

The  most  recent  interior  ballistics  analysis,  completed  by  Lockheed  for 

(9) 

the  Air  Force,  resulted  in  the  Lockheed  637  program*  . This  program  couples 
a comprehensive  ballistics  analysis  to  a limited  grain  geometry  capability.  How- 
ever, the  methodology  is  compatible  for  use  with  the  Hercules  geometry  analysis 
for  more  general  capability.  Therefore,  this  analysis  was  selected  as  best-suited 
to  the  requirements  of  the  present  program.  The  effort  required  was  to  reconstruct 
LINK  3 of  the  original  Hercules  program,  using  the  ballistics  analysis  from  Ref.  9 
as  the  basts.  The  original  designation  for  the  ballistics  subroutine  In  the  Hercules 
program,  MAIN3,  was  retained, 

3. 3. 2. 2 Flow  Equations 

The  equations  for  one-dlmenslonal  flow  used  tr.  the  computer  program  are 
based  on  the  following  assumptions: 

Particles  occupy  negligible  volume  compared  to  the  gas. 

Particle  velocity  and  temperature  are  equal  to  the  respective  gas  values. 

TV  oartlcle-gas  mixture  behaves  as  an  ■equivalent"  perfect  gas. 

Mass  is  generated  from  the  propellant  without  an  axial  component  of  vel- 
ocity; although  Invalid  for  end-burners , the  error  will  be  negligible  be- 
cause of  the  large  flow  area  downstream  of  the  end  face. 

Transient  effects  ure  simply  approximated. 

With  these  assumptions,  the  consecration  equations  are: 

P = eRT  (3-16) 

iM  + ifjA-*;,  (3-17) 
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ra 

e 

(3-18) 

h2g 

•cpTr 

(3-19) 

P = pressure,  llv'ft* 

p *»  gas  density,  lb/if 

R « gas  constant  for  the  products,  P/F 

T »•>  gas  temperature,  °F 

t - tine,  set 

A = flow  area,  ft? 

x - axial  distance,  ft 

mp  = incremental  mass  flow  rate,  lb/ft-sec 

u = gas  velocity,  ft/sec 

g » gravitational  constant,  ft/sec* 

Cp  = gas  heat  capacity,  ft/°  B 
Tp  m flame  temperature,  °P. 

These  equations  are  applied  at  each  computational  increment,  as  defined  by  Hie 
mesh  for  the  grain  geometry  calculation. 

The  cross-sectional  areas  are  supplied  by  the  grain  geometry  calculations. 
The  local  mass  flow  Increments,  m^,  are  determined  by  the  grain  geometry  and 
the  local  burning  rate  (sr»  Section  5.1  for  description  of  burning  rate  calculation). 

At  each  time  increment  a solution  is  obtained  In  the  following  manner. 

The  he-'  .-end  pressure  is  estimated,  and  equations  3-16  to  3-19  are  Integrated 
down  the  port,  step  by  step.  When  the  end  of  the  grain  Is  reached  the  nozzle 
choking  constraint  Is  tested.  The  procedure  Is  Iterated,  with  successive  head- 
end  pressure  estimates  until  the  choking  constraint  is  satisfied.  The  choking 
constraint  is  applied  as  follows.  First,  a throat  stagnation  pressure  Is  defined 
as: 

?t  * T7^m~  (3-20) 


rT 

W(HN) 

C* 

AP(NN) 


throat  stagnation  pressure,  lb/ln2 
mass  flow  through  the  throat,  lb/sec 
flow  characteristic  velocity,  ft/sec 
throat  area,  in? 
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Starting  with  ?„  and  a Mach  number  equal  to  one  at  the  throat,  the  lasntropic  re- 
lations are  used  to  calculate  what  the  static  pressure  at  the  nozzle  entrance  must 
be  1£  the  throat  Us  to  be  choked.  This  pressure,  P(NN),  Is  compared  to  the  static 
pressure  resulting  from  the  Integration  down  the  port,  P(N).  Convergence  (choking) 
is  realized  when  these  two  pressure  are  within  a specified  tolerance,  the  C*  used 
In  equation  3-20  is  the  theoretical  C*  multiplied  by  an  empirical  C*  efficiency 
(Section  4.1)  and  an  insulation  dilution  factor  (Section  S.3). 

The  thermochemical  properties  required  for  the  ballistics  calculation  ere 
obtained  horn  the  theoretical  Performance  Module,  the  nozzle  erosion  rate  re- 
quired to  define  the  instantaneous  throat  area  may  be  input  as  a table,  or  calcu- 
lated internally  based  on  the  simple  correlations  presented  in  Section  S.2. 

3. 3. 2. 3 Other  Mass  Contributions 

Mass  may  enter  the  rocket  motor  from  a gas  generator,  an  igniter,  or  in- 
sulation oblation.  These  are  treated  as  lumped,  rather  than  distributed,  additions; 
hence,  they  influence  the  pressure  magnitude  but  not  the  shape  of  the  pressure 
distribution  down  the  grain  port.  For  lumped  additions,  it  does  not  matter  whether 
the  mass  enters  at  the  head  end  or  the  aft  end  of  the  motor.  In  constructing  the 
Iterations,  tt  was  entered  at  the  aft  end  (which  best  represents  insulation  in  general), 
therefore. 


W(NN)  = f me  dx  + WI  + WIGN  (3-21) 

or 

W(NN)  = W(N)  + n Aj  fj  + WIGN  (3-22) 


W(N)  = 
Bl 

*1  “ 
rI 

WIGN  = 


propellant  flow  rate 
insulation  density 
Insulation  exposed  area 
insulation  ablation  rate 
igniter  or  gas  generator  flow 
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These  masses  are  Input  in  tabular  form.  IgrJter  {or  gas  generated  flow  la  input 
as  a function  of  lime,  Insulation  density  and  burning  rate  are  input  constants, 
while  exposed  in*  ul  at  too  area  is  input  as  a function  cf  time.  Ail  of  these  masses 
are  assessed  to  have  tire  same  thenaochemicaj  effects.  All  are  optional  end  not 
required  input.  It  should  not  be  interred  from  the  Inclusion  of  an  igniter  mass  flow 
that  the  analysis  can  predict  ignition  delay:  it  does  not  do  so.  The  analysis 
ascuaes  that  burning  has  started  to  a motor  initially  at  ambient  pressure,  and  will 
property  fill  the  chamber  {also  empty  the  chamber  at  burnout).  Therefore,  WIGN 
is  an  artifice  and  except  where  ft  Is  a significant  contribution  to  the  total  Impulse 
cf  the  motor,  can  be  included  or  excluded  as  desired.  If  the  user  Is  seriously 
interested  to  predicting  the  ignition  transient,  the  method  of  Reference  (13) 
should  be  employed. 


3.3 ,2.4  Flow  Across  Radial  Slots  (Segmented  Motor) 


Pressure  drop  across  a radial  slot  separating  the  segments  of  a segmented 
motor  is  treated  for  two  cases:  (1)  where  the  gas  velocity  in  the  slot  exceeds  the 
gas  velocity  in  the  adjacent  upstream  port:  U)  where  it  is  less  than  that  velocity. 
In  either  case,  the  state  properties  to  the  slot  are  assumed  equal  to  the  static 
values  cf  the  gas  at  the  end  of  the  upstream  port;  slot  velocity  is  then  determined 
from  continuity  in  the  slot: 


u 


s 


w 


ADD 


(3-23) 


«s  = slot  gas  velocity,  to. /sec 

WADD  “ A0*  rate  or  propellant  from  the  slot  faces,  lb/sec 
As  = slot  width  flow  area,  in? 
p =•  gas  density,  lb/to? 

It  is  implied  that  pressure  Is  constant  in  tho  slot.  Further,  it  is  asserted  that 
burning  rates  are  equal  over  the  entire  slot  face  and  at  the  strand  rate  value. 
Propellant  to  the  slot  does  not  see  radiating  walls,  and  to  practice  the  separation 
will  be  large  enough  to  preclude  erosive  burning. 

If  the  entering  gas  velocity  from  die  slot  exceeds  the  approaching  gas 
velocity  in  the  port,  the  slot  pressure  drop  Is  calculated  as  one  Increment  with  the 
momentum  equation.  It  is  as  though  the  slot  were  physically  absent,  being  filled 
with  propellant  of  high  ma»  Thus  *>i  ra'  lal  momentum  Is  dissipated,  and  a 
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psssinu re  cbop  equivalent  to  accelerating  the  added  flow  to  port  velocity  Is  de- 
termined. 

If  the  entering  gee  velocity  from  the  slot  is  less  than  the  approaching 
gas  velocity  In  the  port,  an  empirical  relation  Is  used  to  find  the  port  area  of 
the  expanded  stream  tube  over  the  slot: 

A , = A + (D*/A  ) (1  + -f>  (3-24) 

pi  pm  x pm  u 


v 

*pm 


stream  tube  area 

local  port  area 

slot  separation  distance 


Knowing  the  area  and  static  pressure  In  the  stxeamtube,  other  parameters  are  de- 
termined by  Iteration.  The  pressure  drop  across  this  streamtube  Into  the  down- 
stream port  Is  determined  with  the  assumption  of  an  Inlet  efficiency: 


rl+l 


[(pu*),  - {6u»)l+1|  (3-2S) 


P1+j  = downstream  static  pressure 
Pj  = upstream  static  pressure 

ETAIN  = Inlet  recovery  efficiency,  typically  0.9  for  shaip-edged  entrance 
geometry 


Optional  capability  Is  provided  to  refine  the  treatment  to  simulate  a slot 
face  which  Is  not  perpendicular  to  the  motor  axis.  Causes  of  a non-normal  face 
Include  variable  Ignition  delays  In  the  slot,  grain  slump,  or  tapered  Inhibitors. 
The  corrections  are  used  to  locate  the  axial  position  of  the  outer  slot  diameter, 
knowing  the  position  of  the  Inner  slot  diameter  at  the  grain  port.  Inputs  relating 
to  these  effects  are  tables  of  the  amount  of  burning  slot  diameter  versus  time, 
and  lateral  burn-back  versus  web. 


The  analysis  of  the  slot  ballistics  Is  coupled  to  the  analysis  of  the  slot 
geometry,  and  overrides  the  normal  grain  geometry  calculations.  A card  input  keys 
the  existence  of  a slot. 


3. 3. 2.5  Delivered  Performance 

The  various  quantities  which  characterize  delivered  performance  are  also 
calculated  with  the  Grain  Design  and  Ballistics  Module.  The  performance  calcu- 
lations within  this  module  do  not  replace  the  ones  previously  described  (Sections 
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2.3  and  3.1).  They  are  supplementary;  to  be  referred  to  only  when  the  Grain 
Design  and  Ballistics  Module  Is  used  in  a stand-alone  mode. 

Delivered  thrust  (vacuum)  Is  given  by 

V 

Fd  = W(N)  -Z-  „F  + Pe  Ae  (3-26) 


W(N)  = propellant  flow  rate,  Ib/sec 

Ve  * exit  velocity,  ft/sec 

fjp  » input  nozzle  efficiency 

Pe  = exit  static  pressure,  lb/in? 

= nozzle  exit  area,  in? 

Exit  velocity  and  pressure  are  determined  from  the  standard  Uentmplc  re- 
lations. The  nozzle  efficiency  may  be  determined  from  the  seml-emplrical  corre- 
lations, which  are  amenable  to  hand  calculation.  The  exact  calculations  (TD2P, 
etc.)  would  not  be  available  in  the  stand-alone  mode.  As  used  here,  nozzle  effi- 
ciency Is  an  input  value.  Thrust  delivered  to  the  Input,  ambient  pressure,  P&,  Is 
given  by 


Ambient  pressure  may  be  a constant  or  a table  to  simulate  a flight  schedule. 

Total  impulse  Is  calculated  as  the  time  integral  of  Fp  . Specific  Impulse 

a 

Is  calculated  on  tn  Instantaneous  basis  by  dividing  Fp  by  W(H),  and  on  a total 
basis  by  dividing  the  Integral  of  Fp  by  the  Integral  of  W(N).  Insulation  weight 
Is  not  .ncluded  because  of  the  convention  to  base  performance  on  propellant  weight 
only;  however,  both  Instantaneous  and  Integrated  Insulation  weight  flow  ere  output 
If  such  correction  is  desired. 


The  loss  of  performance  tn  solid  rocket  motors  due  to  finite  rate  chemical 
kinetics  has  been  estimated  to  range  between  0.2%  for  high  pressure,  low  metal 
loading  systems,  to  as  high  as  5%,  or  mere,  Ic:  low  pressure,  high  metal  loading 
systems,  Bence,  if  the  performance  of  solid  propellant  rocket  motors  Is  to  be 
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accurately  predicted  over  a wide  range  of  configurations  and  operating  conditions, 
the  loss  due  to  chemical  kinetics  must  be  accurately  modeled. 

The  loss  of  performance  which  we  attribute  to  finite  rate  chemical  kinetics 
Is  usually  associated  with  the  difference  between  the  maximum  possible  sensible 
heat  re,ease  (chemical  equilibrium)  and  the  actual  heat  release  due  to  chemical 
reactions  which  occur  In  the  rocket  nozzle.  The  loss  Is  small  If  the  reactions 
proceed  at  a rate  comparable  to  a small  fraction  of  the  stay  time  In  nozzle.  How- 
ever, If  the  recombination  reaction  rates  are  slow  all  of  the  heat  released  by  com- 
bustion of  the  propellant  cannot  be  converted  Into  kinetic  energy.  As  the  ratio  of 
nozzle  stay  time  to  a characteristic  reaction  time  approaches  zero,  the  so-called 
frozen  performance  limit  Is  approached. 

The  selection  of  a method  to  compute  the  kinetics  loss  was  not  as  straight- 
forward as  some  of  the  other  selection  procedures.  As  pointed  out  !n  Reference  1, 
there  are  two  basic  techniques  which  have  been  used  for  calculating  finite  rate 
performance  losses.  The  first  method  seeks  numerical  solutions  to  the  exact 
one  dimensional,  two  phase,  reacting  gas,  equations  of  motion.  The  second 
technique  is  approximate,  and  Is  based  upon  application  of  a “sudden  freeze" 
criteria. 

Equilibrium  flow  is  assumed  up  a given  point  In  the  nozzle,  at  which  the 
flow  Is  then  frozen.  The  sudden  freeze  approximation  Is  more  applicable  to  liquid 
rockets  than  solid  rocket  motors.  Such  methods  were  not  considered  for  the  pre- 
sent program  as  their  accuracy  was  not  deemed  to  be  high  enough.  The  only  ex- 

(19) 

tstlng  progr.-n  of  the  first  type  Is  that  of  Kllegel,  et  al  .As  discussed  below, 
this  program  was  also  considered  Inappropriate  for  the  present  purpose. 

The  one  dimensional,  two  phase,  reacting  gas  program  of  Reference  19 
yields  accurate  one  dimensional  performance  predictions,  but  does  not  directly 
yield  an  estimate  for  the  kinetics  loss  alone.  The  two  phase  flow  loss  Is  presently 
calculated  on  a two  dimensional  basis  using  the  TD2P  program  (Section  3.5). 

Thus,  If  the  one  dimensional,  two  phase,  reacting  gas  program  Is  to  be  used,  the 
two  phase  and  kinetics  losses,  which  are  coupled  in  this  analysis,  must  be  separ- 
ated in  order  to  avoid  a redundancy. 

An  accurate  method  for  separating  the  two  losses  could  be  developed  by 
modifying  the  one  dimensional  equilibrium  program  (ODE)  to  handle  finite  particle 
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lags.  Such  an  approach  was  discarded,  however,  due  to  the  extent  of  the  modi- 
fications required.  Other,  more  approximate  means  of  obtaining  a kinetics  loss 
could  be  developed,  but  not  without  compromising  the  resultant  accuracy.  The 
program  described  In  Reference  19  Is  also  complex,  very  large,  long  running, 
and  not  easy  to  use.  For  all  of  these  reasons  alternate  means  were  considered. 

Since  the  present  effort  was  not  to  Include  any  new  computer  program  de- 
velopment, the  JANN/iF  standard  computer  program  for  calculating  one  dimensional 
kinetics  losses  In  liquid  rockets,  ODK,  was  considered.  The  ODK  program  Is 
widely  used  and  Is  computationally  efficient.  Complete  documentation  of  the  analysis 
and  computer  program  are  contained  In  Reference  5.  In  Its  existing  form  ODK 
was  not  applicable,  however,  since  It  contained  no  provision  for  solid  or  liquid 
phases. 

The  existence  of  a second  phase  In  the  flow  can  result  In  significant  al- 
teration of  the  axial  distributions  of  the  various  flow  quantities  Cl,  B,  u,  etc.), 
and,  hence,  impacts  the  reaction  rates  and  Kinetic  losses.  The  existence  of 
finite  lags  between  the  particles  and  gas  modifies  the  effect  of  two  phase  flow 
on  the  kinetics  loss,  but  In  a second  order  manner. 

Various  modifications  could  be  made  to  the  ODK  program  to  Include  the 
effect  of  two  phase  flow  on  the  kinetics  loss.  In  Its  unmodified  form  the  ODK 
program  numerically  Integrates  the  exact  one  dimensional,  finite  rate  kinetics, 
equations  of  motion  for  a gas  only  system.  A modification  of  the  program  to 
allow  It  to  treat  the  particles  <n  an  exact  manner  would  have  resulted  In  a program 
similar  -o  that  of  Reference  19,  and  would  have  engendered  the  same  problems. 

A modification  based  on  the  adoption  of  the  constant  fractional  lag  approximation 
was  also  considered.  In  this  approximation  the  velocity  and  temperature  lags 
of  the  particles  are  assumed  to  be  finite,  but  constant.  With  this  assumption  the 
one  dimensional,  two  phase,  equations  of  motion  can  be  reduced  to  those  for  an 
"equivalent"  gas  only  system,  thereby  eliminating  the  need  for  solving  a separate 
set  of  particle  equations.  The  required  modifications  to  the  ODK  program,  and 
the  ODE  program  (in  order  to  generate  a reference  equilibrium  solution),  were  de- 
termined to  be  substantially  less  than  those  required  to  Implement  a solution  of 
tho  “exact"  two  phase  flow  equations.  The  effort  required  was,  however,  still 
greater  than  desired. 
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Another  possible  modification,  actually  a degenerate  form  of  the  previous 
one,  was  also  evaluated.  In  the  limit  of  vanishingly  small  lags,  the  constant 
fractional  lag  equations  reduce  to  the  so  called  equilibrium  limit,  l.e.  particles 
and  gas  have  the  same  velocity  and  temperature.  The  system  of  equations  again 
reduces  to  that  for  an  "equivalent"  gas  only  flow,  but  the  program  modifications 
required  to  Implement  a no  lag  model  wen  simpler  and  less  extensive  than  those 
for  the  other  alternatives  examined.  While  this  approach  seemed  attractive,  a 
series  of  calculations  was  performed  In  order  to  estimate  the  effect  of  Ignoring  finite 
lags  before  It  was  adopted.  As  a result  of  these  calculations  It  was  determined 
that  die  zero  lag  assumption  Is  quite  adequate.  It  Is  esttmeted  that  the  use 
of  die  zero  lag  approximation  should  impact  the  overall  performance  prediction 
by  less  than  0.1*  (when  the  kinetics  loss  itself  Is  but  a fraction  of  a percent 
the  Impact  should  be  closer  to  0.01%). 

In  view  of  the  previous  discussion, the  ODK  program  modified  to  treat  par- 
ticles on  an  equilibrium,  no  lag,  basis,  was  adopted  as  the  method  to  be  used  in 
computing  the  kinetics  loss.  A kinetic  efficiency  factor,  njau'  ls  deflned  by 
taking  die  ratio  of  the  Isp  (vacuum)  calculated  by  ODK  to  an  equilibrium  Igp  ob- 
tained by  exercising  a special  option  which  was  incorporated  in  the  ODE  program. 
When  approached  In  this  manner  the  kinetics  loss  is  obtained  directly;  without 
the  need  to  separate  out  redundant  two  phase  flow  losses.  The  only  requirement 
ls  that  the  particles  be  treated  Identically  In  ODK  and  ODE, 

Both  ODE  and  the  modified  ODK  program  assume  the  particles  and  gas  are 
In  thermal  and  dynamic  equilibrium.  The  ODE  program  In  Its  original  form,  allows 
for  the  gas  phase  to  condense  as  It  expands  out  the  nozzle,  l.e.  the  weight 
fraction  of  the  condensed  phase  may  Increase.  The  modified  ODK  program,  how- 
ever, does  not  allow  mass  transfer  between  phases.  : is,  to  be  consistent,  and 
to  preserve  accuracy,  a special  option  was  added  to  the  ODE  program  which  allows 
It  to  "freeze"  the  concentration  of  the  condensed  phases  at  their  chamber  values. 
The  same  dictates  of  consistency  and  accuracy  governed  the  treatment  of  particle 
solidification.  The  ODE  program  allows  the  liquid  phase  to  solidify,  and  properly 
accounts  for  the  resultant  energy  release.  An  analogous  treatment  of  particle  so- 
lidification was,  therefore,  Incorporated  Into  the  ODK  program. 

The  previously  mentioned  definition  of  the  kinetic  performance  efficiency, 
as  calculated,  may  be  written  In  the  following  symbolic  form. 
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wherein  1 denotes  the  specific  Impulse  calculated  using  the  so-called  restricted 

spR£ 

equilibrium  option  to  the  ODE  program.  An  additional  advantage  of  using  this  refer- 
ence method  of  calculating  if.  Is  that  only  thermodynamically  important  species 
sp 

need  be  considered  in  the  ODK  and  ODE-RE  calculations.  Thus,  trace  species, 
and  reactions  Involving  them,  need  not  be  Included  in  the  calculation. 

The  following  is  a brief  summary  of  the  ODK  program  (prior  to  modification). 
Further  details  regarding  the  details  of  the  analysis  and  the  program  Itself  may  be 
found  In  Volume  II  of  tills  report  and  In  Reference  * The  ODK  program  calculates 
the  invlscid,  one-dimensional,  nonequilibrium,  frozen,  and  equilibrium  (the  ODE 
program  Is  Incorporated  Into  ODK),  nozzle  expansion  of  gaseous  combustion  pro- 
ducts. The  program  can  treat  large,  complex,  chemical  systems.  A maximum 
of  ISO  distinct  chemical  reactions  and  40  Individual  species  can  be  included  In 
a given  calculation.  Chemical  reactions  are  input  in  stanaard  reaction  form  and 
are  translated  by  an  Input  processor  Into  a mathematical  form  suitable  for  compu- 
tation. Chemical  reactions  can  be  added  or  deleted  by  adding  (or  deleting)  a 
single  input  card  per  reaction.  Up  to  10  reactants  and  10  products  can  be  Included 
In  each  reaction.  Thera  Is  also  a convenient  method  for  specifying  third  body 
efficiencies.  The  coupled  set  of  nonlinear,  ordinary  differential  equations.  Is 
integrated  using  a very  efficient,  stable,  second  order.  Implicit  method.  Under 
certain  conditions,  when  efficient  to  do  so,  the  program  switches  to  an  explicit 
Integration  method.  The  program  Is  designed  for  englneeiing  use.  Is  well  docu- 
mented and  user  oriented. 

The  modifications  to  the  ODK  program  accomplished  during  the  present 
effort  are  described  In  the  next  section  and  in  Volume  II  of  this  report.  The  reaction 
set  and  reaction  rates,  for  aluminized  propellants,  that  were  employed  In  the  cal- 
culations performed  to  date  are  described  In  Appendix  A.  The  manner  In  which 
these  reactions  and  rates  were  selected  is  also  discussed  therein. 


3.4.2  Modifications  To  The  ODK  Program 

Modifications  to  the  ODK  program  were  required  In  order  to  Incorporate 
Into  It  the  ability  to  treat  two  phasr  —nansions,  with  solidification.  As  discussed 
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In  the  previous  section,  particle-gas  equilibrium  (zero  lag)  was  assumed.  Thus, 
the  gas  and  particle  velocities  and  temperatures  are  taken  to  be  equal.  Other 
assumptions  which  were  employed  are: 

• There  Is  no  phase  change  from  the  gaseous  phase  to  either  the  liquid  or 
solid  state. 

•The  liquid  particles  of  a given  specie  solidify  at  the  melting  point  of  the 
solid  phase. 

•The  particles  occupy  zero  volume  and  exert  no  pressure  on  the  gas. 

As  would  be  expected  for  a program  like  ODK,  a significant  effort  was  re- 
quired in  order  to  identify,  program  and  check  out  the  requisite  modifications.  As 
mentioned  previously,  the  zero  lag  assumption  allows  the  two  phase  flow  to  be 
treated  as  an  "equivalent"  gas  only  flow.  The  existing  equations  employed  In 
the  ODK  program  can  be  made  to  correspond  to  the  equations  for  an  equilibrium 
(no  lag)  partlcle-gcs  mixture  simply  by  treating  the  condensed  phases  „s  gas  phase 
species  with  Infinite  molecular  weights.  Thus,  the  only  the  equation  for  the  gas 
constant  required  special  treatment.  The  gas  constant  Is  defined  by 


where 


R = 2 CjRj 


R,  - — for  gas  phase  species 
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= 0 for  condensed  phase  species 


A method  for  handling  particle  solidification  also  had  to  be  Incorporated 
Into  the  program.  Under  the  zero  lag  assumption,  particle  solidification  must 
be  handled  in  a special  manner,  since  there  Is  no  explicit  heat  transfer  term  to 
govern  the  rate  at  which  a particle  solidifies.  It  is  assumed  that  when  the  gas 
reaches  the  particle  solidification  temperature  die  gas  temperature  remains  con- 
stant, l.e.,  d T/dx  =»  o,  until  the  latent  heat  of  solidification  Is  given  up  by  the 
liquid  particles.  This  constraint  yields  an  equation  for  the  rats  of  change  of  the 
condensed  phase  mass  fraction,  as  follows.  By  noting  the  constancy  of  the  total 
enthalpy,  and  applying  the  restriction  dl/dx=0.  It  Is  easily  shown  that  the  follow- 
ing equation  must  be  satisfied: 
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where  Ng  Is  equal 
relations. 


to  the  number  of  gas  phase  species, 

dcs  dcz 

"ST  = " Hx 


By  using  the  following 

(3-31) 


h -h  = aH  (heat  of  fusion) , v3— 32) 

is  rn 


equation  (3-30)  may  be  solved  for  the  rate  of  change  of  the  liquid  phase  mass 
fraction,  yielding, 
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The  ODK  program  integrates  the  quasl-ono  dimensional  equations  of  motion 
either  on  a pressure  defined,  or  area  defined  basis.  Therefore,  before  equation 
(3-33)  could  be  implemented  the  velocity  derivative  had  to  be  replaced  by  either 
a pressure  or  area  derivative.  Using  the  equations  of  motion  given  in  Reference 
S,  it  can  be  shown  that  equation  (3-33)  can  be  written  in  the  following  alternate 
forms,  for  pressure,  or  area  defined  problems,  respectively. 


~k  - j’r1  T -K  - Bj  pressure  defined  (3-34) 

-S?  = <^npr-  | <7#  * A,'B|  Area  deflned  {3*3S) 


With  dc^dx  and  dc^/dx  defined,  the  remaining  equations  can  be  Integrated. 

This  procedure  is  then  applied  at  each  succeeding  Integration  step  until  the  concen- 
tration of  the  liquid  phase  becomes  equal  to  zero  (l.e.  all  of  the  heat  of  fusion 
contained  in  the  original  mass  fraction  of  liquid  has  been  given  up  to  the  gas  phase). 

Since  this  part  of  the  analysis  Is  new,  and  cannot  be  found  elsewhere,  the 
partial  derivatives  of  the  rate  of  change  of  the  liquid  mass  fraction,  dc^dx  (de- 
noted c ),  viith  respect  to  all  of  the  other  dependent  varlalbes  are  given  here,  for 
complet<5hes:t.  These  partial  derivatives  are  required  by  the  Implicit  integration 
scheme  used  In  ODK. 
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Pressure  Defined  Equations 
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Area  Defined  Equations 
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where  Tx  = dT/dx  = gas  temperature  derivative  In  the  absence  of  solidlilcatlon. 

the  0DK  computer  code  has  been  dimensioned  to  handle  up  to  10  condensed 
phase  species  (or  5 doublets  of  solid  and  liquid  phase)  with  10  discrete  melting 
points.  This  allotment  should  be  adequate  to  handle  any  metalized  propellant 
system  of  Interest. 


3-28 


While  strictly  a modification  to  ODE,  and  not  DDK,  the  manner  In  which 
ODE  was  modified  to  restrict  the  condensed  phase  mass  fract.on  from  changing 
Is  described  here,  as  this  “restricted  equilibrium"  modification  serves  only 
to  provide  a reference  I value  for  the  kinetics  loss  computation. 

For  purposes  of  computing  a reference  (or  base)  I for  the  kinetics  loss 

sp 

calculation  an  option  was  Incorporated  Into  the  ODE  program  which  allows  it  to 
"freeze"  the  concentrations  of  the  condensed  phase  species  at  their  chamber  values 
(except  that  corresponding  liquid  and  solid  phases  are  allowed  to  exchange,  l.e., 
solidification  Is  allowed). 

In  order  to  Implement  this  option,  the  ODE  program  searches  the  list  of 
species  to  see  If  aluminum  (Al),  beryllium  (Be),  boron  (B),  lithium  (LI) , magnesium 
(Mg),  zinc  (Zn),  zirconium  *Zr)  , or  any  compounds  containing  these  elements,  are 
present  In  either  the  liquid  or  solid  phases.  It  any  such  species  are  found,  the 
corresponding  liquid  and  solid  pnases  aie  "enamed,  but  assigned  their  usual 
thermodynamic  properties.  Tats  is  done  at,  that  the  program  will  not  recognize 
them  fi'wft-t  they  are,  tnereby  restricting  these  liquid-solid  pairs  from  reacting 
with  any  species  other  than  themselves.  As  far  as  the  gaseous  species  are  concer- 
ned, these  condensed  phases  no  longer  exist,  hence,  uddltonal  amounts  of  conden- 
sed phase  cannot  be  formed  as  the  flow  expands  out  the  nozzle. 

3,5  Two  Dimensional  Two  Phase  Module 

3.5.1  '-nproach 

Almost  all  rocket  nozzles  of  practical  interest  have  a Joss  tr,  perfonnance 
due  to  a nor.axlal  component  of  velocity  at  the  nozzle  exit  plane.  This  loss  is 
usually  referred  to  as  a divergence,  or  two  dimensional  flow,  loss.  If  the  pro- 
pellant Is  metalized,  there  is  an  additional  loss  when  the  particles  cannot  main- 
tain a state  of  dynamic  and  thermal  equilibrium  with  the  gas  flow.  This  loss  Is 
referred  to  as  the  two  phase  flow  loss. 

The  two  dimensional  and  two  phase  tlow  effects  are  coupled.  From  the 
discussion  In  Section  2.4  It  Is  clour  :hat  any  analytical  methodology  for  predic- 
ting two  phase  flow  loss  should  be  Integrated  with  the  calculation  of  dl’Xigence 
lass.  This  requires  a two-dimensional  flow  Held  program  with  two  phase  capab.llty. 
Both  a survey  conducted  by  JAN'  . «-e  study  reported  In  Reference  1 reco-n- 
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mended,  as  the  best  available  computer  program  of  this  type,  the  Two  Dimensional  j 

Two  Phase  Flew  computer  program  developed  by  Kllegel  and  Nickerson^.  This 
program  was,  therefore,  adopted  for  use  In  the  present  performance  prediction  j 

methodology.  The  basic  assumptions  employed  In  the  Tti2P  program  are:  \ 

• The  flow  Is  at  the  choked  steady  state  condition  for  an  axially  symmetric  ; 

converging  - diverging  nozzle. 

SThe  gas  and  particle  phases  each  have  constant  thermodynamic  properties. 

•There  are  no  mass  or  energy  losses  from  the  system.  i 

«Ihc  gas  is  Inviscld  except  for  Its  Interactions  with  the  condensed  particles.  ] 

•The  volume  occupied  by  the  condensed  particles  Is  negligible.  j 

•The  thermal  (Brownian)  motion  of  the  condensed  particles  Is  negligible.  ^ 

•The  condensed  particles  do  not  Interact.  ; 

•the  condensed  particle  size  distribution  may  be  approximated  by  groups 
of  different  size  spheres.  i 

•The  internal  temperature  of  the  condensed  particles  Is  uniform. 

•Energy  exchange  between  the  gas  and  the  condensed  particles  occurs  only  - 

by  convection. 

•The  only  forces  on  the  condensed  particles  are  viscous  drag  forces. 

•There  Is  no  mass  transfer  from  the  gas  to  the  condensed  phase  during 
the  nozzle  expansion.  j 

6 

The  overall  method  of  solution  consists  of  essentially  two  distinct  parts:  i 

a transonic  sc’utlon;  and  a supersonic  solution.  The  transonic  solution  is  required 
in  order  to  provide  an  Initial  line  for  the  supersonic  solution. 

The  transonic  solution  I s obtained  by  an  approximate  technique.  The 
nozzle  geometry  Is  assumed  to  consist  of  a conical  inlet  section  Joined  smoothly 
to  a constant  radius  of  curvature  threat.  The  flow  in  the  conical  Inlet  portion  of 
the  nozzle  Is  assumed  lo  be  a one-dimensional  sink  flow.  The  flow  in  this  region 
is  coupled  to  a two  dimensional,  Sauer  type,  expansion  solution  in  the  throat  re- 
gion. The  output  from  this  analysis  Is  a set  of  flow  properties  along  a supersonic 
start  line.  The  flow  field  from  this  initial  line  out  to  the  nozzle  exit  is  obtained 
numerically  using  the  method  of  characteristics  to  sorve  the  exact  two  phase  flow 
equations  of  motion. 

This  overall  procediue  works  quite  well  In  most  cases,  however,  the  naiur® 
of  the  transonic  solution  Introduces  certain  restrictions  on  the  problems  which 
can  be  solved.  As  previously  mentioned  the  transonic  analysis  requires  that  the 


3-30 

i 

i 


l 


convergent  portion  of  the  mzzle  be  conical  and  that  the  throat  be  circular,  with 
a tingle  radius  of  curvature.  Nozzles  not  of  this  shape  must  be  “fitted"  to  conform 
to  this  geometry  in  order  to  obtain  a solution.  While  this  can  often  be  accomplished 
without  Introducing  significant  error,  the  "fitting"  procedure  Is  subjective,  which 
is  undesirable,  and  the  size  of  the  error  Introduced  Is  net  known  a priori.  The 
transonic  method  of  solution  will  also  fail  when  the  nondinenslonal  throat  radius 
of  curvature  (radius  of  curvature  divided  by  throat  radius)  Is  less  than  about  l.S. 
iind  when  the  conical  Inlet  angle  exceeds  approximately  45s.  This  fact  results  in 
two  additional  restrictions  on  the  nozzle  geometries  which  may  be  considered. 

It  Is  clear  from  this  discussion  that  a mare  exact,  less  restrictive,  method 
for  solving  the  two-dimensional . two-phase  flow  equations  would  be  desirable. 

With  presently  available  numerical  techniques  such  a method  of  solution  could, 
and  should,  be  developed. 

There  is  one  additional  limitation  of  the  TD2P  program  that  has  yet  to  be 
mentioned.  If  the  exhaust  nozzle  Is  contoured  In  such  a manner  that  the  particles 
impinge  upon  It  before  reaching  the  exit  plane,  the  supersonic  solution  will  fail 
terminate)  at  the  first  point  of  Impingement.  Xt  is  a relatively  simple  matter  to 
modify  the  program  to  continue  the  calculation  In  an  artificial  manner  after  Impinge- 
ment occurs.  It  is  not  quite  so  simple,  but  not  too  difficult,  to  Incorporate  a 
simple,  but  reasonable,  physical  model  of  Impingement  into  the  program. 

In  order  to  use  the  TD2P  program  It  Is  necessary  to  determine  appropriate 
average  values  for  the  perfect  gas  and  particle  phase  properties.  A method  for 
obtaining  such  properties  was  developed  by  XI  leg  el  and  Nlcherson  and  was  for- 
malized In  the  Average  Gas  Properties  (ACT)  computer  program This  program 
uses  the  one -dimensional  equilibrium  (zero  lag)  approximation  to  find  a set  of 
"equivalent"  perfect  gas  properties  such  that  the  solidification  area  ratio  Is  Iden- 
tical to  the  value  computed  by  the  ODE  program.  There  are  other  details  involved 
In  this  computation  which  will  not  be  discussed  at  this  point.  An  updated,  modi- 
fied, version  of  the  ACT  program  has  been  Incorporated  In  the  TD2P  program,  and 
Is  fully  described  in  Volume  □ of  this  report  (see  subroutine  ACT). 

Because  the  two-phase  flow  loss  Is  an  Important  source  cf  performance 
loss  and  Is  a function  of  particle  size,  the  ability  to  predict  particle  size  will 
have  an  important  Impact  upon  the  success  of  the  computer  program  over  a broad 
range  of  conditions.  The  present  method  for  predicting  particle  size  Is  described 
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in  detail  in  Section  5.4.  The  correlation  described  therein  yields  the  mean  particle 
size.  The  distribution  function  employed,  and  the  method  for  selecting  the  dia- 
meter at  each  of  the  particle  groups  is  described  in  the  subroutine  DIST  write  up 
in  Volume  E. 
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3.5.2  Modifications  tr>  the  TD2P  Program 

Numerous  modifications  and  additions  to  the  TD2P  program  were  accomplish- 
ed during  the  present  effort:  The  more  significant  of  these  changes  are  discussed 
below. 

The  TD2P  program  was  originally  written  In  FORTRAN  U,  and  the  common 
linkages  were  completely  unsatisfactory  for  use  as  a module  element.  The  pro- 
gram had  only  blank  common  with  multiply  defined  entries.  These  were  eliminated 
and  data  storage  replaced  by  multiple  labeled  common  statements.  These  modi- 
fications were  necessary  to  reduce  program  size  to  the  55X  requirement  for  the  com- 
bined program.  The  TD2P  over  lay  structure  was  revised  accordingly. 

The  number  of  options  for  describing  the  nozzle  wall  contour  was  increased. 

In  particular,  the  ability  to  use  a cubic  spline  fit  contour,  rather  than  a point  by 
point  specification,  was  added.  A cone,  defined  by  the  nozzle  end  point  only, 
option  was  also  incorporated  into  the  program.  The  maximum  number  of  wall  points 
allowed  was  also  increased  from  60  to  150  to  allow  better  definition  of  the  contour. 

The  TD2F  input  was  revised  to  be  more  convenient  for  the  user.  Certain 
input  items  were  built  into  the  program,  while  others  were  modified  or  eliminated. 

The  partible  size  correlation  presented  In  Section  5.4  was  incorporated 
into  the  program,  as  was  a subroutine  CDISfl  to  provide  a log-normal  distribution 
into  the  specified  number  of  particle  groups. 

(23) 

The  AGP  program  was  modified  and  then  Incorporated  into  the  TD2P  program. 
The  AGP  program  provide!;  a rational  means  for  selecting  perfect  gas  properties 
to  use  in  the  TD2P  calculation.  Previously,  this  program  had  to  be  executed 
separately  and  the  requisite  parameters  input  to  the  TD2P  program.  The  AGP  pro- 
gram was  modified  to  calculate  I on  a cne-dlmenslonal,  zero  lag,  basis,  to 

sp 

serve  as  a reference  value  for  the  two-dimensional,  two-phase  loss  calculation. 

The  manner  in  which  the  AGP  program  calculates  specific  heats,  and  the  way  it 
treats  solidification,  was  also  modified. 


3-32 


[ 


I 

i 

i 


j 

The  TD2F  prog  ram  was  also  modified  to  provide  certain  quantities  to  the  { 

roundary  Layer  Lons  Moduli.  These  quantities  are  the  internally  tabulated  wall 
geometry,  the  Mach  number  distribution  along  the  wall  streamline,  the  two-diaen-  ij 

siosa.,  two-pnase,  mass  flow  rate,  and  the  total  pressure  and  total  temperature  ' 

at  the  wall  initial  11ns  point.  The  latter  two  quantities  are  provided  to  allow  the  ]j 

TBL  calculated  wall  pressure  and  velocity  distributions  to  closely  approximate  the 
corresponding  TD2P  value: . 

The  object  of  the  HJ2P  Module  Is,  of  course,  to  calculate  a performance 
loss  efficiency.  The  manner  In  which  the  program  war  modified  to  do  this  Is  ] 

described  In  the  next  sub-section.  \ 

| 

3.5.3  Performance  Loss  Calculation 

The  TD2P  Module  provides  throe  types  of  quantities  to  the  Master  Control  ‘ 

Module  for  use  in  the  overall  delivered  performance  calculation.  These  quantities 
are  the  nozzle  discharge  coefficient,  CD:  and  two-dimensional,  two-phase,  specific 

Impulse,  I , and  specific  Impulse  efficiency,  nlrwo.  as  functions  of  area 

ratio.  PTD2P  ^D2P 

The  discharge  coefficient,  C^,  is  obtained  as  a direct  result  of  the  TD2P  ’ 

transonic  solution.  It  is  used  In  obtaining  the  combustion  efficiency  from  the  c*  J 

efficiency  - equation  (3-4),  and  in  computing  the  ambient  pressure  correction  to 

delivered  I - equation  (2-4).  £ 

sp  v 

The  specific  impulse  calculated  by  the  TD2P  program  I , Is  used  in 

spTD2P 

the  performance  calculation,  but  never  on  a absolute  basis.  At  every  wall  point 

in  the  TD2P  characteristics  solution  a two-dimensional,  two-phase  flow,  I 

efficiency,  tlyQjpif)  1*  calculated  by  ratiolng  the  TD2P  I to  a cme-dtoenslonal. 

zero  lag,  1 computed  at  the  same  area  ratio.  This  one-dlmenslonal,  reference 

Ijp,  Is  obtained  from  the  subroutine  used  by  the  modified  AGP  program  (AGP  was  | 

Incorporated  into  the  TD2P  program)  using  the  same  Ideal  gas  properties  as  employed  ) 

In  die  TD2P  program.  By  using  a ratio  calculated  In  this  manner.  Instead  of  using  5 

the  absolute  IJp  number  from  TD2P,  the  effect  of  neglecting  real  gas  effects  Is  j 

minimized.  Also,  extra  cere  was  taken  to  accurately  calculate  the  reference  I 1 

sp  J 

during  solidification  so  as  to  avoid  the  use  of  different  assumptions  in  the  two  | 

calculations.  The  equations  used  in  calculating  the  one  dimensional  reference  I j 

sp  3 

are  given  below.  ] 
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The  following  equations  were  iterated  using  the  method  of  secants  to  calcu- 
late the  ca»  dimensional  reference  I for  a given  expansion  ratio,  , . The  three 
sen  of  equations  reflect  the  regions  of  interest.  These  are  before  solidification 
of  the  condensed  phase  begins,  during  solidification,  and  after  solidification  ends. 
In  the  first  and  last  regions,  temperature  was  used  as  the  Independent  variable, 
while  during  solidification,  the  ratio  of  the  amount  of  solidified  condensed  phase 
to  total  condensed  phase  was  used. 


Before  solidification  begins;  t*. 


where. 
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where. 


a«„ 


“p  “g  RT 
m 

gas  velocity  at  beginning  of  solidification 

expansion  ratio  at  beginning  of  solidification 
weight  flow  of  solid  condensed  phase 
weight  Cow  of  the  gas 

weight  flow  of  the  total  condensed  phase  Clqutd  plus  solid) 
melt  temperature  of  the  solid  phase 


3-34 


I 


* 'end  of  solidification 


U2 

. 


u = gas  velocity  at  the  end  of  solidification 
3es 

Bs  - (1  + ip/tigl/U  + ^ -JW) 

9 p.9 

In  all  regions  the  one-dimensional  zero  lag  vacuum  I , I , is  cal- 

sp  spIDOL 

ciliated  as  follows: 
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The  efficiency,  *»tx>2P'  6191x341165  only  the  combined  two-dimensional  two- 
phase  flow  effect  on  performance,  as  the  TD2?  and  reference  I calculations  are 
Identical  in  all  other  respects.  When  the  nozzle  throat  erodes,  nozzle  exit  area  ratio, 
and,  hence,  tjjpjP'  ls  tlm6  dependent.  The  program  calculates  a time  averaged 
2-D,  2 phase,  efficiency  In  the  manner  shown  In  equation  (3-3). 


The  specific  impulse  calculated  by  the  TD2P  program  (as  a function  of 
area  ratio)  is  also  used  to  compute  the  erosion  loss.  While  this  loss  could  have 
been  estimated  based  on  one-dimensional  results  (using  the  ODE  program).  It  was 
felt  that  a two-dimensional  calculation  would  be  superior.  As  defined  herein,  the 
erosion  loss  Is  equal  to  the  ratio  of  the  2-D,  2-phase  I at  the  time  averaged 

sp 

expansion  ratio,  to  the  specific  Impulse  calculated  at  the  initial  expansion  ratio. 


i.e.  i 


I 

^erosion  = 1 


spTD2P 

spTD2P 
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It  should  be  noted  that,  here  again,  I values  are  used  on  a relative, 

rather  than  an  absolute,  basis.  This  effi8?£icy  ls  not  reported  separately,  but  Is 
combined  with  the  time  averaged  two-dimensional  two  phase  efficiency,  »7td2P' 
to  yield  a single  value,  ' S6e  equMioo  (3-2). 

In  the  present  analysis,  the  fact  that  the  two  phase  flow  loss  tends  to 
decrease  with  Increased  throat  diameter  (erosion)  was  Intentionally  Ignored,  since 
It  was  felt  that  the  added  loss  due  to  throat  roughness  would  compensate  for 
this  effect. 


3.6  Boundary  Laver  Loss  Module 
3.6.1  Approach 

IVro  dimensional  two  phase  flow  solutions  Indicate  that  the  bulk  of  the 
particulate  matter  gets  channeled  towards  the  centerline  of  the  motor  due  to  initial 
effects  In  the  throat  region  of  the  nozzle.  As  a result  of  this  phenomena  gas 
phase  only  boundary  layer  calculations  can  usually  be  made  without  fear  of  Intro- 
ducing significant  error  due  to  the  neglect  of  the  particulate  matter. 

The  usual  boundary  layer  loss  mechanisms  which  must  be  accounted  for  In 
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rocket  engine  perfonnance  predictions  are  the  energy  loss  through  the.  walls  due 
to  heat  transfer,  and  the  momentum  loss  due  to  sheer  stress  at  the  wall.  As  dis- 
cussed In  Reference  24,  these  losses  can  either  be  found  by  properly  integrating 
the  stress  tensor  along  the  wall  (pressure  plus  skin  friction)  or  by  calculating  the 
momentum  deficit  In  the  boundary  layer  at  the  nozzle  exit  plane.  The  latter  method 
Is  the  simplest  and  easiest  to  use,  and  Is  die  recommended  JANNAF  procedure 
Following  this  procedure  die  boundary  layer  loss  may  be  calculated  as 


iFtt-(2»rpU»e  cos  a)e 


1- 


j£_P 
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where  r,  p,  U and  P are  the  nozzle  radius,  density,  velocity  and  pressure  at  the 
wall  exit  station,  respectively,  6*  and  8 are  die  boundary  layer  displacement  and 
momentum  thicknesses,  and  a Is  the  angle  between  die  slope  of  the  nozzle  contour 
and  the  centerline. 


It  should  be  pointed  out  that  this  equation  gives  the  complete  boundary 
layer  loss.  The  reduction  In  fluid  momentum  due  to  both  shear  and  heat  transfer 
losses  to  the  wall  is  properly  accounted  for. 

The  boundary  layer  edge  conditions  that  ero  required  In  this  formula  are 
supplied  by  the  TD2P  Program.  This  allows  the  effect  of  two  dimensional  flow  In 
the  nozzle  to  be  Incorporated  Into  the  boundary  layer  solution.  This  Is  the  same 
procedure  as  adopted  by  JANNAF  In  the  liquid  propellant  perfonnance  methodology 
(l.e.  the  TDK  - TBL  Interface). 


Ihe  boundary  layer  thicknesses  5*  and  8 may  be  calculated  In  several  ways, 
categorized  below  In  order  of  Increasing  accuracy: 


Simple  correlational  techniques 
Integral  methods 

Finite  difference  and  related  numerical  methods 


The  correlational  techniques  are  normally  based  on  semi-emplrlcal  cor- 
relations of  experimental  data  for  flat  plate-zero  pressure  gradient-flows.  In 
many  cases  the  relations  are  based  on  Incompressible  results  and  are  applied  to 
compressible,  variable  property  flows,  through  die  use  of  compressibility  trans- 
formations and  reference  temperatures  (enthalpies).  Normally,  one  of  several 
forms  of  Reynolds  analogy  Is  ap-"<!<5  *o  relate  skin  friction  to  heat  transfer.  As 
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a result  of  their  simplicity,  these  methods  have  proven  to  be  of  great  utility  in 
preliminary  design  calculations;  however,  final  engine  design  calculations  are 
normally  based  on  more  accurate  techniques. 

the  majority  of  rocket  engine  heat  transfer  and  boundary  layer  loss  cal- 
culations are  currently  carried  out  using  integral  methods  such  as  the  TBL  program 
described  In  Reference  7.  The  integral  methods  are  based  upon  solutions  of  various 
integral  forms  of  the  differential  equations  of  motion.  Ihese  methods  give  two 
ordinary  differential  equations  for  die  momentum  thickness,  g,  and  energy  thick- 
ness, j*,  (one  equation  for  g only,  when  the  energy  equation  Is  replaced  by  the 
Crocco  relation  relating  the  temperature  profiles  to  velocity,  as  Is  the  case  In 
TBL).  These  equations  cannot  be  solved  without  additional  assumptions  for  the 
form  or  shape  factor,  the  surface  shear  stress,  r , and  heat  transfer,  qw< 
and  the  boundary  conditions  such  as  wall  temperature  and  pressure  gradients. 

The  accuracy  of  Integral  methods  Is  limited  by  the  types  of  relations 
which  are  adopted  for  computing  T^,  qw,  etc.  The  majorityof  the  integral  methods 
that  have  been  applied  to  rocket  nozzle  flows  proceed  on  the  basis  of  all  cr  most 
of  the  following  assumptions: 

Tw  can  be  calculated  from  correlations  based  on  zero  pressure  gradient 

flows. 

qw  can  be  calculated  using  Reynolds  analogy. 

The  shape  factor,  H=«*/e,  can  be  determined  using  power  law  velocity  pro- 
files (usually  a l/7th  power  Is  taken). 

The  Crocco  relation  can  be  used  to  calculate  the  temperature  profiles. 

Reference  temperature  or  enthalpy  methods  such  as  Eckert's  can  be  used 

to  account  for  variable  properties  in  the  boundary  layer. 

All  of  these  assumptions  are  violated  to  various  degrees  in  variable  pro- 
perty flows  In  cool  wall  rocket  motors;  however.  Integral  methods  of  the  above 
type.  Including  TBL,  give  results  which  In  many  cases  compare  fairly  well  with 
experimental  data.  In  the  supersonic  portion  of  the  nozzle,  TBL  heat  transfer 

predictions  are  generally  within  2056  o:  the  measured  values,  while  in  the  general 

(26) 

region  of  the  throat,  discrepancies  of  50%  or  more  have  been  noticed'  The 
large  discrepancies  between  the  Integral  theories  of  TBL  type  and  the  convergent 
section-throat  data  can  be  attributed  to  the  presence  of  large  prossure  gradients 
In  the  throat  region  which  seriously  loiote  some  of  the  assumptions  upon  which 
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the  analysis  Is  based,  the  occurence  of  "lamlnarlsatton*  of  the  boundary  layer 
due  to  die  pressure  gradients,  and  die  sensitivity  of  the  calculations  tz.  assump- 
tions made  about  the  boundary  layer  flow  In  the  chamber. 


While  of  major  Importance  to  local  heat  transfer  calculations  these  dls-  j 

crepancles  in  die  throat  region  are  not  serious  from  the  standpoint  of  performance 
calculations.  For  most  rocket  motors  T8L,  or  related  Integral  methods,  should  ’ 

be  able  to  calculate  boundary  layer  losses  to  within  about  20%.  Considering 
that  boundary  layer  losses  are  usually  about  1%,  and  rarely  exceed  2%  (small 
motors,  of  short  firing  duration  being  one  of  the  exceptions),  the  total  impact  of  j 

such  a 20%  error  on  overall  performance  loss  should  be  restricted  to  several  tenths  .5 

of  a percent;  an  amount  that  should  prove  tolerable  In  most  cases. 

There  are  more  accurate  methods  avail  able  for  calculating  boundary  layer 

performance  loss,  such  as  the  Mass  Addition  & 'indary  Layer  Program^)  and  the 
(28) 

BUMP  Program'  , These  methods,  are  based  on  finite  difference  or  other  num- 
erical solutions  of  the  turbulent  boundary  layer  partial  differential  equations. 

While  the  advances  achieved  In  recent  years  allow  such  solutions  to  be  achieved 

In  reasonable  amounts  of  computer  time,  the  trade-off  between  lnceresed  accuracy 

and  Increased  operational  complexity  (as  well  as  run  time)  did  not  appear  to  favor 

the  uae  of  a finite  difference  type  approach  In  the  SPP  Program.  This  conclusion  • 

follows  from  the  usually  restricted  size  of  the  boundary  layer  loss,  which  logically 

precludes  a msjor  effort  to  calculate  It. 

The  previous  discussion  outlines  why  It  is  felt  that  Integral  methods,  and 
TBL  In  o-tlcular,  are  adequate  for  the  task  at  hand.  More  accurate  methods  such 
as  MABL  and  BUMP  are  not  warranted  for  general  solid  rocket  motor  performance 
prediction.  However,  in  certain  cases  the  Increased  generality  of  programs  like 
these  may  bo  the  only  way  -o  obtain  an  adequate  prediction  for  boundary  layer  per- 
formance loss.  If  such  cases  can  be  Identified,  a program  like  MABL  can  serve, 
not  as  an  Integral  pert  of  the  performance  prediction  methodology,  but  as  an  adden- 
dum to  It;  to  be  used  as  required.  On  the  other  hand,  there  will  be  classes  of  motors 
for  which  data  and  or  correlations  exist  which  will  adequately  be  able  to  characterize 
the  boundary  layer  loss.  This  will  be  especially  true  for  the  class  of  motors  for 
which  boundary  layer  losses  are  small.  In  such  cases  the  use  of  any  analytical 
boundary  layer  solution  may  be  obviated. 
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The  SPP  Program  has  been  structured  to  account  for  all  of  these  possibilities. 
The  TBL  program  was  selected  as  the  standard  method  for  computing  the  boundary 
layer  loss,  and  was  Incorporated  into  the  program  as  a module.  Two  other  methods 
for  detenu ining  the  boundary  layer  loss  have  been  explicitly  provided  for.  A simplified 
semi-empirical  boundary  layer  loss  correlation  can,  on  option,  be  utilized  in  lieu 
of  the  TBL  program.  This  correlation  is  described  in  Section  4.2.4,  The  boundary 
layer  loss  can  also  be  directly  input  to  the  SPP  Program.  In  such  a case,  the  input 
value  of  the  loss  is  used  in  deference  to  all  others.  In  this  manner,  should  the 
circumstances  warrant  it,  a more  exact  boundary  layer  method  may  bo  utilized  to 
calculate  the  boundary  loss  in  an  external,  uncoupled,  fashion.  The  computed 
value  may  then  be  Input  to  the  SPP  program. 

The  TBL  program  does  not  allow  for  a direct  assessment  of  the  time  de- 
pendent nature  of  the  boundary  layer  loss.  Por  motors  of  all  but  very  short  firing 
duration,  this  is  of  little  consequence  to  Hie  accuracy  of  the  overall  performance 
prediction.  Since  wall  temperature  is  an  input  to  the  TEL  program,  the  variation 
of  the  boundary  layer  loss  with  wall  temperature  (and  therefore  with  time)  may  bo 
evaluated  by  performing  a series  of  boundary  layer  solutions  wherein  the  nozzle  wall 
temperature  distribution  is  varied.  To  perform  this  calculation  in  a rational  man- 
ner a thermal  analyser  type,  transient,  heat  conduction  code  should  be  externally 
coupled  to  the  boundary  layer  solutions  to  provide  the  required  wall  temperature 
distributions. 

3.6.2  Modifications  to  the  TBL  Program 

Very  iittla  modification  was  required  to  Implement  the  TBL  program  as  a 
module  In  the  SPP  program.  Most  of  the  modifications  were  related  to  the  Incor- 
poration of  additional  labeled  commons  to  allow  TEL  to  communicate  with  ether 
modules,  as  appropriate.  Many  of  the  Items  that  must  normally  be  Input  to  TBL 
are  instead  obtained  through  this  linkage.  For  most  problems,  the  peed  to  specify 
the  bulk  of  the  remaining  Input  items  was  obviated  by  presett.ng  them  In  program 
data  statements.  Thus,  very  little  user  specified  Input  Is  normally  required  to 
execute  TBL  as  an  element  of  the  SPP  performance  prediction  methodology. 

The  only  ether  modification  of  substance  to  TBL  was  the  incorporation  of  a 
calculation  of  the  boundary  layer  induced  performance  loss  decrement,  *1 

spm 
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This  quantity  was  calculated  as  follows: 


sp. 


TBL 


A FBL 
mTD2P 
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where  jF^  Is  given  by  equation  (3-4S),  and  ls  the  no2zle  mass  flow  rate 

calculated  In  the  TD2P  module. 

The  TBL  calculation  Is  carried  out  at  the  average  chamber  pressure,  and 
the  initial  expansion  ratio.  The  boundary  layer  loss  is  assumed  to  be  Independent 
of  erosion  Induced  effects  for  the  present  purposes. 


3.7  Inter-Module  Interfaces 

In  order  for  the  various  programs  contained  within  the  SPP  program  to  function 
together,  automatically,  the  program  modules  must  be  able  to  communicate  with 
each  other,  as  required.  A perfunctory  description  of  the  Inter-module  Interfaces 
ls  given  here  to  provide  the  reader  with  an  overall  Idea  of  how  the  elements  of  the 
SPP  program  are  linked  together.  More  detailed  descriptions  of  the  program  linkage 
structure  may  be  found  In  Volumes  II  and  HI  of  this  report. 

The  following  five  simplified  diagrams  are  designed  to  Indicate  only  the 
nature  of  the  Information  which  flows  out  of  the  five  basic  computational  modules, 
and  Its  ultimate  destination.  In  actuality,  all  of  the  Inter-module  communications 
are  processed  through  the  Master  Control  Module.  Each  block  In  these  diagrams 
represents  a module.  The  following  abbreviations  are  used  for  convenience: 

Master  Control  Module:  MC 
Theoretical  Performance  Module:  ODE 
Grain  Design  and  Ballistics  Module:  GD&BC 
Xlnetlcs  Module:  ODK 
Two-Dimensional  Two-Phase  Module:  TD2P 
Boundary  Layer  Module:  TBL 
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4.  PERFORMANCE  LOSS  CORRELATIONS 


4.1  Combustion  Efficiency 
4.1-1  Approach 

The  combustion  efficiency  of  aluminized  solid  propellants  Is  ordinarily  close 
to  100  percent.  Departures  from  Ideal  characteristic  velocity  generally  result 
from  wall  losses  (as  can  exist  with  an  exposed  steel  chamber  In  subscale  test 
hardware),  from  particle  lags  In  the  throat  region,  and  from  dilution  by  ablatives 
(generally  corrected  out  In  data  reduction).  These  thr>  . , 'fects  will  be  further 
commented  upon  subsequently.  Exceptions  to  the  ef  '<*>  . combustion  of  aluminum 
are  found  with  low  burning  rate  propellants,  operation  at  low  chamber  pressure, 
and  tests  In  combustors  of  exceptionally  low  L*.  Interest  In  spacecraft  motors, 
upper-stage  motors  and  sustalner  operation  has  provided  substantial  data  showing 
the  loss  of  combustion  efficiency  at  low  chamber  pressures  and  burning  rates^4^. 

Therefore,  the  correlation  approach  began  with  a theory  that  alum',  am  combus- 
tion efficiency  Is  a function  of  propellant  burning  rate.  Burning  rate  bears  directly 
on  those  processes  which  Influence  the  competition  between  extensive  metal 
agglomeration  and  rapid  metal  Ignition  at  the  propellant  burning  surface^1  ^ . Burning 
rate  directly  influences  the  combustion  zone  temperature  gradient,  the  particle 
heating  rate,  i—1  the  gas  velocity  normal  to  the  propellant  surface.  Reported 
effects  of  propellant  lormuiatlon  ore  relatable  to  changes  In  burning  rate^4,  ^ 
cs  long  a'  ‘he  oxidizer  and  binder  provide  enough  energy  for  me.al  Ignition  (always 
true  in  ale  ..i.ilzed  propellants  ci  oractlca!  Interest,  not  always  true  for  other  metals). 
Likewise,  reported  effects  of  pressure  are  relatable  to  burning  ‘ate.  The  motor  L* 
does  not  appeal  to  be  Important  so  long  as  It  Is  greater  than  approximately  50  In., 
according  to  experiments  measuring  complete  combustion  dlrectly(16),the  Implication 
Is  that  complete  combustion  is  largely  determined  by  what  happens  at  or  near  the 
propellant  surface.  Motor  C*  efficiency  data  were  accumulated  Ir  several  groups 
for  the  purpose  of  constructing,  verifying  and  refining  the  correlation  with  burning 
rate,  and  finally  to  demonstrate  the  completed  correlation  on  an  aprlorl  basis. 

Statistical  analysis  Is  a usef  .l  diagnostic  tool  to  provide  Insight  Into  the 
nature  of  poorly-understood  phenomena,  but  has  proven  misleading  when  used  to 
generate  correlations  for  refined  predictions.  The  present  computr-r  program  devel- 
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opmeot  was  predicated  upon  the  premise  that  factors  affecting  the  specific  impulse 
efficiency  of  aluminized  propellants  are  not  poorly  understood,  and  that  existing 
theories  and  models  can  be  evaluate*  and  combined  tc  achieve  the  desired  accuracy. 
Further,  statistical  analysis  is  not  warranted  when  the  data  base  is  insufficient 
compared  to  the  potential  scope  of  application.  While  die  amount  of  BAXES  motor 
data  Is  considerable,  tee  fact  teat  it  is  a standard  test  motor  becomes  restrictive; 
potential  governing  parameters  will  vary  together  in  a corresponding  manner  and  so 
cannot  be  meaningfully  separated.  It  ts  also  true  teat  tee  great  majority  of  the 
BATES  data  are  at  a standard  test  condition,  and  cover  a fairly  narrow  range  of  bum 
ing  rates  and  aluminum  concentrations.  Therefore,  tee  use  of  statistical  analysis 
was  limited  to  a supplementary  role. 

Data  groups  that  were  studied  are  listed  in  Table  1. 

4.1.2  The  Basic  Correlation  with  Burning  Kate 

The  basic  correlation  of  C*  efficiency  with  burning  rate  is  shown  In  Figure  4-1. 
These  data  consist  of  Groups  I and  II,  excluding  data  for  propellants  having  aluminum 
concentrations  below  10  percent  or  greater  than  24  percent.  Each  point  represents 
an  average  for  a series  of  flrii,gs  of  a given  propellant,  in  a given  motor,  at  a par- 
ticular test  condition.  Where  BATES  data  for  a given  propellant  existed  in  g re  ps, 
each  group  was  averaged  separately.  It  is  observed  that  th  bulk  of  tho  data  exist 
between  burning  rates  of  0.3-0. 6 in./sec.  and,  within  teat  range,  the  correlation 
Is  not  meaningful  at  first  glance.  It  is  only  when  a broad  rar  ,*e  of  burning  rates 
are  examined  that  tee  correlation  appears  significant. 

It  is  observed  teat  C*  efficiency  falls  rapidly  below  burning  rates  of  0.25 
in./sec.  This  effect  persists  In  tee  full-scale  motors  which  are  induced  in  tnat 
range.  Therefore,  C*  efficiency  Is  not  expected  to  depend  upon  scale.  The  pro- 
blem with  the  steep  behavior  at  low  burning  rate  Is  teat  the  uncertainty  In  predict- 
ability becomes  large  In  that  regime.  At  high  burning  rate,  the  efficiency  appears 
to  approach  100  percent. 

4.1.3  BATES  Motor  Reproducibility 

The  basic  correlation  with  burning  rate  was  examined  in  tee  light  of  tee  repro- 
ducibility of  the  BATES  motor  C*  efficiency  data,  particularly  in  view  of  the  above- 
mentioned  confinement  of  most  of  the  data. 
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Table  4-1 


DATA  USED  TO  DEVELOP,  REFINE  AND  VERIFY 
C*  EFFICIENCY  CORRELATION 


GROUP  I DATA  - CONSTRUCTION  OF  CORRELATION 


a)  Propellants  In  6 x 11.4  Rohm  & Baas  Motors 

LPC-343-1401 

LPC-343-1402 

LPC-343-1403 

LPC-58QA 

LPC-S&2A 

LPC-586A 

I.T-C-62SC 

NOTS  X-32 

RHP-112 

TPH-1011 

LPC-36S-1SAL  (NAS  1-9182) 

NASA/Langley  (NAS1-10956) 


b)  Full-Scale  Motors 


ASTROBEE-F  (Susialner) 

Aerojet 

AVANH 

Lockheed 

Company  records 

ATS  (Etvi- Sumer) 

NASA/JPL 

Ref  ( 29 ) 

ATS  (Internal- Burner) 

nasa/jpl 

Ref  ( 29 ) 

HARPOON  (Lockheed  Testwelght) 

Lockheed 

Company  records 

1S6-S 

Lockheed 

Ref  ( 30 ) 

GROUP  H DATA  - BATES  MOTORS,  VERIFICATION  AND  REFINEMENT1 
a)  Conventional  Propellants 


Ref  ( 17 ) 

Ref  (17) 

Ref  ( 17 ) 
Company  Records 


Ref  (14) 
Ref  ( 18) 


RHP-112 

LPC-S43A 

RHP- 113 

LPC-543D 

RHP-161 

LPC-580A 

RHP- 162 

RDS-S01 

RHP-163 

RDS-507 

TFG-3016D 

RDS-510 

TPH-10RS 

UTF-11475 

TPH-1011 

CSSA 

TPH-1066 

VBB 

TPH-8163 

SPIS-31 

AHP-2969 

SPIS-32 

ANB-306F 

ANB-3105 

AAB-331B 

Table  4-1  (Continued) 


b)  Advanced  or  Experimental  Propellant* 


ARCHT-S015 

RHU-105 

ARCHT-9018 

DLH-258 

ARCHT-9021 

ANB-3406 

ARCHT-9024 

ANP-3413 

ARCHT-9027 

ARCEF-147 

ARCHT-9030 

LXA-100 

UTT-1S908 

SPIS-26 

UTX-14617 

SP1S-27 

UTX-14682 

UTX-14631 

UTP-1S151 

VLU 

VLZ-IJ 

VMO 

' m DATA  - DEMONSTRATION  OF  CORRELATION 

2.75  FFAR 

Aerojet 

Ref  (31) 

2.75  FFAR 

Hercules 

Ref  (32) 

SRAM 

Lockheed 

Ref  (33) 

HYDAC 

Lockheed 

Ref  (34) 

HI-EX  (Experimental) 

Lockheed 

Ref  ( 35) 

MINUTEMAN  III  3rd 

Aerojet 

Ref  (36) 

FW-4S 

UTC 

Ref  (37) 

SURVEYOR  RETRO 

Thlokol 

Ref  i 38) 

ANTARES  I 

Hercules 

Ref  ( 39 ) 

ANTARES  H 

Hercules 

Ref  ( 40 ) 

SPARROW 

Rocketdyne 

Ref  (41) 

SAM-D 

Thlokol 

Ref  ( 42 ) 

-°ARTAN-2nd 

Thlokol 

Ref  ( 43 ) 

156-5 

Lockheed 

Ref  (30) 

FIGURE  4-1  BASIC  CORRELATION  OF  C EFFICIENCY 


BURNING  RATE  (IN'SEC) 


Within  the  range  of  burning  rates  0. 3-0.6  In/sec,  the  correlation  exhibits 
a maximum  deviation  of  1.2  percent  and  an  average  deviation  of  0.6  percent.  All 
groups  of  SAKS  data  were  examined  for  reproducibility.  The  maximum  deviation 
within  a group  was  found  to  be  as  high  as  2.0  percent,  and  the  average  deviation 
within  a group  was  found  to  be  as  high  as  1 .4  percent.  These  results  were  not 
flash-ln-the-pan.  One  group  consisting  of  eight  nominally  Identical  firings  exhibited 
an  average  deviation  of  1.2  percent,  and  average  deviations  covering  all  groups 
were  typically  greater  than  0.5  percent.  The  accumulated  average  deviation  was 
0.8  percent.  Thus,  the  accuracy  of  the  correlation  could  be  said  to  be  comparable 
to  the  reproducibility  of  the  data.  This  was  not  particularly  satisfying  because 
unsound  correlations  might  turn  out  just  as  good  within  this  range  of  burning  rates, 
forcing  reliance  on  the  data  outside  of  this  range  to  support  the  proposed  correlation. 

Some  attention  was  devoted  to  the  reason  for  the  data  scatter  In  the  BATES 
motors.  The  principal  uncertainty  In  the  C*  measurement  Is  caused  by  nozzle  erosion; 
C*  requires  relating  the  measured  pressure  Integral  with  an  average  throat  area, 
a factor  not  present  In  the  specific  Impulse  measurement.  All  other  factors.  In- 
cluding Instrumentation  and  propellant  weight,  are  also  present  In  the  specific 
Impulse  measurement.  Average  deviations  In  specific  impulse  efficiency  ranged 
from  0. 1-0.3  percent,  clearly  and  significantly  superior  to  the  C*  error.  Reported 
nozzle  erosion  rates,  and  throat  averaging,  are  based  upon  pre-fire  and  post-fire 
measurements,  they  do  not  reveal  the  nature  of  the  throat  behavior  in  the  course 
of  the  firing.  A cause  of  variability  In  throat  erosion  behavior  Is  the  repetitive  use 
of  the  same  nozzle  Insert  for  economy.  It  was  confirmed  by  AFRPL  that  this  was  the 
general  practice  and,  within  a given  group  of  data,  one  can  observe  this  taking  place 
from  the  Initial  throat  sizes.  The  effect  Is  aggravated  In  high  erosion  situations. 
Therefore,  unless  this  data  scatter  In  C*  Is  acceptable.  It  would  be  recommended 
that  fresh  or  noneroding  nozzles  be  used  on  each  test  for  Improved  accuracy. 

The  Improved  accuracy  of  the  specific  Impulse  measurement  afforded  an 
alternative  opportunity  to  look  for  burning  rate  effects  through  specific  Impulse. 

4.1.4  Effect  of  Burning  Rate  on  Specific  Impulse  Efficiency 

Additional  confirmation  of  the  effect  of  burning  rate  on  performance  appears 
In  the  specific  Impulse  efficiency  data,  which  are  more  accurate  than  the  C*  efficiency 
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data,  and,  therefore,  lend  more  credence  to  the  effect.  This  Is  shown  for  Group 
□a  propellants  In  Figure  4-2.  Specific  Impulse  efficiency  Is  plotted  versus  aluminum 
concentration,  as  Is  often  done.  Each  datum  point  here  Is  an  average  of  an  over 
all  series  of  tests  for  each  propellant.  The  number  placed  next  to  the  datum  point 
Is  the  average  burning  rate.  Note  that,  at  each  aluminum  concentration.  Increasing 
burning  rate  Increases  specific  Impulse  efficiency.  This  Is  believed  to  reflect  C* 
efficiency  changes.  Therefore,  wide  variations  in  burning  rate  can  explain  why 
specific  Impulse  efficiency  cannot  cleanly  be  plotted  versus  aluminum  concentration. 

The  darkened  points  In  Figure  4-2  denote  one  type  of  propellant  In  which 
aluminum  concentration  was  systematically  changed,  with  burning  rate  approximately 
constant.  All  of  the  of  the  other  Ingredients  are  the  same.  Thus,  the  darkened 
points  would  basically  reflect  the  expected  change  In  nozzle  efficiency  with  Increa- 
sing aluminum  concentration. 

Assuming  the  burning  rate  correlation  to  be  valid  guide  to  C*  efficiency, 
the  nozzle  efficiency  was  then  computed  from  this  correlation  and  the  specific 
impulse  efficiency  data.  Results  are  plotted  versus  aluminum  concentration  In 
Figure  4-3.  Note  that  the  spread  of  data  points  Is  considerably  less  In  Figure 
4-3  than  In  Figure  4-2.  This  reflects  normalization  of  the  influence  of  burning 
rate. 

4.1.5  Effect  of  Binder  on  Specific  Impulse  Efficiency 

A curious  oddity  appears  In  Figure  4-3.  Except  for  the  propellant  In  which 
aluminum  - Mcentratlon  was  systematically  varied.  It  cannot  be  concluded  that 
efficiency  decreases  from  15  percent  aluminum  to  19  percent  aluminum.  Rather, 
the  data  Indicate  that  efficiency  Is  higher  at  16-18  percent  aluminum  than  at  15 
percent.  Upon  closer  examination  of  these  conditions.  It  was  found  that  the  anomaly 
could  be  explained  by  a binder  effect.  Data  from  15  percent  to  19  percent  aluminum 
could  be  systematically  organized  by  binder,  ar.d  the  more  efficient  binders  were 
concldentally  at  the  higher  aluminum  concentrations.  Within  the  Group  Ila  pro- 
pellants, the  order  In  the  direction  of  increasing  efficiency  was  polyurethane,  NC/ 
NG,  CTP3  and  HTPB,  PBAA  and  PBAN. 

This  binder  effect  was  Incorporated  as  a correction  on  the  correlation  of  C* 
efficiency  with  burning  rate.  By  correcting  C*  efficiency  for  binder,  the  Figure 
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FIGURE  4-2  CORRELATION  OF  SPECIFIC  IMPULSE  EFFICIENCY  WITH 


PER-CENT  ALUMINUM 


FIGURE  4-3  CORRELATION  OF  SPECIFIC  IMPULSE  EFFICIENCY  WITH 
ALUMINUM  CONCENTRATION,  NORMALIZED  FOR  BURNING  RATE 


:A3N3IOIJJ3  3TZZ0N 


PER-CENT  ALUMINUM 


4-3  data  should  collapse  Into  place  and  clearly  show  the  expected  dependence 
of  specific  Impulse  or  nozzle  efficiency  on  aluminum  content.  This  is  shown  in 
Figure  4-4,  which  incorporates  all  Group  n data  and  similar  corrections  for  Group 
nb  binders. 

The  value  of  the  C*  correlation  In  specific  impulse  efficiency  correlation 
is  summarized  as  follows: 

a)  If  specific  impulse  efficiency  in  BATES  motors  is  correlated  i .iy  with 
respect  to  aluminum  concentration,  the  standard  deviation  Is  0.8  per- 
cent and  the  probability  of  accuracy  within  2 percent  is  94  percent. 

b)  If  the  burning  rate  correlation  is  taken  Into  account,  the  standard  de- 
viation is  0.4  percent  and  the  probability  of  accuracy  within  2 percent 
is  better  than  99  percent. 

c)  If  both  binder  and  burning  rate  effects  are  taken  into  account,  the  stan- 
dard deviation  is  0.2  percent,  and  the  probability  of  accuracy  within 

2 percent  is  better  than  99  percent;  this  Is  comparable  to  the  reproduci- 
bility of  the  specific  Impulse  data. 

Thus,  the  BATES  motors  were  successfully  correlated  within  the  limits  of 
the  data  received. 

4.1.6  Effects  of  Extremes  of  Aluminum  Content  on  C*  Efficiency 

The  data  correlations  do  not  Include  propellants  having  aluminum  concentra- 
tions in  excess  of  24  percent.  Such  propellants  were  tested  in  BATES  motors  as  a 
matter  / academic  interest,  but  are  too  fuel-rich  to  be  of  practical  interest.  As 
might  be  expected,  values  of  C*  efficiency  fell  below  the  Figure  4-1  correlation, 
metal  agglomeration  was  undoubtedly  extensive.  On  the  other  hand,  values  of 
nozzle  efficiency  increased.  The  Increase  In  nozzle  efficiency  could  be  explained 
by  a reduced  total  particle  concentration  resulting  from  the  fuel-richeness.  In  any 
event,  it  was  concluded  that  should  propellants,  of  practical  Interest,  having  high 
aluminum  concentration  ever  be  developed,  they  would  probably  obey  the  correlation. 
Therefore,  the  correlation  is  not  modified  for  high  aluminum  content. 

C*  efficiency  data  for  low  burning  rate  propellants  having  aluminum  con- 
centrations below  10  percent  fell  above  the  Figure  4-1  correlation.  Aside  from  this 
fact,  a correction  would  be  warranted  by  the  premise  that  combustion  inefficiency 
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FIGURE  4-4  CORRELATION  OF  SPECIFIC  IMPULSE  EFFICIENCY  WITH 
ALUMINUM  CONCENTRATION,  NORMALIZED  FOR  BURNING  RATE  AND  BINDER 


stems  from  the  aluminum;  therefore,  some  means  should  be  afforded  to  approach 
the  unmetalllzed  propellant  condition.  A correction  factor  In  the  form  of  an  equa- 
tion was  added  to  accomplish  this. 
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4.1.7  Effect  of  Scale  on  C*  Efficiency 

The  C*  correlation  should  be  unaffected  by  scale  Inasmuch  as  Group  I 
data  Include  large  motors.  However,  there  Is  one  factor  related  to  scale  which 
does  need  to  be  addressed.  A large  motor  is  likely  to  be  fully  Insulated,  whereas 

a subscale  test  motor  (Including  BATES  motors)  Is  likely  to  have  exposed  walls. 

(44) 

Measurement  of  heat  loss  In  BATES  motors  Indicates  that  a 1 percent  loss  may 
be  assigned  to  the  chamber.  A fully  Insulated  motor  would  have  negligible  heat 
loss  because  the  thermal  energy  in  the  insulation  is  recovered  by  ablation  until 
the  last  Instant  of  the  firing.  Therefore,  In  constructing  the  final  correlation,  all 
subscale  data  were  shifted  upwar  J by  one  percent,  and  a one  percent  loss  factor  was 
Included  In  the  correlation  for  an  uninsulated  motor. 

4.1.8  Treatment  of  Discharge  Coefficient  Effect 

In  addition  to  combustion  and  heat  losses,  C*  efficiency  also  embodies 
a two  phase  two  dimensional  discharge  coefficient  effect.  The  effect  of  Cp  on 
performance  is  calculated  within  the  TD2P  module,  so  the  danger  that  a performance 
loss  would  be  counted  twice  In  the  specific  impulse  calculation  was  presented. 

To  assure  that  this  would  not  happen,  the  following  approach  was  taken. 

!he  value  of  C*  efficiency  determined  from  the  correlation  was  accepted 
as  valid,  being  based  upon  motor  data.  In  the  overall  performance  calculation  the 
CD  related  portion  of  the  C*  efficiency  Is  Isolated,  and  removed.  The  remainder 
Is  a combustion  efficiency  factor  (and  heat  loss).  This  procedure  is  detailed  in 
Section  3.1.  In  the  Grain  Design  and  Ballistics  Module  the  C*  efficiency  Is  used 
directly  In  the  calculation  of  chamber  pressure. 

4.1.9  Verification  of  Co  relation 

The  final  C*  efficiency  correlation  consists  of  a table  look-up,  with  cor- 
rection factors  In  equation  form,  and  Is  presenter  ir>  Table  4-2.  A comparison  with 
motor  C*  efficiency  data  from  Group  III  is  presented  in  Table  4-3 . 
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Table  4-2 


C*  EFFICIENCY  FUNCTION 


”C*  = [ K + TT  (10°  ' 10  ] x b x 0 

percent  aluminum;  a = a for  a £ 10,  a •-  10  for  a > 10 
1.00  for  fully  Insulated  motor,  0,990  for  uninsulated  motor 
binder  constant 


NF,  C = 1.003 
PGA/NC,  C 1.008 
PBAA,  C = 1.006 
PBAN,  C = 1.006 
HTPB,  C = 1.003 
CTPB,  C = 1.000 
NC,  C = 0.998 
PU,  C = 0.992 


burning  rate  constant 


k 

r 

k 

SI. 4 

0.50 

98.6 

93.1 

0.60 

98.9 

94.0 

0.70 

99.1 

94.6 

0.80 

99.2 

95.1 

0.90 

99.3 

95.6 

1.00 

99.4 

96.0 

1.20 

99.6 

96.4 

1.40 

99.7 

96.7 

1.60 

99.8 

97.0 

1.80 

99.9 

97.7 

98.2 

a 2.00 

100.0 

'1 

i 


>100.0,  = 100.C 


Table  4-3 
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Group  III  motors  were  selected  bacause  of  certain  features  bearing  upon 
specific  impulse  predictability.  In  the  C*  context,  thev  also  afford  a wide  range  of 
burning  rates,  aluminum  concentrations  and  hinder  types.  The  range  of  variability 
of  3.3  percent  lr.  C*  efficiency  provides  a reasonable  basts  for  testing  the  method; 
clearly,  the  test  has  little  meaning  If  the  range  Is  comparable  to  the  accuracy. 

It  Is  observed  that  the  average  deviation  of  the  predictions  Is  0.4  percent. 
Although  the  predictions  are  uniformly  distributed  between  high  and  low  values, 
the  overpredlctions  tend  to  be  larger  In  magnitude.  This  would  be  consistent  with 
the  uncertainty  of  accounting  for  insulation  dilution  m the  actual  data.  Insulation 
dilution  was  not  considered  in  these  predictions  {this  effect  Is  discussed  In  Sec- 
tion S.3).  Because  these  motors  were  not  used  In  developing  the  correlation, 
and  the  accuracy  Is  a small  fraction  of  the  range  of  variability.  It  Is  concluded 
that  the  Table  4-2  correlation  is  satisfactory. 

4.2  Specific  Impulse  Efficiency 
4.2.1  Approach 

The  developme  it  of  simple  expressions  for  each  source  of  specific  Impulse 
loss  affords  a means  for  the  rapid  and  convenient  assessment  of  motor  delivered 
performance.  Such  expressions  therefore  provide  an  Important  and  useful  supple- 
ment to  the  analyr-col  computer  program  subroutines  which  deal  with  these  losses 
In  a more  exact  and  detailed  mat  er. 

In  order  to  enhance  the  scope  of  applicability  of  the  simplified  method,  a 
seim-anai  . ral  rather  than  statistical  correlative  approach  was  adootec.  The  dis- 
advantages of  the  statistic,.’  approach  were  discussed  In  t'  c context  of  combus- 
tion efficiency  correlations.  In  correlating  specific  impulse  efficiency,  there  is 
the  added  consideration  that  tt  is  a function  or  n.  merous,  Interrelated  processes. 
Depending  upon  the  propellant  and  the  motor  design,  different  physical  mechanisms 
may  dominate  the  loss.  Therefore,  correlations  of  delivered  specific  Impulse  or 
total  nozzle  efficiency  In  terms  of  alam.num  concentration  and  scale  parameters 
will  be  deficient  to  the  extent  that  the  .ndlvldu3l  loss  components  do  not  vary 
systematica’'  with  these  variables.  Ac.  appirac’'  that  addresses  *he  Individual 
loss  components  is  preferred,  and  a sec i-analytlcal  approach  Is  bettei  able  to 
d~  so. 


1 i 
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The  resulting  expressions  were  derived  from  analytical  or  seml-emplrlcal 
egressions  appearing  in  the  literature,  or  by  fitting  the  results  of  the  analytical 
computer  programs  used  In  the  SPF  program.  A group  of  motors  was  selected  that 
would  serve  to  highlight  various  loss  components,  and  furnish  a basis  for  their 
verification.  The  final  expressions  were  compared  to  the  measured  total  loss  for 
all  of  these  motors . 

The  motors  selected  for  this  exercise,  and  the  reasons  for  their  selection, 
are  listed  In  Table  4-4.  These  motors  represent  extremes  of  the  motor  design  vari- 
ables that  influence  delivered  performance.  The  nozzle  efficiencies  of  the  produc- 
tion motors  are  plotted  versus  aluminum  concentration  In  Figure  4-5,  and  are  com- 
pared to  the  Figure  4-4  BATES  correlation.  It  .s  observed  that  there  is  no  correla- 
tion with  aluminum  content.  It  is  also  Interesting  to  note  that  most  of  the  produc- 
tion motors  deliver  lower  performance  than  the  BATES  motor  even  though  most  are 
larger.  Therefore,  factors  other  than  aluminum  content,  or  size  per  se,  must  be  im- 
portant. These  factors  do  not  vary  appreciably  In  BATES  motors  used  as  a standard 
propellant  reference,  but  must  be  considered  in  predicting  motor  performance  In 
general.  The  factors  which  distinguish  production  motors  are  contained  within  the 
expressions  for  the  loss  components. 

The  seml-emplrlcal  efficiencies  described  In  Sections  4,2.2  through  4.2.7 

are  In  terms  or  per  cent  nozzle  thrust  coefficient,  Cp,  losses  Instead  of  fractional 

I efficiencies, 

sp 

4.2.2  Divergence  Loss 

The  simple  half -angle  correction,  modified  for  contoured  nozzles^  * , is 
employ  Thus,  the  details  of  the  source  flow  from  the  throat  region,  and  the 
coupling  of  the  two-dimensional  particle  lags,  are  not  taken  Into  accojnt  here. 

This  Is  a potential  source  cf  error  in  the  simplified  analysis.  There  Is  insufficient 
Information  to  isolate  these  coupled  effects,  either  experimentally,  or  from  para- 
metric TD2P  solutions,  to  create  a Justifiable  simplified  expression.  One  dlt  .ulty 
In  performing  a parametric  analyses  Is  that  the  full  shape  of  the  nozzle  has  an  Impor- 
tant Influence  on  the  coupling  of  the  two-dimensional  and  two-phase  flow  effects. 

The  expression  used  for  divergence  loss  is  as  follows: 


"DIV-  * ™|»  -*«<— 7s4-)] 

a - half  angle  of  conic  i.  nozzle.  Included  angle  of  contoured  nozzle 
Spy1  exit  angle  of  co-.i  . i -ozzle 


a-tS 


FIGURE  4-5  ABSENCE  OF  PRODUCTION  MOTOR  PERFORMANCE  CORRELATION 


4.2.3  Kinetics  Loss 

The  kinetics  loss  expression  originally  proposed  simply  stated  that  this 
loss  is  1/3  of  the  difference  between  theoretical  shifting  and  theoretical  frozen 
performance  at  the  particular  expansion  ratio.  It  came  from  a study  of  various  pro- 
pulsion systems,  including  liquid  engines  and  fluorlnated  propellants  under  condi- 

(45) 

tions  that  would  emphasU^  this  loss  . However,  subsequent  comparison  with 
results  of  the  kinetics  calculations  from  the  ODK  program  indicated  that  a pressure 
correction  should  be  provided  to  reduce  the  predicted  loss  at  higher  pressures. 

The  final  form  of  the  kinetics  loss  expression  is  therefore: 


, = 33.3  1 - 


theoretical  frozen  Isp 
theoretical  shifting  Isp 


,200  , (4-2) 

• '*7~  T>200psl 


A thenr.ochemical  calculation  is  necessary  to  exercise  this  expression,  and  care 
should  be  taken  that  the  frozen  result  is  for  the  same  expansion  ratio  as  the  shifting 
result  (some  computer  programs  furnish  optimum  results  at  different  expansion 
ratios).  The  difference  between  shifting  and  frozen  performance  will  vary  with 
propellant  chemistry,  pressure  and  expansion  ratio. 


4.2.4  Boundary  Layer  Loss 

Boundary  layer  loss  is  expressed  in  the  form  of  a time-dependent  heat  transfer 
expression,  with  a correction  term  for  nozzle  expansion  ratio: 


’Ibi  =ci'^t|1  *2exp  (-c2p°-Vd°-2)J[ 


1 + 0.016  (( - 9)  j (4-3) 


P = pressure,  psi 
D(  = throat  diameter,  in 
t * time,  sec 
( = expansion  ratio 

Ordinary  Kozzle  : 


Steel  Nozzle: 


Cj  = 0.00365 
C2  = 0.0009J7 

C.  * 0.00506 

c2  * 0 


-.-vrv 
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The  time-dependence  Is  In  exponential  form  to  represent  transient  heat-up. 
Wall  heating  with  time  Is  an  Important  factor  to  consider  In  motors  of  short  bum 
time  {e.g.,  less  than  4 seconds).  The  heat  transfer  variables  are  combined  in  the 
well-known  form  for  pipe  flow.  The  dependence  on  expansion  ratio  represents 
the  effect  of  increased  nozzle  surface  area. 

The  bracket  term  containing  the  time-dependence  states  that  the  heat  trans- 
fer with  a cold  nozzle  (t=0)  will  be  a factor  of  3 greater  than  with  a fully  heated 
nozzle  (t -«*>).  This  factor  comes  from  comparing  a typical  equilibrium  wall  temp- 
erature-gas temperature  difference  with  the  difference  between  ambient  temperature 
and  the  gas  temperature.  It  is  consistent  with  the  differences  between  hot  wall 
and  cold  wall  heat  losses  derived  from  computer  calculations,  and  plotted  in  Ref. 

1 . The  time  constant  (C^)  comes  from  an  analysis  of  the  transient  heating  of 

a standard  BATES  motor.  The  coefficient  (Cj)  was  Judiciously  selected  based  upon 
a study  of  the  following:  a direct  measurement  of  heat  loss  in  a BATES  motor  (Ref. 

44  );  computer  calculations  of  boundary  layer  loss  (Ref.  1 );  measurement  and 
analysis  of  nozzle  heat  transfer  coefficients  (Refs.  46,47);  deductions  of  the 
boundary  layer  losses  In  standard  3KS  and  1GKS  motors  by  knowing  the  total  measured 
loss  with  unmetalllzea  propellant  (Ref.  48,49)  and  assuming  that  divergence  loss  and 
kinetics  loss  are  correctly  described  by  the  equations  above. 

The  constant  preceding  the  expansion  ratio  term  was  deduced  from  the  HI-EX 
motors  as  lollows.  The  motor  data  revealed  a decrease  m nozzle  efficiency  with 
increased  expansion  ratio.  A series  of  TD2P  calculations  provided  results  showing 
an  Increase  In  nozzle  efficiency  (decrease  In  2- D lag  losses)  at  high  and  Increasing 
expan.  . , n ratios*.  Assuming  the  expression  for  kinetics  loss  to  be  correct,  the 
remainder  >■  as  assigned  to  differences  In  boundary  layer  loss  with  expansion  ratio. 

A second  set  of  constants  is  provided  for  a "steel  nozzle,"  which  Is  primarily 
addressed  to  the  Hercules  2.75  motor.  These  constants  assume  a constant  cold 
wall  (Cj-O) , and  a higher  temperature  potential  (higher  Cj).  They  should  be 
limited  to  a motor  employing  a solid  steel  nozzle  of  relatively  thick  walls. 

*An  Important  practical  application  of  this  result  Is  that  propellants  of  higher 
metal  content  and  total  solids  ma,  .<  considered  for  high  expansion  ratio  than 
the  optimum  at  low  or  standard  e>-ans!on  ratios. 
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4.2,5  Ttoo-Phase  How  loss 


V \ 


The  expression  for  two-phase  flow  loss  Is  more  complicated,  attempting  to 
take  Into  account  different  regimes  of  concentration  effects  and  Interdependencies 
of  particle  size  and  throat  size  effects.  The  form  of  the  expression  and  the  values 
of  the  various  constants  were  obtained  by  fitting  detailed  solutions  appearing  in 
References  1,  50-52,  and  from  parametric  TD2P  solutions  obtained  during  the  pre- 
sent effort. 

The  two-phase  flow  loss  Is  given  by: 

. iVl 


”TP  a ^o.is^o.od  D Cg 

£ = mole  fraction  of  condensed  phase,  moles/100  gm 

Dp  = partie’e  size,  u 
P = pressure,  psl 
r = expansion  ratio 
D = throat  diameter.  In. 


If  r*0.09,  C4  = 

O.S 

Dt  '•  1:  C3  = 

9.0, 

C5 

= 

1.0, 

C6 

= 

1.0 

l<Dt<2:  C3  = 

9.0, 

C5 

1.0, 

C- 

o 

= 

0.8 

Dt-.2:  and  Dp  < 4: 

C3 

= 

13.4, 

cs 

= 

0.8, 

= 

0.8 

4 s Dp  <8: 

C3 

= 

10.2, 

C5 

= 

0,8, 

t 

0.4 

Dp.P: 

=3 

= 

7.58, 

S 

= 

O.S. 

C6 

0.33 

If  {<0.09,  = 

1.0 

Dt<l:  C3  = 

30.0, 

~5 

- 

1.0, 

C6 

= 

1.0 

l<Dt<2:  C3  = 

30.0, 

C5 

- 

1.0, 

C, 

o 

S 

0.8 

Dj  •,  2 and  Dp  < 4: 

C3 

44.5, 

C5 

= 

0.8, 

C5 

= 

0.8 

4<Dp-s8: 

S 

= 

34.0, 

S 

e 

0.8, 

C6 

0.4 

Dp -8: 

C3 

= 

25.2, 

C5 

= 

0.8, 

C6 

- 

0.33 

Dp  = 0.45-1  P1/3  r1/3  1 1 - exp 
L*  - motor  characteristic  length.  In, 


-0.004L*)  1+0 


,0451^  j 
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Hie  expression  branches  first  as  a function  of  exhaust  condensibles.  It 
then  branches  as  a function  of  throat  diameter,  and  then  with  particle  diameter,  for 
large  throat  diameters.  In  general,  the  dependencies  upon  particle  size  and  throat 
size  Increase  as  these  quantities  get  smaller.  There  Is  a slight  dependence  upon 
expansion  ratio.  Effects  of  pressure,  exhaust  condensibles,  throat  size  and  particle 
size  are  coupled.  In  that  particle  size  Is  a function  of  the  other  three;  the  result  Is 
that  the  loss  will  Increase  slightly  with  pressure,  strongly  with  exhaust  condensibles, 
and  go  through  some  orttmum  respecting  throat  size. 

The  expression  tor  particle  size  combines  theories  of  particle  growth  by 
condensation  in  the  chamber  and  collisions  in  the  nozzle,  as  discussed  In  Sec.  5.4. 

4.2.6  Submergence  boss 

The  expression  for  submergence  loss  Is  In  the  form  suggested  by  Kordtg  and 
Fuller  J , with  constants  modified  based  upon  the  data  appearing  In  Ref.  ( 1 ). 

P,  0.8  c0.4 

hSUB^  0-0584  ZfJ  (4'6) 

“t 

P « pressure,  pst 

£ = mole  fraction  of  condensed  phase,  moles/100  gra 

A*  = nozzle  entrance  area/nozzle  throat  area 
S = length  of  subT.crgcr.ee/iength  cf  internal  motor 
Dt  = throat  diameter,  in- 

4.2.7  bosses  Not  Corslderec 

Expressions  were  not  Included  for  the  following:  flow  separation  In  highly 
over -expanded  nozzles;  shock  losses  stemming  from  rough  nozzle  surfaces  caused 
by  nor, uniform  erosion  oi  ^palling,  expansion  ratio  losses  stemming  from  erosion  or 
nonopttmum  expansion  by  design,  "the  last  Item  Is  not  really  a loss,  and  may  be 
accommodated  by  comparing  predicted  and  measured  performance  at  the  actual 
overage  expansion  ratio.  It  has  s.gntficgace  lit  spriorl  prediction  where  the  actual 
throat  erosion  u not  Known.  Ar-el/sos  oi  shock  losses  and  nonuniform  erosion  are 
not  yet  veli-dsveiopce  tor  curler,  rrtdlctron.  Separa.iec  requires  considerable  over- 
expansion,  not  uxoiy  except  dur  , portions  of  Ignition  and  tall-off  which  are 
generally  :econpeqvonttai  awe  : ’ - total  Imculse. 


4.2.8  Verification  of  Simplified  Loss  Predictions 
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Results  of  applying  the  simplified  methodology  to  the  selected  gropu  of 
motors  are  listed  In  Table  4-5.  The  simplified  performance  predictions  were  ob- 
tained In  the  following  manner.  The  Individual  nozzle  loss  per  cent  efficiencies 
given  by  equations  (4-1)  to  (4-4),  and  (4  S)  were  combined  to  yield  an  overall 
nozzle  Cp  efficiency, 

^CF  100  - (^bl  + "DIV  + nf3N  r nSUB  + 

In  performing  this  calculation  average  values  for  pressure,  expansion  ratio  and 
L*  had  to  be  estimated.  Particle  size  was  obtained  from  equation  (4-5).  A C* 
efficiency  was  obtained  from  Table  4-2,  and  combined  with  the  efficiency 
to  yield  an  overall  delivered  I efficiency, 

T)j  = t?c*  tjcf  (4-8) 

sp 

The  measured  values  of  I efficiency  were  obtained  by  dividing  the  measured 
sp 

delivered  specific  Impulse  (delivered  to  ambient  pressure)  by  the  theoretical 
specific  impulse  (at  ambient  back  pressure)  calculated  at  the  estimated  average 
expansion  ratio. 

It  is  observed  that  all  predictions  are  within  l.S  percent  of  data,  and  that 
the  average  deviation  Is  0.5  percent  over  a range  of  variability  of  5.9  percent. 

The  overpredlctions  and  underpredictions  of  efficiency  appear  to  be  equally  distri- 
buted In  number,  but  the  overpredlctlans  tend  to  be  greater  in  magnitude. 

Th  SATES  motor  results  are  added  to  the  production  motor  results  In  order 
to  comment  upon  the  performance  of  NF  propellants  relative  to  standard  propellants. 
These  results  indicate  that  the  supeilor  performance  of  the  NF  propellants  In  3ATES 
motors  was  an  artifice  of  the  particular  propellants  and  the  BATES  motor.  The  NF 
propellants  were  of  substantially  higher  burning  rate.  Improving  C*  efficiency  Ir. 
addition  to  the  Inherent  Improvement  from  the  binder  itself.  The  higher  burning  rate 
necessitated  a larger  throat  size,  which  In  Itself  unproved  the  ncczle  efflclen  y 
1 percent  due  to  lower  two-phase  flow  losses.  Thus,  the  Indtcatea  2 percent 
Improvement  would  not  be  prea'lctec  to  be  maintained  m a full-scale  motor  where 
bum  rate  and  throat  size  would  be  fixed  by  design.  The  Inherent  Improvement 
In  combustion  efficiency  could  also  ue  compromised  by  an  increased  kinetics 
loss  with  fluorine  at  high  expansion  rv.:-. 
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Range  of  variability  = 5.9% 
Average  deviation  = 0.5% 


S.  OTHER  CONSIDERATIONS 
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5.1  Burning  Rate 

5.1.1  Approach 

Motor  burning  rates  cannot,  as  yet,  be  predicted  on  a strictly  a priori 
basis.  Therefore,  burning  rate  predictions  must  be  based  on  experimentally  mea- 
sured strand  burning  rates  for  the  propellant  of  interest.  The  problem  of  calculating 
motor  burning  rate  then  reduces  to  accounting  for  the  factors  which  cause  It  to 
differ  from  the  strand  buring  rate 

At  a given  set  of  conditions  (i.e.  pressure,  temperature)  strand  burning 
rate  Is  a constant.  Motor  burning  rapes  are  design  dependent,  however,  being  a 
function  of  the  local  convective  and  radiative  heating  at  the  burning  siiface^' 

The  differences  between  motor  burning  rates  and  strand  homing  rates  are  often  large 
enough  to  Impact  pressure  and  thrust  to  a greater  extent  than  the  desired  tolerance 
level  of  five  percent;  hence,  these  differences  must  be  accounted  for.  The  sophis- 
tication of  the  analysis  used  to  convert  strand  data  to  motor  bum  rate  Is  dictated 
by  a compromise  between  accuracy  and  the  information  readily  available  for  an 
a priori  calculation. 

One  cf  the  more  comprehensive  coupled  analyses  of  motor  burning  rate 
appears  in  Sef.rence  13.  That  analysts  was  constructed  for  the  special  pijposes 
of  predicting  Ignition  delay.  Ignition  spikes,  the  success  or  failure  of  command 
termination,  and  delivered  Impulse  during  transient  start-up  and  shut-down  lr, 
multiple  su,k,-restatt  operation.  For  those  purposes,  a comprehensive  combustion 
and  energy  exchange  analysis  was  required,  Including  an  In  depth  heat  transfer 
solution  for  the  solid  phase.  The  accurate  prediction  of  steady-state  ballistics 
followed  as  a by-product  of  the  accurate  prediction  of  transient  ballistics.  The 
disadvantages  of  this  analysis  are  Its  complexity,  size  and  run  time.  The  run 
time  and  size,  with  a limited  number  of  axial  Increments,  are  comparable  to 
goals  established  for  the  entire  Grain  Design  and  Ballistics  Module;  coupling  It  to 
the  much  finer  grid  required  by  the  gram  design  solution  was  considered  imprac- 
tical for  a priori  predictions.  Therefore,  the  use  of  this  analysis  should  be 
limited  to  the  detailed  study  cf  transient  ballistics  for  which  It  was  created. 

The  present  program  Is  chiefly  ccncor  *»d  v ith  steady  or  quasi-steady  ballistics. 
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The  approach  that  was  taken  was  to  construct  a closed-iorm  expression 
consisting  of  an  erosive  term,  a radiation  term,  and  a transient  term: 


r = rQ  (1  + fir)  (I  +r‘) 


r = strand  rate.  In. /sec. 
o 

fir  = steady-state  corrections 
r'  = transient  correction 


5.1.2  Erosive  Burning 

An  extensive  review  of  erosive  burning  was  performed  by  LPC  In  associa- 
tion with  the  development  of  a computer  program  to  predict  the  ballistics  of  nozzle- 
19  54) 

less  rocket  motors  ' .All  of  the  erosive  burning  laws  were  found  to  be  deficient 
In  some  respect,  but  the  best  and  most  comprehensive  equation  appears  to  be  that 
of  Lenoir  and  Roblllard: 

C1G°’8  r2°sr 

ar  - e*p  <-  -V— > (5-2> 

Cj,  c2  = constants 

G = local  port  mass  flux,  lb/ln?  -sec. 

x = local  axial  position.  In. 

o = propellant  density,  lb/ln“ 


In  this  v ork,  the  r appearing  In  the  exponential  was  changed  to  rQ  to  remove 
the  iterative  step:  It  Is  a minor  correction,  within  the  uncertainty  In  the  values  of 
Cj  anc.  The  dimensionless  c2  was  left  at  53  per  Lenolr-Roblllard.  The  constant 
Cj  may  be  expressed  In  analytical  form,  but  then  includes  other  corstants  of  more 
or  less  uncertainty.  The  value  of  this  constant  was  theiefore  established  from 
selected  motor  data  wherein  erosive  burning  would  be  expected  to  be  significant: 
0.0068  In  Inch"’®-  pound"’®  - second’®  units. 


5.1.3  Radiation  Contributions 

Radiation  can  be  a significant  source  of  combustion  augmentation  in  large 
(13) 

motors  . In  the  case  of  unmet -tlilzed  propellant,  radiation  will  be  imparted  by 
hot  nozzle  and  Insulation  surface;  as  well  as  by  the  hot  gases  flowing  inside 
the  motor.  In  the  case  of  metal,..  «.l  propellant,  the  particle-gas  cloud  shields 
the  propellant  from  the  hot  inc.  >ut  is  Itself  a source  of  significant  radia- 

tion heating . 
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Radiation  contributions  to  motor  burning  rate  have  been  studied  using  the 
sophisticated  transient  analysis  computer  programs  of  reference  (13)  which  dealt  with 
motors  containing  unalumlnl^ed  propellant  and  Reference  (55)  dealt  with  motors 
containing  aluminized  propellant.  Closed  form  expressions  are  constructed  as 
ratios  of  radiation  flux  to  combustion  flux:  a)  in  the  case  of  unaluminized  propellant. 


rF 

0.1 

1 ^ 

1.125  T*f  - 

)(Aj/Ab)  + 

(A^/Ag) 

+ 

l-exp(-xP) 

(Tp  - T* ) 

(5-3) 


b)  in  the  case  of  aluminized  propellant 

[l-exp(-\  P}]  (Tp  - Tj) 

l.Bps^p 


(5-4) 


where: 


= Stefan- Boltzmann  constant,  3.52  x 10 


-14 


BTU 


ln.a  -sec-6  K* 

average  view  factor  from  walls,  based  on  Ref.  (13)  computations 
thermochemical  flame  temperature,  °K 

propellant  surface  temperature,  °K=850  for  AP  composite  propel- 
lants, 523  for  double-base  propellants 

fourth  power  of  ratio  of  insulation  wall  temperature  to  propellant 
flame  temperature,  taken  from  typical  Ref  (13)  output 

exposed  insulation  area,  in? 

propellant  burn  surface  area,  in? 

projected  annular  nozzle  area,  radiating  surface,  in? 


0.1 

TF 
Ts 

0.125 
A- 

an 

0.410  = nozzle  analogy  to  0. 1 25  above 
= emlsslvlty,  (lb. /in?)  * 


X 

P 

Cs 

1.8 


BTU 


pressure,  lb/ln? 
propellant  heat  capacity, 

conversion  factor 


and 


X = 


/14.69 


(5-5) 


0.001  + 0.04  (Percent  Alumlnum)/100 
from  Reference  (55). 

Results  obtained  for  motors  / heroin  the  radiation  contribution  was  signi- 
ficant Indicates  that  equations  (5-  *'  ' -r0  yield  good  results,  without  having 

to  apply  any  additional  "correction  - 
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5,1.4  Transient  Term 


Effects  of  P on  Instantaneous  burning  rate  were  approximated  using  the  ex- 
pression derived  by  Paul  and  Von  Elbe: 

r‘  = 27=*P  (5-6) 

o 

n = pressure 

a = thermal  diffuslvity,  ln?/sec 
P = rate  of  change  of  pressure 

5.1.5  Other  Effects 

Effects  of  acceleration  on  burning  rate,  brought  about  by  spin-stabilization 
of  motors,  were  not  included  in  the  program.  If  desired,  such  effects  can  be  in- 
corporated at  a later  date. 

Burning  rate  augmentation  along  bondlines  has  been  encountered  in  certain 
motors.  The  phenomenon  has  been  attributed  to  propellant  nonuniformity.  No 
attempt  was  made  to  account  for  it  because  It  is  to  be  avoided  in  the  future. 

Recent  analyses  have  begun  to  address  two-dimensional  erosive  burning. 

In  these  analyses  bum  rate  differences  between  the  tips  and  valleys  of  stars  and 
keyhole  slots  are  accounted  for  in  addition  to  longitudinal  differences.  Although 
progress  a'  ing  these  lines  is  encouraging,  the  incorporation  of  two-dimensional 
erosive  burning  into  the  program  would  have  required  further  new  developments, 
and,  therefore,  was  beyond  the  scope  of  the  present  effort. 

Conditioning  temperatuie  is  accounted  for  through  program  input,  l.e., 
the  input  strand  burning  rates  should  correspond  to  the  conditioning  temperature 
of  interest. 

5.1.6  Strand  Data  as  Input 

In  recognition  of  the  fact  that  the  burning  rate  pressure  exponent,  n.  Is 
not  a constant,  and  that  conditioning  temperature  effects  vary  with  pressure  and 
temperature,  It  was  decided  to  Inn.;:  strand  burning  rate  data  in  tabular  form. 
Logarithmic  interpolation  Is  used  within  this  table. 
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S . 2 . 1 Approach 

Changes  In  nozzle  throat  area  with  ttir.e  mainly  Impact  the  chamber  pressure. 
Thrust  and  specific  Impulse  will  be  affected  to  a lesser  extent,  but  the  effect  can 
become  significant  If  the  throat  area  changes  appreciably. 

Because  nozzle  erosion  Is  a time-dependent  aerothermochemlcal  effect.  It 
Is  not  sufficient  to  correlate  erosion  by  steady-state  heat  transfer,  parameters 
which  would  Include  pressure,  throat  size  and  temperature.  The  time-dependence 
of  nozzle  heating  rates  must  be  accounted  for,  particularly  for  short  bum  time, 

A parameter  related  to  the  oxidation  potential  of  the  combustion  products  Is  also 
required.  Combustion  of  contemporary  propellant  formulations  produces  exhaust 
gases  which  contain  several  chemical  species  that  are  reactive  with  the  free  carbon 
contained  In  throat  materials.  It  follows  that  these  additional  considerations  will 
vary  with  the  formulation  of  the  propellant  and  the  type  of  throat  material. 

The  prediction  of  erosion  rates  In  carbonaceous  and  graphitic  throat  materials 
has  been  the  subject  of  several  computer  program  developments  combining  chemical 
kinetic  theory  and  heat  transfer.  These  efforts  are  very  sophisticated  and  Include 
experimental  work  as  well  as  theoretical  analysis.  However,  when  applied  to  motor 
firings,  accuracies  could  net  confidently  be  expected  to  be  better  than  30  percent 
without  a posteriori  adjustment  of  constants. 

According  to  the  original  program  plan  the  nozzle  erosion  rate  was  to  be 
specified  by  means  of  a user  determined  Input  table.  During  the  development  of 
the  computer  program  It  was  agreed  that  provision  for  some  means  of  calculating 
the  erosion  rate  would  be  desirable,  since  not  all  potential  users  could  be  expected 
to  have  the  capability  to  Independently  determine  It.  Rather  than  burden  the  SPP 
program  with  additional  complexity.  It  was  decided  to  Incorporate  an  erosion  rate 
capability  on  a seml-emplrlcal  basis.  Towards  this  end,  the  CMA^^  and  ASTHMA^7^ 
computer  programs  were  run  for  a matrix  of  conditions  consisting  of  two  propellants 
(zero  aluminum  and  16  percent  aluminum),  three  pressures  (500  psi,  1000  psl, 

2000  psl),  and  three  materials  (carbon  phenolic,  ATT  graphite  and  pyrolyt'c  graphite). 
For  each  material  the  calculated  erosion  rates  were  fit  as  a tabulated  function  of 
temperature  and  parameters  characterizing  heat  transfer  and  combustion  products. 
Temperature,  in  turn,  Is  derived  from  tve  propellant  flame  temperature  and  e tlme- 
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dependent  heat  transfer  expression.  Final  adjustment  of  constants  followed 
application  to  motors  representing  each  material.  The  exact  nature  of  the  erosion 
rate  as  a function  of  time  has  beon  retained.  Thus,  when  warranted,  the  present 
correlations  may  be  overridden  by  data,  or  the  results  of  another  analysis,  deemed 
more  appropriate  for  the  problem  at  hand. 


5,2.2  Form  of  the  Nozzle  Erosion  Correlation 

Nozzle  erosion  rate  was  correlated  In  the  following  form: 


^ - i;  ac  ch 

(5-7) 

where  is 

a time  dependent  film  coefficient  having  the  form. 

c,/c  • 0.8  „ 

CH  " TCI  hr)  "O' 

Dt  * M 

(5-8) 

c3 

= a constant.  Including  conversion  factors 

Dt 

= throat  diameter.  In. 

m 

= propellant  weight  flow  rate,  lb/sec 

\ 

= throat  area  In? 

X 

= gas  thermal  conductivity,  BTU/ln.-sec.-°R 

4 

= gas  viscosity,  lb/ln,-sec 

= gas  heat  capacity,  ft/"R 

* t 

= erosion  rate,  in. /sec 

°rp 

= material  density,  lb/ft3 

8C 

= dimensionless  erosion  rate  parameter 

CH 

= heat  transfer  parameter,  lb/fts-sec 

12 

= conversion  factor 

For  each  of  the  three  materials  considered,  equation  (5-7)  yields  the  relation. 

\ = 

0.173  Cjj  Carbon  Phenolic 

(5-9) 

K = 

0.109  S'c  CH  ATI  Qraphlte 

(5-10) 

^t  = 

0.086  8J.  CH  P>Tolytlc  Graphite 

(5-11) 
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(56,57) 

Hie  Aerotherm  computer  programs  calculate  g'c  as  a function  of  wall 

temperature,  for  a given  propellant  and  material.  The  results  were  found  to  be 
dependent  on  temperature  and  the  corrosivity  of  the  propellant  combustion  products, 
but  were  essentially  Independent  of  pressure.  Thus,  the  erosion  parameter  ^ , was 
correlated  In  the  following  form, 

fi'c  = c4(D  l (3)  (5-12) 

where  c^(T)  Is  a tabular  function  of  temperature,  and  f(8)  is  a function  of  the  pro- 
pellant products  corrosivity  parameter,  9.  The  temperature  in  equation  (5-12) 
should  be  wall  temperature,  however,  wall  temperature  was  converted  to  flame  temp- 
erature by  assuming  the  equilibrium  wall  temperature  to  be  70  percent  of  the  flame 
temperature.  The  flame  temperature  is  obtained  from  a thermochemlcal  equilibrium 
solution  (ODE  program).  The  time  dependence  of  the  wall  temperature,  expressed 
as  a flame  temperature.  Is  obtained  from  a standard  form  of  transient  heather  “re- 
pression: 

Tg  = To  + <Tf  * T0>  j 1 * <>XP  -CS(Y,0’8t  <5‘13) 

At  time  equal  to  zero  the  temperature  is  equal  to  the  initial  temperature,  ToP>K), 
while  at  longer  times,  the  temperature  approaches  the  propellant  flame  tempera- 
ture, TpPK) . The  constant  c^  is  found  by  linearly  Interpolating  In  a stored  table, 
with  T^  as  the  argument.  There  Is  a separate  table  for  each  of  the  three  materials 
considered. 


Th 

meter  9: 


-orroslvlty  of  the  combustion  products  is  characterized  by  the  para- 


M 5 


Mg  = molecular  weight  of  gaseous  combustion  products 

Mq  = molecular  weight  of  :..utal  oxide 

A - concentration  of  metal  oxide,  moles/100  gm 

^ th 

X,  = mole  fraction  concentration  of  the  1 oxidizing  species. 

1 There  arc  5 oxidizing  species  considered,  H.O,  CO., 
1/2  Oz,  C,  OH  Z 2 
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The  parameter  8 Is  computed  by  the  ODE  program  and  transmitted  to  the  Grain 
Design  and  Ballistics  Module,  where  It  Is  used.  Based  on  the  matrix  of  calcu- 
lations that  were  carried  out,  the  following  relations  were  obtained  for  f(fl)  (see 
eq.  5-12): 

f(g)  = 1 Carbon  Phenolic  (5-15) 

j 2 (no  dependence  on  0) 

f(8)  = (^-)  ' ATJ  Graphite  (5-16) 

f(3)  ~ (j^24  Pyrolytic  Graphite  (5-17) 


Thus,  to  summarize,  for  each  material  the  erosion  rate,  rt.  Is  found  from  the  one 
of  the  equations  (5-9)  - (5-11).  CH  Is  given  by  equation  (5-8)  and  0c  by  equation 
(5-12)  (with  the  subsldary  relations  (5-13)  - (5-17)).  The  constants  c^/Cp,  In 
equation  (5-8),  and  c$  in  equation  (5-13)  were  adjusted  based  on  comparisons  with 
a limited  amount  of  erosion  data  for  each  of  the  three  materials.  The  final  values 
selected  for  these  constants  are: 


c3/cp  = 1.682 
c$  = 3.368 


(5-18) 


Other  nozzle  materials  can  be  added  to  the  computer  program  by  repeating  the  pre- 
viously described  procedure  for  the  materials  of  interest. 

For  operating  pressures  below  300  psl,  the  phenomenon  of  nozzle  deposition 
rather  than  erosion  Is  likely  to  be  encountered  In  aluminized  propellants.  Depo- 
sition -z  not  been  studied  to  the  degree  of  erosion.  Limited  data  suggest  a 
rate  of  -0.001  In. /sec.  for  every  50  psl  below  300  psl.  This  is  offered  as  a guide 
in  constructing  an  Input  table  for  use  with  low-pressure  motors. 


Insulation  Degradation 


5.3.1  Approach 


Insulation  has  a minor  effect  on  the  Interior  ballistics  of  motors.  However, 
it  Is  enough  to  warrant  accosting  for  this  effect,  but  not  enough  to  conduct  de- 
tailed analyses.  Therefore,  a simple  methodology  was  adopted. 
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S.3.2  Performance  Degradation 


i ! 


From  the  standpoint  of  program  accuracy  and  efficiency  multiple  thermo- 
chemical  calculations,  as  the  Insulation  becomes  an  Increasing  fraction  of  the 
total  mass  flow,  would  not  be  Justified.  Instead,  one  matrix  of  thermochemical 
calculations  was  performed  varying  aluminum  content  from  0 to  25%  and  Insula- 
tlon/propellant  ratio  from  0 to  0.1.  The  following  expressions  were  fit  to  the 
results: 

-jjv  = ll  -[0.5  - 0.025  (%A1)]  R (5-19) 

0 = |l  for  %A1  i 20 

= 1 -[0.7  - 0.035(%A1)]  R (5-20) 

spo  = 1 for  %A1  > 20 
T 

= ll  - 1.9  R (5-21) 

Fo  = lo.Sl  for  Ri  0.1 

where  R denotes  the  lnsulatlon/propellant  weight  flow  ratio,  and  the  subscript  o 
refers  to  values  with  no  Insulation.  The  performance  degradation  due  to  Insulation 
Is  considered  only  In  the  Grain  Design  and  Palllstlcs  Module  and  is  not  transmitted 
to  the  Master  Control  Module. 

5.3.3  Determination  of  the  Insulation  Weight  Flow  Ratio 

Coupling  the  Insulation  exposure  to  the  grain  regression  pattern  would  not 
be  Justifies  ■>/  the  magnitude  of  the  Insulation  effect.  Rather,  a table  Is  provided 
for  a priori  Input  of  exposed  area  versus  time.  Ablation  rate  and  Insulation  density 
are  single  Inputs.  Therefore,  at  any  time,  the  flow  rate  of  Insulation  Is  determined 
and,  by  comparison  with  the  propellant  flow  rate,  the  value  of  the  ratio  is  established. 
The  thermochemical  parameters  used  m the  ballistics  analysis  are  then  degraded 
accordingly. 

5.4  Aluminum  Oxide  Particle  Size 
5.4.1  Approach 

The  two-phase  flow  loss  Is  an  Important  source  of  performance  loss  and 
Is  a strong  function  of  particle  size,  * the  ability  to  predict  particle  size  will 
have  an  Important  Impact  upon  the  s„:  e.  ,1  the  computer  program  over  a broad 
range  of  conditions. 


j 
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Considerable  particle  size  data  have  been  reported  In  the  literature.  Al- 
though the  collection  methods  and  the  relevance  of  the  results  have  been  open  to 
question  (58),  the  existing  data  must  be  accepted  in  lieu  of  anything  better.  A 
number  of  mathematical  models  of  oxide  particle  formation  and  growth  have  also 
been  published  In  the  literature. 

The  approach  that  was  taken  was  to  construct  a seml-emplrlcal  expression 
for  the  mean  particle  size.  The  form  of  the  expression  Is  analytical  in  origin, 
ano  the  constants  result  from  a fit  of  data  covering  a wide  range  of  conditions. 

A log-normal  distribution  Is  constructed  about  the  mean,  and  selected  sizes  repre- 
senting distribution  Intervals  are  used  in  the  TD2P  module. 

5.4.2  Experimental  Data  and  Pure  Empiricism 

(59) 

Aluminum  oxide  particles  were  first  collected  and  measured  by  Brown 
and  Sehgal^®®!  Brown's  work  provided  the  first  Indication  of  their  rounh-order-or- 
magnltude  size,  and  Sehgal  showed  the  size  to  be  a function  of  pressure.  Cheung 
and  Cohen^61*  extended  Sehgal's  work  to  show  dependendlces  no*,  only  upon  pressure 
but  also  upon  residence  time  and  aluminum  concentration;  no  effect  of  propellant 
aluminum  particle  size  was  detected.  Since  that  time,  the  Air  Force  has  Inves- 
tigated particle  size  In  a variety  of  motors which  has  formed  a basis  for  Air 

(44) 

Force  performance  Investigations  . Having  acquired  data  over  a wide  range  of 
motor  sizes,  the  Air  Force  chose  a simple  correlation  based  on  nozzle  throat  diameter. 

The  problem  with  a correlation  based  on  throat  diameter  is  that  It  does  not 
addrer  'he  controlling  physical  phenomena.  This  is  a drawback  of  empirical 
performance  correlations  based  on  size  parameters,  in  general.  The  reason  that 
such  correlations  work  in  an  approximate  sense  is  that  many  of  the  parameters 
tend  to  vary  together  (o.g. , throat  size,  motor  1* , propellant  weight  or  weight 
flow,  bum  time,  oil  tend  to  Increase  in  larger  motors).  However,  when  Improved 
accuracy  Is  desired.  It  is  necessary  to  accommodate  variables  outside  of  the  cor- 
relation range  Including  combinations  of  these  variables  which  are  exceptions  to 
general  trends  (e.g.,  a decrease  Ir.  L*  with  increasing  motor  size)  and  therefore 
sources  of  inconsistency  or  error.  Because  the  scope  of  computer  program  appli- 
cation may  encompass  a wide  range  of  variables,  singularly  and  in  combination, 
the  throat  size  correlation  would  not  be  expected  to  produce  the  desired  accuracy. 
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The  throat  diameter  correlation  was  constructed  from  data  over  a relatively 
narrow  range  of  pressure  and  aluminum  content,  but  over  a wide  range  of  sizes, 
wherein  L*  and  throat  diameter  tended  to  vary  together.  Effects  of  pressure  and 
particle  concentration,  as  well  as  L*,  must  be  Included  If  the  scope  of  propellants 
and  operating  conditions  Is  to  be  adequately  covered.  Cheung  and  Cohen  actually 
found  a decrease  In  particle  size  with  Increasing  throat  diameter  In  their  data 
because  they  increased  throat  diameter  In  order  to  reduce  pressure  (and  average 
L*)  In  fixed  hardware.  Similarly,  nozzle  efficiencies  reported  from  BATES  motor 
data^44'^'®^,  again  In  a fixed  size,  tend  to  infer  a smaller  particle  size  with 
decreasing  pressure  (increased  throat  size  and  smaller  mean  L*).  The  BATES 
data  also  indicate  a smaller  particle  size  with  decreasing  aluminum  content,  and 
with  changes  In  binder  type  that  would  result  In  decreasing  Al^O^  concentration. 

The  phenomena  of  oxide  particle  formation  and  growth  must  be  addressed  In  some 
manner  In  order  tc  account  for  these  detailed  effects. 

S.4.3  Models  and  Semi  Empiricism 

Mathematical  models  of  oxide  particle  growth  have  been  reviewed  by  Brown  ^S5^. 
The  Cheung  and  Cohen  paper  Included  a model  describing  particle  growth  by  con- 
densation and  agglomeration,  assuming  that  all  particle  growth  takes  place  within 
the  motor  cavity.  They  used  this  model  to  correlate  their  data  and  explain  the 
trends  observed.  They  then  extrapolated  to  large  boosters  and  predicted  the  sizes 
actually  measured  from  these  motors  two  years  later.  The  principal  effect  was 
residence  *t.re,  relatable  to  the  motor  L*.  The  model  also  predicts  particle  size 
to  be  dependent  upon  pressure  to  the  one-third  power  and  Al^Ch,  concentration  to 
the  one-third  power.  The  residence  time  dependence  Is  more  complicated  because 
condensation  and  agglomeration  have  different  functional  time-dependencies. 

Assuming  that  these  time  dependencies  can  be  combined  In  a single  term  represen- 
ted by  L*  to  a power,  there  results  the  following  seml-emplrical  expression: 


Dp~  P1'"3  $1/3  t*m  (5-22) 

This  expression  was  fit  to  the  Air  Force  data,  and  resulted  In  a very  good 
correlation  as  shown  In  Figure  5-1  (m=0.8).  The  correlation  was  superior  to  one 
simply  based  upon  throat  diameter.  However,  when  other  data  were  considered. 
Including  lag  loss  Implications  In  large  motors  having  small  nozzles,  the  L*  de- 
pendence was  concluded  to  be  too  . . 
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FIGURE  5-1  PARTICLE  SIZE  CORRE1  *TION  WITH  MOTOR  L * , CORRECTED  FOR 
VARIATIONS  IN  PRE;.  U^E  AND  ALUMINUM  CONTENT 
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Another  class  of  models  assumes  that  particle  growth  occurs  chiefly  In  the 
nozzle  by  collisions  resulting  from  velocity  gradients.  Such  a model  would  not 

have  an  Independence,  nor  does  It  address  pressure  and  concentration  effects 

(66) 

from  the  standpoint  of  physical  chemistry.  The  model  of  Crowe  and  Wllloughly'  ' 

Is  representative  of  this  class.  A seml-emplrlcal  correlation  of  the  data  was  con- 
structed based  on  this  model,  and  was  found  to  be  poor.  Therefore,  It  was  con- 
cluded that  a model  of  particle  growth  should  not  emphasize  growth  In  the  nozzle. 
This  Is  fortunate  from  the  standpoint  of  the  complexity  of  the  lag  loss  analysis. 

Finally,  a correlative  expression  was  constructed  combining  the  elements 
of  the  Cheung-Cohen  and  Crowe -Willoughby  analyses.  This  assumes  that  particle 
growth  occurs  in  both  the  chamber  and  the  nozzle,  but  each  to  a lesser  extent  than 
predicted  by  either  theory  alone.  It  also  serves  the  purpose  of  reducing  the  In- 
dependence of  growth  In  the  chamber  under  the  circumstances  where  It  was  found 
to  be  excessive.  The  expression  fit  to  the  data  Is  as  follows: 

Dp  = 0.454  P1/3  e1/3[l'- exp(-0.004L*)]  (1+0. 045  D{)  (5-23) 

Dp  In  microns 
P in  lb/ln? 

£ In  moies/100  gm 
L*  In  Inches 

D(  = throat  diameter.  In  Inches 

Thus,  a limiting  diameter  Is  reached  In  the  chamoer  at  large  1.*,  depending  upon 
pressure  c.‘  i concentration,  and  additional  grow*h  occurs  In  the  nozzle  depending 
upon  throat  size.  A comparison  with  the  data  Is  shown  tabulated  In  Table  5-1. 

The  average  deviation  is  18  percent,  which  represents  a variance  In  specific  Impulse 
efficiency  of  less  than  0.5  percent. 

Additional  studies  of  particle  formation  and  growth  are  in  progress  In  the 
context  of  particle  damping  of  combustion  Instability.  Observed  trends  support 
the  theory  that  substantial  growth  must  occur  In  the  motor  cavity.  Results  of 
these  studies  may  be  useful  In  future  modification  of  the  above  correlation. 
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Table  S-l 

COMPARISON  OF  PREDICTED  AND  MEASURED  PARTICLE  SIZES 


Range  of  variability  = 12. 8U 
Average  deviation  = 17.7% 
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6.  COMPARISONS  WITH  MOTOR  FIRING  DATA 


While  the  SPP  program  was,  by  design,  based  almost  entirely  on  pre- 
existing, proven,  analytical  and  semi-empirlcal  techniques,  there  was  no 
guarantee  that  the  Integration  of  all  of  these  methods  was  properly  implemented. 
Nor  could  It  be  known,  ahead  of  time,  tf  the  combined  techniques  would  be 
capable  of  meeting  the  accuracy  goals  which  had  been  set  (i.e.,  c-llvered 
Isp  to  + 2%,  delivered  thrust  and  total  impulse  to  + S%):  for  a wide  range  of 
motors , propellants  and  operating  conditions . 

In  order  to  validate  the  present  computer  program,  a series  of  compari- 
sons with  full  scale  motor  firing  data  was  carried  out.  To  date,  four  program-data 
comparisons  have  oeen completed.  The  first  three  comparisons  were  basically 
used  as  vehicles  to  check  out  the  operation  of  the  computer  program.  The  three 
check  out  cases  were  selected  to  provide  a reasonable  range  of  motor  sizes,  pro- 
pellants and  operating  conditions.  The  motors  considered  are  listed  below, 
together  with  some  of  the  characteristics  that  led  to  their  selection . 

Thlokol  Extended  Delta:  complicated  three  dimensional  grain, 
relatively  large  size,  high  expansion  ratio. 

AeiQlet  2.75:  small,  short  bum  time,  multiple  nozzles,  high 
throt  t erosion. 

Lockheed  bRAM:  end  burner,  unusual  nozzle  shape,  high 
pressure,  low  expansion  ratio. 

As  a result  of  their  diversity,  these  motors  provided  a reasonable 
thorough  basis  for  evaluating  the  performance  of  the  program.  The  fourth  com- 
parison with  motcr  data  was  canied  out  after  the  previously  mentioned  comparisons 
were  completed.  This  was  a so-called  sealed  envelope  comparison.  The  Air 
Force  selected  the  motor  to  be  considered  and  transmitted  to  us  only  the  infor- 
mation required  to  execute  the  program.  The  motor  chosen  for  this  comparison 
was  the  C-4  Third  Stage.  It  represents  a motor  of  much  current  Interest,  having 
a high  energy  propellant,  complex  grain  geometry  and  reasonably  high  expansion 
ratio.  The  completely  a priori  performance  predictions  for  this  motor  were  then 
delivered  to  AFRPL  and  evaluated. 

The  results  of  each  of  these  four  comparisons  are  Individually  discussed 
be'ow.  The  results  for  all  of  these  cases  are  then  summarized  in  Section  6.5. 
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6.1  Extended  Delta  Motor 


I 


i 

i 


(67) 


The  present  results  for  the  Extended  Delta  motor  were  obtained  by 
executing  the  complete  performance  prediction  methodology  on  a completely 
coupled  basis.  For  reasons  discussed  subsequently,  this  was  the  only  one 
of  the  first  three  check-out  cases  for  which  this  was  done. 

The  theoretically  calculated  delivered  specific  Impulse  for  this  motor 
was  found  to  be  0.46%  less  than  the  measured  value.  The  simplified  semi- empi- 
rical performance  prediction  correlations  under  predicted  the  measured  value  by 
0.79%.  The  results  of  the  ballistics  calculations  are  discussed  below. 

The  predicted  pressure-time  curve  for  the  Extended  Delta  motor  is 
compared  to  the  actual  pressure  trace  in  Figure  6-1.  It  is  observed  that  the 
qualitative  and  quantitative  features  of  the  data  are  well-predicted.  Pressure 
is  slightly  underpredicted  for  the  initial  portion,  and  slightly  overpredicted 
for  the  latter  portion.  The  action  time  is  somewhat  overoredicted. 

The  predicted  action  time  average  pressure  was  534  psia , compared 
to  the  actual  result  of  545  psia,  for  an  error  of  2.0%.  Average  thrust  was 
underprr  dieted  by  3.1%.  The  action  time  was  overpredicted  by  4.3%,  in  part 
because  of  a propellant  weight  error  attributed  to  the  fact  that  the  web 
dimensions  had  to  be  measured  from  a motor  drawing.  Correcting  for  the 
pressure  and  assigned  web  errors,  burning  rate  was  underpredicted  by  2.4%. 

This  burning  rate  error  .n  part  explains  the  low  initial  pressure  ana  longer 
action  time.  The  pressure  becomes  high  near  the  end  of  the  burn  in  part 
because  the  throat  erosion  was  underpredicted  by  18.7%  and  in  part  because 
of  the  icreased  input  web.  As  a result  of  compensating  effects  (i.e.,  thrust 
underpredlution,  weight  flow  overprediction)  the  predicted  total  impulse  was 
only  0.5%  above  the  measured  value. 

A burning  rate  error  of  2.4%,  if  uncompensated  by  errors  in  the  other 
parameters  governing  motor  pressure,  would  produce  a 3.4%  pressure  error  for 
a 0.3  exponent  and  a 6.1%  pressure  error  for  a 0.6  exponent.  Therefore,  an 
error  even  in  the  course  of  strand  burning  rate  measurement  could  make  it 
difficult  to  predict  pressure  to  within  5%.  Sources  of  compensating  (or  aggra- 
vating) error  would  be  in  the  grain  geometry,  nozzle  erosion  and  C*  efficiency 
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Predicted  end  ‘ctual  Pressure  Traces  For  Extonded  Delta  Motor 
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computations.  That  there  Is  a grain  geometry  erroi  Is  evident  lr.  Figure  6-1  In  that 
the  predicted  pressure  dip  precedes  the  actual  pressure  dip  even  though  the 
predicted  burning  rate  was  low.  It  follows  that  the  predicted  pressure  dip 
should  have  occurred  later  than  actual.  However,  the  prediction  of  the  basic 
ciiaracter  of  the  pressure  curve  shows  that  the  grain  design  program  Is  properly 
computing  the  geometry  In  general.  In  this  case,  the  error  In  the  throat  erosion 
partially  compensated  for  the  error  In  the  burning  rate. 

(13l 

6.2  SRAM  First  Pulse 


The  primary  reason  for  selecting  SRAM  as  a validation  case  was  to 
provide  a means  to  check  out  the  grain  geometry  and  ballistics  computations 
for  an  end-bumer.  The  motor  also  has  a unique  nozzle  geometry  and  r Igh 
operating  pressure.  A full  SRAM  firing  consists  of  two  pulses.  Calculations 
were  carried  out  only  for  the  first  pulse. 

The  current  grain  design  analysis  does  not  contain  a model  for  *he  meniscus 
burning  phenomenon  encountered  In  the  SRAM  motor.  Thu".,  the  ballistics  pre- 
dictions for  the  SRAM  motor  had  to  be  corrected  for  the  > A meniscus  burning. 
Neglecting  meniscus  burning,  the  computed  pressure  and  thrust  levels  were 
expected  to  be.  and  were,  low.  The  computer  output  was  corrected  for  meniscus 
burning  by  raising  the  pressure  at  each  instant  of  time  and  by  reducing  the  value 
ol  each  time  in  accord  with  the  Increase  in  effective  burning  rate.  Thus,  the 
fact  that  a given  KN  is  encountered  sooner,  was  accounted  for.  The  basis  for 
the  pressure  correction  at  each  instant  was  the  previously  known  effect  of  the 
phenomenon  on  the  grain  regression.  Although  the  important  result  was  the 
successful  operation  of  the  computer  program  for  an  end-bumer,  the  corrected 
pressure-ttrne  trace  is  shown  compared  with  data  in  Figure  6-2. 

The  shape  of  the  pressure  trace  is  well-predicted.  The  predicted  action 
time  average  pressure  was  low  by  1.5%.  The  predicted  motor  burning  rate  as 
corrected  was  low  by  1.3%,  and  was  compensated  by  an  underprediction  of  throat 
erosion  of  2.5%.  Action  time  average  thrust  was  underpredicted  by  0.8%.  The 
underprediction  of  pressure  was  compensated  by  a 0.9%  overprediction  of 
delivered  specific  Impulse,  which  is  well-within  the  program  goal. 

The  specific  Impulse  prediction  for  this  case  was  not  based  on  the  motor 
operating  conditions  output  by  theGraln  Design  and  Ballistics  Module,  due  to  the 
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neglect  of  the  previously  mentioned  me:  "is  burning  phenomenon.  Average  motor 

pressure  and  erosion  rate,  based  on  the  data  were  Input  to  the  program,  together 
with  an  average  motor  L*,  and  the  ODE,  TD2P  and  TBL  Modules  were  executed. 

The  Kinetics  Module  was  not  utilized  since  the  loss  is  Insignificant  at  the  SRAM 
operating  pressure  (the  simplified  kinetics  loss  correlation  predicts  a 0.1%  kinetics 
loss  for  this  case) . 

6.3  Aerojet  2.75  FFAR(31) 


The  ballistics  prediction  for  the  Aerojet  2.75  inch  motor  (Ref  31) 
encountered  an  anomaly.  As  shown  in  Figure  6-3,,  there  Is  no  resemblance 
between  the  predicted  and  actual  pressure-time  trace.  A disturbing  circumstance 
was  an  oscillograph  trace  for  a ”2.75  Inch  motor"  which  had  been  previously 
furnished  to  Ultrasystems  by  AFRPL,  and  which  showed  progressive  burning. 

Tils  piece  of  data  will  be  discussed  further  subsequently.  In  view  of  this 
startling  result,  an  Investigation  was  conducted  as  to  the  nature  of  the  problem. 

After  checking  the  inputs  to  the  computer  program,  a discussion  was  held 
with  members  of  the  LPC  Engineering  Department  who  had  2 . 75  inch  motor  exper- 
ience. It  was  reported  that  there  is  no  reasonable  way  in  which  a circular-port 
grain  could  show  a neutral  pressure  trace,  and  that  inhibiting  the  aft  portions 
would  have  little  effect  on  the  progressivity. 

Nt  rtheless.  Ref  (31)  was  checked  against  the  engineering  drawing  of 
the  grain  and  there  appeared  a basis  for  suggesting  a modification  to  the  Input 
description  of  the  Inhibited  surfaces.  A test  case  was  run  with  this  modified 
input,  but  showed  little  change  in  the  progresslvlty.  However,  because  a 
different  ansi  ' did  result,  it  was  possible  to  conclude  that  the  grain  geometry 
subroutine  was  at  least  iecognizlng  inhibited  surfaces  and  probably  correctly 
accounting  for  them. 

The  next  step  was  a conversation  with  Aerojet.  It  was  established  that 
the  neutral  character  of  the  pressure  trace  was  reproducible  over  an  extensive 
test  history.  The  explanation  offered  was  a combination  of  super  erosive  burning 
and  considerable  nozzle  throat  erosion.  However,  this  explanation  proved 
unsatisfactory  for  the  following  reasons: 
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Figure  6-3  Predicted  and  Actual  Pressure  Traces  For  Aerojet  2.75-Inch  Motor 
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a)  The  computer  program  adequately  predicted  the  reported  measured 
throat  erosion  (11.3%  error) . It  was  verified  that  a single  nozzle 
had  not  mistakenly  been  utilized  In  the  calculation  (the  motor  has 
four  nozzles) . 

b)  An  approximate  calculation  established  that  the  rate  of  throat 
erosion  would  have  to  be  Increased  by  a factor  of  2 In  order 
to  produce  a neutral  pressure  trace.  Discussion  with  LPC 
engineers  Indicated  that  a different  throat  material,  such  as 
an  asbestos-based  material,  would  be  required  to  produce 
such  a rate. 

c)  The  computer  program  closely  predicted  the  reported  measured 
motor  burning  rate,  when  corrected  for  the  difference  between 
predicted  and  measured  average  pressure  (1.4%  error). 

Therefore,  the  factor  increase  between  motor  and  strand  rates 
was  satisfactorily  predicted. 

d)  The  Initial  port-to-throat  diameter  ratio  of  the  motor  is  2 , which 
is  too  high  a value  for  exceptional  erosive  burning.  Moreover, 
the  aft  end  port,  which  would  be  expected  to  exhibit  the  highest 
eroslvlty , was  Inhibited. 

e)  An  approximate  calculation  established  that,  to  produce  the 
measured  average  burning  rate  and  a neutral  pressure  trace  with 
the  measured  throat  erosion,  the  initial  erosive  burning  would 
have  to  exceed  that  produced  In  nozzleless  rocket  motors  and 
the  burning  rate  G-dependence  would  have  to  be  equal  to  the 
pressure-dependence  (like  a hybrid  rocket  fuel  containing  some 
AP).  Both  of  these  requirements  are  unrealistic. 

f)  The  long  tall-off  exhibited  by  the  motor  Is  indicative  of  some 

anomaly.  In  seeking  to  explain  the  anomaly,  discussions  were 
held  with  LPC  engineers  regarding  the  possibility  of  unbonds. 
Inhibitor  failure  or  cracks.  An  unbond  was  considered  unlikely 
in  view  of  the  encapsulation  of  the  propellant  in  this  particular 
design.  Inhibitor  failure  would  not  significantly  alter  the  pro- 
gresslvlty,  as  mentioned  earlier.  The  possibility  of  a reproducible 
cracking  of  the  grain  was  then  assessed. 
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An  LPC  stress  analyst  was  aware  of  a reproducible  crack  occurring  at  the 
interface  between  an  Inhibited  surface  an  an  uninhibited  surface.  This  occurred 
In  an  Aerojet  Mlnuteman  motor.  The  motor  had  been  stress-analyzed  in  Its  Initial 
configuration,  but  not  when  a fillet  was  formed  from  burning  at  the  Interface: 

fillet 
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''-uninhibited  V 
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The  problem  manifested  itself  as  a motor  failure  some  time  into  the  bum,  and  there- 
fore, warranted  detailed  scudy.  (Such  Is  not  the  case  with  the  2.7S  inch  motor  In 
question,  which  to  all  Intents  and  purposes,  performs  as  required.)  Under  separate 
contract,  UTC  performed  a stress  analysis  of  the  Mlnuteman  Motor  taking  the  fillet 
Into  account  and  the  problem  was  subsequently  Identified  and  resolved. 


A series  of  curves  was  then  constructed  for  an  Idealized  cylindrical 
2.7S  Inch  motor  and  are  shown  In  Figure  6-4.  Die  first  curve  shown  for  the  absence 
of  port  Inhibition,  or  complete  Inhibitor  failure  on  the  poet.  The  second  is  for  a 
partial  Inhibition,  or  a part  failure.  The  third  Is  for  port  inhibition  per  the  Aerojet 
design  without  failure.  Die  fourth  curve  takes  ihe  Aerojet  design  without  Inhibitor 
failure,  and  adds  a diagonally  circumferential  crack  beginning  at  the  fillet  end 
propagating  at  a 45°  angle  to  the  wall.  In  each  case,  the  Is  adjusted  for  the 
measured  throat  erosion  by  assuming  a constant  rate  with  web  burned. 


Each  Kn  curve  Is  normalized  to  the  Initial  Kjj  of  the  Aerojet  design 
without  a crock  or  Inhibitor  failure.  It  Is  verified  that  the  port  inhibition  does 
not  significantly  alter  the  progresslvlty:  the  design  inhibition  is  about  10%  more 
progressive  than  no  port  inhibition.  However,  the  assumed  Ideal  crack  does 
produce  a near-neutrallty.  Moreover,  the  Initial  Kjg  was  determined  to  be  consis- 
tent with  the  Initial  pressure  of  the  motor  assuming  the  predicted  erosive  burning 
to  be  valid  and  not  a super  erosive  burning. 

The  ballistics  prediction  for  the  Aerojet  motor  is  qualitatively  compared 
to  the  pressure  trace  furnished  by  AFRPL  In  Figure  6-5 , The  AFRPL  motor  has  the 
3ame  grain  design  as  the  Aerojet  motor,  but  differs  from  the  Aerojet  motor  as  to 
propellant,  nozzle  configuration,  and  complete  lack  of  grain  Inhibition.  The  follow- 
ing significant  aspects  may  be  noted:  The  AFRPL  motor  show3  progresslvlty:  It 

•Ratio  of  propellant  bum  area  to  nozzle  throat  area. 


Figure  6-4  Diagnostic  KN  Curves  For  2. 75  Motor 
Including  Measured  Throat  Erosion 


Fractional  Web  Burned 


Figure  6-5  Comparloan  Of  Predicted  Pressure  Trace  For  Aerojet 
2.75  Motor  With  Actual  Trace  For  AFRPL  2.75  Motor 
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was  verified  that  its  progressivity  is  consistent  with  its  change  In  Km ; the  sharpness 
of  its  tail-off  is  closer  to  the  predicted  Aerojet  motor  than  the  actual  Aerojet  motor. 
The  AFRPL  motor  shows  less  progressivity  than  the  predicted  Aerojet  motor  because 
of  its  complete  lack  of  grain  inhibition.  It  is  apparent  that  the  AFRPL  motor  does 
not  exhibit  any  super  erosive  burning  with  high  G-dependence . and  does  not 
exhibit  any  anomaly. 

For  all  of  the  foregoing  reasons,  it  is  concluded  that  the  Aerojet  2.75 
inch  motor  exhibits  some  reproducible  anomaly,  that  this  anomaly  is  probably  a 
crack  emanating  from  the  inhibitor  Interface,  and  therefore,  it  is  impossible  to 
predict  the  observed  ballistics.  On  the  other  hand,  the  ballistics  of  the  AFRPL  2.75 
inch  motor  would  probably  have  been  successfully  predicted  were  that  the  test  case. 

Since  a satisfactory  ballistics  solution  for  the  2.75“  motor  could  not  be 
obtained,  quantitative  comparisons  with  the  measured  pressure  and  thrust  versus 
time,  and  total  Impulse,  ate  not  indicative  of  the  accuracy  of  the  program,  and  are 
not  given.  A valid  comparison  with  the  delivered  specific  impulse  was  obtained, 
but  in  a somewhat  non-aprlori  manner.  Using  one  of  the  many  operational  modes 
of  the  SPP  program,  the  measured  average  chamber  pressure  and  erosion  rates  were 
input,  together  with  an  estimate  for  the  average  motor  L*.  The  Theoretical  Per- 
formance (ODE),  Two-Dimensional  Two-Phase  (TD2P),  and  Boundary  Layer  (TBL), 
Modules  were  then  executed.  At  the  high  pressure  for  this  motor  the  execution 
of  the  Kinetics  Module  (ODK)  is  not  warranted  as  the  loss  is  quite  small,  but  the 
amount  of  computer  time  that  would  be  required  is  relatively  large.  The  simplified 
kinetics  loss  correlation  should  be  adequate  at  pressures  of  1000  psl,  or  more 
for  conventional  aluminized  propellants. 

Preceding  In  the  aforementioned  fashion,  the  predicted  delivered  I was 
found  to  be  Just  0.3%  above  the  measured  value  for  this  motor.  The  simplified 
performance  prediction  correlations  yielded  a value  0.93%  less  than  the  measured 
value  for  this  case.  Both  of  these  results  are  quite  satisfactory. 

6.4 Sealed  Envelope  Prediction;  C-4  Third  Stage 

This  prediction  was  a special  challenge  because  it  was  a blind  prediction 
of  a motor  containing  a complex  grain  geometry  and  a high-exponent  propellant. 

As  would  be  expected  for  such  a case,  the  complete  capability  of  the  program  was 
utilized  In  the  calculation.  Th"  specific  Impulse  prediction  for  this  case  was  also 
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excellent.  The  theoretically  determined  delivered  specific  Impulse  was  only  0.4% 
higher  than  the  reported  motor  firing  value,  while  the  simplified  performance  pre- 
diction correlations  yielded  a value  of  0.6%  on  the  low  side.  , ; 

The  predicted  pressure  trace  Is  compared  to  the  actual  motor  data  In  Figure  • 

6-6.  The  result  Is  similar  to  the  Extended  Delta  result  In  that  pressure  is  under-  . j 

predicted  during  the  Initial  portion  of  the  firing,  and  overpredicted  at  later  times. 

Because  the  pressure  exponent  Is  high,  the  extent  of  the  pressure  departures  Is 
greater  than  In  the  Delta.  The  grain  geometry  calculations  appear  to  be  better 
here  In  that  the  pressure  peak  and  dip  are  predicted  to  occur  at  a later  time,  which 
Is  consistent  with  the  predicted  slower  Initial  burning.  The  compensating  errors 
produce  a near-correct  prediction  of  action  time,  and  the  tall-off  Is  well-predic- 
ted. 

Action  time  average  pressure  was  underpredtcted  by  0.7%,  the  best  result 
of  all  test  cases.  Burning  rate  was  underpredicted  by  1.3%.  The  throat  erosion 
was  underpredicted  by  41%  of  the  average  throat  erosion;  fortunately,  the  mag- 
nitude of  the  throat  erosion  was  small  and  hence,  was  negligible  In  performance 
prediction.  It  Is  thought  that  this  error  In  erosion  rate  stems  from  extrapolation 
of  the  seml-emplrlcal  throat  erosion  constants  to  a new  high  temperature  regime, 
and  should  be  subject  to  adjustment  In  future  work.  Action  time  average  thrust 
was  underpredtcted  by  1.3%,  Total  Impulse  was  underpredicted  by  2.  !%•  this 
resulted  from  an  error  in  computed  propellant  weight,  on  the  low  side.  The  ln.’uts 
for  the  grain  geometry  had  to  be  prepared  from  a small,  relatively  crude  grain 
drawing,  and  strand  burning  rate  was  provided  to  us  at  only  one  condition.  In 
view  of  th.s,  the  accuracy  of  these  ballistics  predictions  Is  quite  good. 

A summary  of  the  results  for  this  case  may  be  found  in  Table  6-1.  A break- 
down of  the  individual  losses,  and  a comparison  between  the  theoretical  and  sim- 
plified predictions  Is  presented  In  Table  6-2. 

6.S  Summary  of  Results 

The  results  of  the  four  previously  described  comparisons  with  motor  firing 
data  are  summarized  In  Table  6-1.  E t for  the  anomalous  pressure  behavior 
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TABLE  6-1  SUMMARY  OF  RESULTS 


Motor  used  In  developing  or  evaluating  the  correlation 


Figure  6-6  Predicted  And  Actual  Pressure  Traces  For  Blind  Prediction  Test  Case 
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observed  in  the  2.75"  motor  (discussed  in  Section  6.3)  all  of  the  predicted  results 
were  well  within  the  accuracy  goals  for  the  program. 
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A breakdown  of  the  individual  losses  for  the  two  cases  m which  the  com- 
plete methodology  was  executed  (Extended  Delta  and  C-4  Third  Stage)  is  pre- 
sented in  Table  6-2.  The  results  ot  both  the  theoretical  and  simplified  methods 
are  indicated.  In  this  way  the  relative  magnitudes  of  the  different  losses  car.  be 
assessed.  The  erosion  loss,  normally  tncluded  in  (see  equation  3-2)  has 

been  listed  separately  for  the  same  reason.  In  order  to  facilitate  comparison  be- 
tween the  theoretical  and  simplified  2D-2phase  losses,  the  simplified  2D  and  2 
phase  losses  have  been  combined  (the  theoretical  loss  cannot  be  separated  into 
Us  components). 

As  a result  of  the  nozzle  throat  geometry  restrictions  of  the  TD2P  module, 
the  actual  nozzle  geometry  could  not  be  analyzed  for  the  three  high  performance 
motors:  Extended  Delta,  SRAM,  and  C4  Third  Stage.  The  transonic  flow  solution 
is  the  TD2P  module  is  limited  to  a simple  circular  arc  th.oat  geometry  and  the 
ratio  of  the  throat  radius  of  curvature  to  the  throat  radius  must  be  greater  than 
1.5.  The  two-phase  flow  loss  is  usually  the  largest  single  loss  for  highly  al- 
uminized propellants,  and  the  magnitude  of  this  loss  is  sensitive  to  the  throat 
geometry.  Therefore,  the  inability  of  the  TD2P  module  to  treat  the  actual  rozzle 
geometry  is  cause  for  some  concern.  An  effort  to  remove  this  limitation  would  be 
warranted.  In  the  mean  time,  the  magnitude  of  the  error  resulting  from  the  use 
of  a rationally  selected  approximate  nozzle  geometry  should  be  ascertained. 


! 
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BroaV-down  Of  Individual  Losses 


7.  CONCLUSIONS  AND  RECOMMENDATIONS 
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7.1  Conclusions 

The  goal  o £ this  effort  was  to  develop  a computerized  performance  predic- 
tion methodology  that  was  capable  of  predicting  delivered  specific  Impulse  to  +2%, 
and  delivered  thrust  and  total  Impulse  to  +5%,  for  a wide  range  of  solid  rocket 
motor  geometries , propellmt,  and  operating  conditions. 

The  Solid  Rocket  Performance  Program  (SPP)  described  herein,  and  in  the 
other  two  volumes  of  this  report,  was  designed  to  achieve  the  goal.  Six  existing 
•computer  programs  we  re  integrated  with  a central  logic  to  form  a strong  analytical 
foundation  for  the  program.  To  supplement  the  theory,  where  necessary,  and  to 
increase  the  flexibility  of  the  program,  a number  of  existing  and  newly  developed 
semi-empirical  correlations  were  Incorporated  into  the  program. 

The  program  is  modular  in  structure,  and  allows  the  user  to  exercise  the 
various  modules  individually,  all  together,  or  in  any  predetermined  combination, 
as  required.  For  those  cases  where  speed  and  simplicity  are  paramount  (e.g.  large 
parametric  studies}  the  analytical  portions  of  the  program  may  be  bypassed  completely 
in  favor  of  the  performance  loss  correlations  which  have  been  incorporated  into 
the  program.  Alternatively,  one,  or  more  of  these  correlations  may  be  employed 
In  lieu  of  the  corresponding  analytical  calculation  when  the  particular  performance 
loss  ii  elthe'  known  to  be  small,  or  adequately  characterized  by  the  correlation. 

Ballistics  and  performance  predictions,  obtained  with  the  SPP  Program,  were 
compare  with  firing  data  for  four  widely  varying  full  scale  motors.  Except  for  the 
anomalous  . .essure-tlme  behavior  of  one  of  the  motors  Can  explanation  for  which  is 
proffered),  allot  the  predicted  results  were  wi.iln  the  accuracy  goals  originally 
set  for  the  program.  In  each  case  the  difference  between  the  predicted  and  measured 
delivered  specific  impulse  was  less  tnan  one  percent.  This  Included  one  case  which 
was  predicted  to  within  0.4%  on  a sealed  envelope  basi3.  The  simplified  performance 
loss  correlations  also  gave  good  results  for  each  case.  Conclusive  statements  re- 
garding thr,  accuracy  and  range  of  validity  of  the  SPP  program  cannot  be  made  until 
additional  verification  efforts  are  conducted.  The  following  tentative  conclusions 
have  been  drawn: 
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»The  accuracy  of  burning  rate  prediction  Is  comparable  to,  and  will  be  limited 
wby,  the  error  In  measuring  strand  burning  rates. 

Action  time  average  thrust  and  pressure  appear  to  be  predictable  to  within 
•sit;  compensating  errors  In  the  various  parameters  entering  Into  the  compu- 
tations of  pressure  and  thrust  being  helpful  In  this  regard. 

^Delivered  specific  Impulse  appears  to  be  predictable  to  within  2%,  and  pro- 
bably to  within  lit  In  most  cases  which  do  not  Involve  unusual  circumstances. 

0The  simplified  method  of  computing  specific  Impulse  efficiency  Is  adequate 
for  Its  Intended  purpose. 

As  a result  of  the  limitations  of  the  present  analysis  (see  Section  7.2)  the 
accuracy  of  the  calculated  results  may  be  Influenced  by  compensating  errors.  For 
Instance,  the  geometric  limitation  of  the  TD2P  analysis  that  requires  the  actual 
nozzle  shape  to  be  approximated  In  many  cases  (Including  3 of  the  present  cases), 
the  possible  errors  In  particle  size  distribution,  the  neglect  of  throat  roughness, 
etc.,  all  Impact  the  calculated  two-phase  flow  loss.  No  attempt  to  quantify  the 
magnitude  of  the  possible  errors  resulting  from  these  and  other  assumptions  was 
made  during  the  course  of  this  Investigation.  Further  verification  efforts,  para- 
metric studies  and  future  experimental  and  analytical  investigations  will  undoubt- 
edly be  Illuminating  In  regards  to  the  adequacy  of  the  individual  assumptions  and 
the  program  as  a whole. 

The  SPP  program  treats  the  complete  solid  rochet  motor  performance  problem, 
I.e.,  grain  geometry,  ballistics  and  efficiency,  within  the  context  of  a single, 
automated,  entity.  As  a result,  accurate  analytical  predictions  of  solid  rocket 
motor  performance  can  be  obtained  In  a single  computer  run.  It  Is  hoped  that  the 
convenience  of  the  present  approach  will  help  stimulate  more  widespread  applica- 
tion of  theoretical  techniques  In  solid  rocket  motor  development  programs. 


Recommendations 


The  SPP  program  has  many  virtues;  it  also  has  its  limitations.  Recognition 
of  these  limitations  serves  a dual  purpose:  first,  to  allow  the  program,  as  it  exists, 
to  be  more  fruitfully  employed;  secondly,  to  yield  a better  perspective  on  the  nature 
of  the  Improvements  and/or  extensions  of  the  program  that  should  be  considered. 

The  following  Is  a list  of  the  more  significant  limitations  of  the  SPP  program. 
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1.  The  TD2P  currently  terminates  If  particle  Impingement  occurs  In  the 
nozzle. 

2.  The  two-dimensional  two-phase  transonic  analysis  Is  inadequate  In  a 
number  of  respects:  the  analysis  will  fall  for  normalized  radii  of  cur- 
vature <4.S  and  for  Inlet  angles  the  wall  geometry  Is  limited  to 
a conical  Inlet  and  single  radius  of  curvature  throat. 

3.  The  TD2P  program  also:  neglects  real  gas  effects;  cannot  handle 
shock  waves;  does  not  explicitly  allow  no  particle  cases  to  be  computed 
(the  program  can  be  “fooled"  Into  doing  so),  prohibits  gas  phase  con- 
densation and  particle  growth  In  the  nozzle. 

4.  The  present  boundary  layer  analysis  (TB1)  Is  not  deemed  to  be  accurate 
enough  when  boundary  layer  losses  are  relatively  largo. 

5.  A method  for  determining  nozzle  wall  temperature  has  not  been  Integra- 
ted Into  the  methodology. 

6.  The  nozzle  erosion  correlations  In  the  program  are  limited,  and  not  as 
accurate  as  would  be  desired. 

7.  Throat  roughness,  and  erosion  Induced  discontinuities  are  not  addressed. 

8.  The  analysis  Is  Inadequate  In  transient  dominated  cases. 

9.  The  grain  design  analysis  is  Inefficient  for  "simpler"  grain  geometries. 

10.  The  effects  of  acceleration  on  burn  rate,  and  nozzle  flow  separation 
have  not  been  considered. 

11.  Three  dimensional  flow  effects  have  not  been  considered. 

The  State-of-the-ort  In  solid  rocket  performance  prediction  was  recently 

(3) 

reviewed  at  the  last  JANNAF  Performance  Standardization  Working  Group  meeting  . 
As  3 result  of  this  meeting  seven  deficiency  areas  were  Identified.  These  areas 
Included  some,  or  all,  of  the  Items  listed  under  numbers  2,  3,  7 and  11,  above. 

In  addition,  they  Included  more  adequate  particle  size  determination,  and  better 
bum  rrtc  modeling.  So-called  "deficiency  statements"  were  recently  prepared  by 
several  of  the  comml.iee  members,  In  which  the  nature  of  these  limitations,  and 
methods  for  alleviating  them  are  ...plored  in  detail.  Based  on  these  statements, 
and  our  own  predilections,  it  Is  recommended  that  any  future  efforts  to  modify  the 
SPP  program  should  begin  with  items  1,  2 and  7,  as  listed  previously.  Item  1 Is  a 
relatively  simple  matter  to  take  care  of,  Items  2 and  7 would  require  significant 
effort,  but  are  entirely  possible  to  overcome  with  currently  available  numerical 
techniques. 
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Most  of  the  other  Items  In  this  list  represent  limitations  which  are  encount- 
‘ ' ered  less  frequently  (e.g. , large  boundary  layer  losses,  nozzle  separation,  tran- 

sient dominated  motors,  etc.),  or  ones  which  would  require  Inordinate  effort  to 
i overcome  (e.g.,  3-D  effects).  Some  of  these  limitations  can  currently  be  over- 

' come  through  the  use  of  existing  analyses  on  an  uncoupled  (to  SPP)  basis,  (e.g., 

: ] boundary  layer  and  tl  nrat  erosion) . Which  of  these  remaining  limitations  will  be 

f j addressed  will  most  likely  be  a function  of  the  predilections  of  future  users,  and 

- [ the  requirements  of  future  motors. 

In  passing,  one  other  recommendation  should  be  made.  Ihe  existing  know- 
ledge regarding  the  kinetics  of  typical  metallzed  solid  propellant  exhaust  species 
is,  as  pointed  out  In  Appendix  A,  meager,  at  best.  In  most  cases  this  Is  of  little 
consequence  In  conventional  aluminized  propellants,  as  the  kinetics  loss  tends  to 
be  small.  At  low  pressure,  however,  the  kinetics  loss  can  be  significant,  and  un- 
less additional  work  towards  Identifying  the  controlling  reactions,  and  their  rates. 
Is  forthcoming,  substantial  predictive  errors  can  be  expected  In  such  cases. 


i 


8.  REFERENCES 


1.  NASA  SP  8039,  "Solid  Rocket  Motor  Performance  Analysis  and  Prediction, " 

May  1971. 

2.  Cohen,  N.S.,etal,  “Solid  Rocket  Performance  Prediction  Techniques", 

AIAA/SAE  10th  Propulsion  Conference,  AIAA  Paper  No.  74-1200,  Oct.  1974. 

3.  JANNAF  Performance  Standardization  Working  Group,  Solid  Performance  Sub- 
committee, Minutes  of  July  30-31,  1974  meeting  at  NASA  Johnson  Space  Cen- 
ter, distributed  by  CPIA. 

4.  Gordon,  S.,  and  McBride,  B.J.,  "Computer  Program  for  Calculation  of  Com- 
ples  Chemical  Equilibrium  Compositions,  Rocket  Performance,  Incident 
and  Reflected  Shocks,  and  Chapman-Jouguet  Detonations,"  NASA  SP  273 
(1971). 

5.  Nickerson,  G.  R. , Coats,  D.  E. , and  Bartz,  J.  L. , “The  Two-Dimensional 

Kinetic  (TDK)  Reference  Computer  Program, " Engineering  and  Programming 
Manual,  Ultrasystems,  Inc.,  Dec.  1973,  prepared  for  Contract  No.  NAS9-  5 

12652,  NASA  JSC.  ' 

6.  Nickerson,  G.  R. , and  Kllegel,  J.  F. , “Axisymmetrlc  TWo-Phase  Perfect 
Gas  Performance  Programs,  Volume  1 Engineering  and  Programming  Descrip- 
tion” perpared  for  NASA  MSC  under  Contract  NAS9-4358,  April  1967. 

7.  Welngold,  H.  D. , “The  ICRPG  Turbulent  Boundary  Layer  (TBL)  Reference 

Program,"  prepared  for  ICRPG  perf.  Std.  Working  Group,  July  1968.  ( 

8.  Barron,  J.  G,,  Jr.,  Cook,  K.  S. , and  Johnson,  W.  C. , "Grain  Design  4 

and  Internal  Ballistics  Evaluation  Program,  Program  No.  64101,"  Hercules 

Powder  Co.,  Bachus  Works,  July  1967,  AD  818321.  * 

9.  Harry,  D.  P.,  Price,  C.  F.,  Small,  K.  R. , and  Taylor,  D.  E.,  "User's  i 

Manual  for  Nozzleless  Rocket  Motor  Internal  Ballistics  Comouter  Program, " 

AFRPL  TR-73-20,  March  1973. 

10.  Air  Force  Rocket  Propulsion  Laboratory . Theoretical  I Program,  User's 

Manual.  sp 

11.  Air  Force  Rocket  Propulsion  Laboratory,  Theoretical  I Program,  Listing 

ans  Sample  Cases.  sp 

12.  McBride,  B.  J.,  and  Gordon,  S, , "Fortran  IV  Program  for  Calculation  of 
Thermodynamic  Data,  ” NASA  TND  4097,  August  1967. 

13.  "Ballistic  Predictions  for  Mass-Augmented  Solid  Rocket  Motors",  AFRPL- 
TR-71-133,  Lockheed  Propulsion  Company,  Redlands,  CA.,  (Dec  1971). 

14.  Cohen,  N.  S.,  and  Derr,  R.  L.,  “Metal  Combustion  Efficiency  Improve- 
ment," 7th  TANNAF  Combustion  Meeting  (CPIA  Publication  204,  Vol  II, 
pp.  1-12,  Feb.  1971). 


8-1 


IS 


“Combustion  Instability  of  Aluminized  Propellants,"  NWC-TP-5060,  U.S 
Naval  Weapons  Center,  China  Lake,  CA.  (April  1971). 


16.  “Combustion  Tailoring  Criteria  for  Solid  Propellants,"  AFRPL-TR-69-190, 
Lockheed  Propulsion  Company,  Redlands,  CA.  (May  1969). 

17.  "Propellant  Tailoring  Techniques  for  Controllable  Rocket  Motors,"  AFRPL- 
TR-71-25,  Vol  n,  Lockheed  Propulsion  Company,  Redlands,  CA.  (Mar.  1971). 

18.  Northam,  G.  3.,  and  Sullivan,  E.  M.,  "Evaluation  of  Mg-Al  Eutectic  to 
Improve  Combustion  Efficiency  in  Low  Burning  Rate  Propellants, " 9th 
fANNAF  Combustion  Meeting  (CPLA  Publication  231,  Vol  I,  pp  1-10,  Dec. 
1972")'. 

19.  Kllegel,  J.  R. , et  a!..  "One-Dimensional  Two-Phase  Reacting  Gas  Non- 
equlllbrlum  Performance  Program,"  Vol.  1 Engineering  Analysis,  Report 
MSC-11780,  TRW  Systems  (Aug.  1967). 

20.  Svehla,  R.  A.,  "Estimated  Viscosities  and  Thermal  Conductivities  of  Gases 
at  High  Temperatures,"  NASA  TR-132,  1962. 

21.  Bird,  R.  B.,  Stewart,  U.  E.,  and  Ltghtfoot,  E.  N, , Transport  Phenomena. 

John  Wiley  & Sons,  1960. 

22.  Mason,  E.  A.,  and  Saxena,  S.  C.,  Physics  of  Fluids.  Vol.  1,  no.  5, 
pp.  361-369,  1958. 

23.  Nickerson,  G.  R. , et  al,  "Average  Gas  Properties  Computer  Program," 
Dynamic  Science  Report  No.  TR-024,  Dec.  1969. 

24.  Alber,  I.  E.,  "Comparison  and  Prediction  of  Computer  Program  Results 
for  Rocket  Engine  Performance  Prediction,  Chapter  V,  Boundary  Layer 
Friction  and  Heat  Transfer,"  Dynamic  Science  Report  No.  SN-82,  perpared 
for  1CRPG,  April  1968. 

25.  Plepcr,  I.  L.,  "1CRPG  Liquid  Propellant  Thrust  Chamber  Performance  Eval- 
uation Manual,"  CPXA  178,  Sept.  1968. 

26.  Back,  u,  H.,  Cuffel,  R.  F.,  and  Mossier,  P.  F.,  “Laminarlzation  of  a 
Turbulent  Boundary  Layer  in  Nozzle  Flow-Boundary  Layer  and  Heat  Transfer 
Measuiements  with  V/all  Cooling,  " ASME  Journal  of  Heat  Transfer,  Aug. 

1970,  (to  be  published). 

27.  Levine,  J.  N. , "Transpiration  and  Film  Cooling  Boundary  Layer  Computer 
Program,  Vol.  1.  Numerical  Solution  of  the  Turbulent  Boundary  Layer 
Equations  with  Equilibrium  Chemistry,”  Final  Report,  Contract  NAS7-79!, 
NASA  JPL,  June  1971. 

28.  Evans,  R.  M.,  and  Morse,  H.  L.,  “Interim  User's  Manual  for  Boundary 
Layer  Integral  Matrix  Procedure, " Version  J,  Aerotherm  Report  UM-74-41, 
prepared  for  NASA  MSC,  Contract  No.  NAS8-29667. 


8-2 


• i 

l * 


e, 

j 


l 


? 


29. 

30. 


31. 

32. 


33. 


34. 


35. 


36. 


37. 


38. 


39. 


40. 


Static  Test  Report,  ATS-EBM-1,  Jet  Propulsion  Laboratory,  Pasadena,  CA 
(May  1969). 

"156-Inch  Diameter  Motor  Liquid  Injection  TVC  Program,  " AFRPL-TR-66- 
109,  Lockheed  Propulsion  Company,  Redlands,  CA.  (1966);  "Rocket  Motor 
Manual"  (CPIA  Publication  M-l,  Unit  464). 

Aerojet  2.75  Final  Report,  AFRPL-TR-69-90,  (C). 

Bullard,  D.,  and  Kolden,  J.,  “Design  and  Development  of  2.75-Inch  Mark 
4 Mod  8 Rocket  Motor,  " Report  NAVWEPS  8647,  NOTS-TP-3688,  U.S. 
Weapons  Center,  China  Lake,  CA.  (April  1965);  "Rocket  Motor  Manual" 
(CPIA  Publication  M-l,  Unit  505). 

Broksbank,  R.  M.,  "Performance  of  Two  Lockheed  AGM-69A  Solid  Propellant 
Rocket  Motors  for  the  Short  Range  Attack  Missile  (SRAM)..."  Report  AEDC- 
TR-70-259,  Arnold  Air  Force  Station,  TN  (Dec.  1970). 

“Report  on  Static  Test  of  Hydac  Motor,  ” Motor  H22-3  (July  1964),  Motor 
H2S-8  (Oct.  1964),  Motor  H25-9  (Dec.  1964),  Lockheed  Propulsion  Com- 
pany, Redlands,  CA.;  “Rocket  Motor  Manual"  (CPIA  Publication  M-l,  Unit 
451). 


"Design  and  Development  of  the  High  Expansion  Ratio  Rocket  Motor,  “ 
AFRPL-TDR-64-31,  Lockheed  Propulsion  Company,  Redlands,  CA.  (Oct. 
1964);  "Rocket  Motor  Manual"  (CPIA  Publication  M-l,  Unit  400). 

Domal.  A.  F. , et  al;  "Research  and  Development  Testing  of  the  Mlnuteman 
III  Stage  III  Solid  Propellant  Rocket  Motor  ..."  Report  AEDC-TR-68-195, 
Arnold  Air  Force  Station,  TN  (Oct.  1968);  "Preliminary  Flight  Rating  Test 
of  the  Mlnuteman  III  Stage  III  Solid  Propellant  Rocket  Motor  ..."  Report 
AEDC-TR-68-249,  Arnold  Air  Force  Station,  TN  (Jan.  1969);  "Qualification 
Testing  of  the  Aerojet  Mlnuteman  III  (LGM-3CG)  Stage  in  Solid  Propellant 
Rocket  Motor  ..."  Report  AEDC-TR-70-89,  Arnold  Air  Force  Station,  TN 
(July  1970);  “Rocket  Motor  Manual”  (CPIA  Publication  M-l,  Unit  457). 

Merryman,  H.  L,,  "Performance  of  a UTC  FW-4S  Solid  Propellant  Rocket 
Motor  ..."  Report  AEDC-TR-68-253,  Arnold  Air  Force  Station,  TN  (Nov. 
1968);  "Rocket  Motor  Manual"  (CPIA  Publication  M-l,  Unit  522). 


Domal,  A.  F. , et  al,,  "Simulated  High  Altitude  Testing  of  a Fllghtweight 
TE-364-1  Solid  Propellant  Retrograde  Rocket  Motor  for  the  Surveyor  ...” 
Reports  AEDC-TR-64-260,  65-21,  65-41  and  65-50,  Arnold  Air  Force 
Station,  TN  (Jan-Apr  1965);  "Rocket  Motor  Manual”  (CPIA  Publication  M-l, 
Unit  361). 

Morash,  R.  T. , "Final  Report  on  Rocket  Motor  ...  Antares  I,"  Report  A8L/ 
R-S3,  Allegheny  Ballistics  Laboratory,  Hercules,  Inc.,  Cumberland,  MD. 
(Nov.  1963). 


"Rocket  Motor  Manual" 


(CPIA  Publication  M-l,  Unit  421). 


1 


L, 


■ ... 


i ' 

i : 


A 


3-3 


p 


| 


* i 

•I 


& 


41.  "Sparrow  in  Atm  7F  Pre-Flight  Rating  Test  Report,"  Raytheon  Co.,  Bedford, 
MA.  (April  1968);  AD  514  402L;  "Rocket  Motor  Manual"  (CPIA  Publication 
M-l,  Unit  476,  477). 

42.  "Rocket  Motor  Manual"  (CPIA  Publication  M-l,  Unit  510). 

43.  "Rocket  Motor  Manual"  (CPIA  Publication  M-l,  Unit  467). 

44.  Beckman,  C.  W. , "Solid  Propellant  Impulse  Scaling  Prediction  Techniques,  " 
AFRPL-TR-71-7,  Air  Force  Rocket  Propulsion  Laboratory,  Edwards,  CA.  (Feb. 
1971). 

45.  "Study  of  a Fluid-Controlled  Solid  Rocket  Motor  for  a Mars  Orblter,  ” Report 
LPC-929Fon  Contract  NAS7-701,  Vol.  II,  Lockheed  Propulsion  Company, 
Redlands,  CA.  (Feb.  1969);  also  cf.  “Fluorine-Hydrogen  Performance  Eval- 
uation," Report  NASA  CR-72038,  Rocketdyne  Division  of  Rockwell  Ir.  er- 
natlonal  Corp. , Canoga  Park,  CA.  (April  1967). 

46.  Coluccl,  S.  E.,  “Experimental  Determination  of  Nozzle  Heat  Transfer  Co- 
efficient," TP106SRP,  Aerojet  Solid  Propulsion  Company,  Sacramento,  CA. 
(May  1960). 

47.  Grover,  S.  S.,  "Analysis  of  Nozzle  Heat  Transfer  Coefficient,”  TM  113 
SRP,  Aerojet  Solid  Propulsion  Company,  Sacramento,  CA.  (April  1960). 

48.  Collection  of  Internal  memoranda  from  L.  J.  Gordon,  J.  P.  Goughlln  and 
R.  M.  Smith,  Aerojet  Solid  Propulsion  Company,  Sacramento,  CA.  (1960- 
1964). 

49.  “Development  of  a Specific  Impulse  Scaling  Technique  for  Smokeless  Pro- 
pellants," Technical  Report  S-286,  Rohm  & Haas  Co.,  Huntsville,  AL. 

(Ncv.  1970). 

50.  Smith,  P.  W.  and  Stetz,  G.  A.,  "Impulse  Scaling  Prediction,"  AFRPL-TR- 
66-297,  Air  Force  Rocket  Propulsion  Laboratory,  Edwards,  CA.  (Nov.  1966); 
also,  cf.  Ref  (3). 

51.  Kllegel,  J.  R. , "Gas  Particle  Nozzle  Flows."  9th  International  Combustion 
Symposium  (Academic  Press,  New  York,  1963)  811-826. 

52.  Kllegel,  J.  R.  and  Nickerson,  G.  R.,  "Flew  of  Gas-Particle  Mixtures  In 
Axially  Symmetric  Nozzles, " ARS  Progress  In  Astronautics  and  Rocketry. 

Vol.  6:  Detonation  and  Two-Phase  Flow  (Academic  Press.  New  York.  1962) 
173-194. 

53.  Kordlg,  J.  W.  and  Fuller,  G.  H. , "Correlation  of  Nozzle  Submergence 
Losses  In  Solid  Rocket  Motors,”  J.  AIAA  j>,  175-177  (Jan.  1967). 

54.  "Nozzleless  Rocket  Motor  Internal  Ballistics  Computer  Program,"  AFRPL- 
TR-73-19,  Lockheed  Propulsion  Company,  Redlands,  CA.  (Mar.  1973). 


8-4 


55.  “Hydroquench  Thrust  Termination  Analysts , " AFRPL-TR-73-62,  Lockheed 
Propulsion  Company,  Redlands,  CA.  (Aug.  1973). 

56.  "User’s  Manual  - Aerotherm  Charring  Material  Thermal  Response  and  Ablation 
Program,  Version  3,  " AFRPL-TR-70-92,  Aerotherm  Corp.,  Mt.  View,  CA. 

(April  1970). 

57.  “User's  Manual  - Aerotherm  Axi-Symmetric  Transient  Heating  and  Material 
Ablation  Computer  Program  (ASTHMA  3),"  AFRPL-TR-72-24,  Aerotherm  Corp. 
Mt.  View,  CA.  (Jan.  1972). 

58.  Courtney,  W.  G. , "Review  of  Particle  Size  Determination  Techniques," 
AFRPL  Two-Phase  Flow  Conference.  AFRPL-TR-67-223,  Vol.  I,  Air  Force 
Rocket  Propulsion  Laboratory,  Edwards,  CA.  (Mar.  1967). 

59.  Brown,  B.  and  McCarty,  K.  P.,  "Particle  Size  of  Condensed  Oxide  from 


60.  Sehgal,  R. , "An  Experimental  Investigation  of  a Gas-Particle  System, " 
TR-32-238,  Jet  Propulsion  Laboratories,  Pasadena,  CA.  (Mar.  1962). 

61.  Cheung,  H.  and  Cohen,  N.  S. , "Performance  of  Solid  Propellants  Containing 
Metal  Additives,"  J.  AJAA  3j  250-257  (Feb.  1965). 

62.  Smith,  P.,  et  al.,  "Summary  Results  of  Particle  Size  Measurements, " 

AFRPL  Two-Phase  Flow  Conference.  AFRPL-TR-67-223,  Vol.  I,  Air  Force 
Rocket  Propulsion  Laboratory,  Edwards,  CA.  (Mar.  1967). 

63.  Collins,  R.  G.,  "The  AFRPL  Ballistic  Test,  Evaluation  and  Scaling  Program," 
AFRPL-TR-68-198,  Air  Force  Rocket  Fropulslon  Laboratory,  Edwards,  CA. 
(May  1965). 

64.  Stetz,  G.A.,  etal.,  "The  AFRPL  Ballistic  Test,  Evaluation  ar.d  Scaling 
Program,"  AFRPL-TR-68-198,  Air  Force  Rocket  Propulsion  Laboratory, 
Edwards,  CA.  (Nov.  1968). 

65.  Brown,  B. , "Oxide  Particles  In  Solid  Rocket  Exhausts,"  ARFPL  Two-Phase 
Flow  Conference,  AFRPL-TR-67-223,  Vol.  I,  Air  Force  Rocket  Propulsion 
Laboratory,  Edwards,  CA.  (Mar.  1967). 

66.  Crowe,  C.  T.  and  Wllloughly,  P.  G. , "A  Study  of  Particle  Growth  in  A 
Rocket  Nozzle,"  J.  AIAA  5,  1300-1304  (July  1967). 

67.  Thiokol  Chemical  Corp.,  Elkton  Division,  Elkton,  Maryland  "TE-M-364-4 
Extended  Delta  Rocket  Motor"  Final  Report  NASA  Contract  NAS7-678,  1972. 

68.  Dalnes,  W.  L.,  Mayberry,  J.  L.,  Lur.d,  R.  K.,  and  Abel,  R.,  "Prediction 
of  Thrust  Losses  Occurrlnq  in  Solid-Propellant  Rocket  Motors,"  AIAA^SAE 
10th  Propulsion  Conference,  AIAA  Paper  No.  74-1201,  Cct.  1974. 

69.  Lavery,  J.R.,  Abel,  R.,  Lund,  R.  K. , and  Dalnes,  W.  L. , "Experimental 
Investigation  of  the  Effects  of  Internal  Nozzle  Contouring  on  Nozzle  Perfor- 
mance," AIAA/SAE  10th  Propulsion  Conference,  AIAA  Paper  No.  74-1202, 

Oct.  1974. 


Combustion  of  Metallized  Solid  Propellant,"  8th  International  Combustion 
Symposium  (Williams  and  Wilkins,  Baltimore,  MD.  1962),  814-823. 


^1 


(' 

i 


Appendix  A - Reaction  Rate  Screening 

The  reaction  rate  screening  done  under  this  contract  should  by  no  means 
be  considered  complete  either  from  an  absolute,  or  from  a relative  basis.  What 
was  attempted  was  to  develop  a preliminary  reaction  rate  set  which  would  allow 
for  the  accurate  prediction  of  the  losses  In  performance  of  a solid  rocket  motor 
which  are  due  to  finite  rate  gas  phase  chemical  kinetics. 

The  criteria  by  which  the  reaction  rate  screening  was  done  was  to  select 
the  minimum  number  of  species  and  reactions  which  would  produce  an  acceptably 
accurate  value  of  Isp  efficiency  ils  efficiency  Is  defined  as 

spKIN 

"KIN  ~ T 


where  I = I as  calculated  In  the  presence  of  finite  rate  chemistry 
spKIN  air  calculated  by  the  ODK  module. 

I =1  as  calculated  by  the  restricted  equilibrium  option  of 
spRE  tffi  ODE  module 

However,  at  the  time  the  reaction  set  screening  was  done,  the  restricted 
equilibrium  option  for  the  ODE  module  was  not  completed.  This  was  unfortunate 
since  the  difference  In  Igp  as  calculated  by  the  two  options  In  the  ODE  module 
can  exceed  I % and  thus  Influence  the  number  of  species  and  reactions  which  could 
dropped  from  reaction  set. 

The  basic  steps  in  developing  »ha  master  reaction  set  for  alumlnumlzed 
propellants  and  the  screening  of  that  reaction  set  are  outlined  below. 

1.  The  reaction  set  of  Reference  A1  (AP-PBAN  with  no  Al)  was  merged 

with  those  reactions  in  Reference  A2  (which  Included  aluminized  species). 
This  process  yielded  a Initial  master  reaction  rate  set  of  31  species 
and  126  reactions. 

2.  The  Initial  reaction  rate  set  was  compared  with  the  results  of  chemical 
equilibrium  calculations  and  was  found  to  be  deficient. 

3.  A short  literature  search  was  conducted  for  measured  rates  concerning 
aluminized  species.  No  measured  rate  data  was  found. 

4.  Estimates  were  made  for  rates  Involving  6 more  species  and  30  more 
reactions  using  the  procedure  reported  In  Reference  Al  (the  same  pro- 
cedure was  used  In  Reference  A2  to  estimate  the  reaction  rates  of 
aluminized  species). 


A-l 


j 

I 


; 


I 


a 


i 

i 


5.  Using  the  Generalized  Kinetics  Analysis  Program  (GKAP)  of  Reference 
A3,  the  total  master  reaction  set  was  screened.  This  computer  program 
SCREENED  the  reaction  in  an  absolute  sense.  That  is,  only  reactions 
which  had  negligible  influence  on  the  fluid  flow  or  rate  of  production 
of  a chemical  species  everywhere  in  the  nozzle  expansion  were  drop- 
ped. The  total  reaction  set  was  reduced  only  minimally  by  this  first 
step  in  the  screening  process. 

6.  The  next  step  in  screening  procedure  was  to  reduce  the  number  of  species 
and  hence  the  number  of  reactions  by  making  an  ODE  run  omitting  species 
thought  do  dbly  to  be  unimportant  and  comparing  the  results  to  an  ODK 
calculation  using  the  reduced  reaction  s ' . This  procedure  reduced 

the  reaction  set  to  23  species  and  73  reactions. 

While  the  above  screening  procedure  reduced  both  the  number  of  species 
and  reactions  considered  in  performance  calculations  by  the  SPP  code,  it  is  by  no 
means  considered  adequate.  Too  many  species  and  reactions  were  Included  in 
this  set,  which  are  due  to  the  comparison  of  two  Inconsistent  physical  models  of 
the  expansion  of  a chemically  reacting  gas  through  a rocket  motor  nozzle.  Thus, 
the  screening  done  to  date  should  not  be  considered  accurate  until  measured  data  are 
taken  on  the  reactions  considered  in  this  set  and/or  performance  losses  due  to 
finite  rate  gas  phase  chemical  kinetics  can  be  quantified. 

The  final  reaction  set  selected  for  conventional  aluminized  propellants  is 
shown  in  Table  2-13  of  Volume  III.  The  following  notes  concerning  the  comment 
field  of  this  reaction  set  are  Intended  to  help  the  reader  understand  where  these 
reactions  and  rate  data  came  from. 

• The  first  part  of  each  section  (third  body  and  binary  exchange  reactions) 
is  from  Reference  Al.  For  example,  BAULCHO 968)1.1,  refers  to  the  mea- 
surements made  by  Baulch  in  1968  and  reported  in  LEEDS  1 report, 

•The  second  part  Is  from  Reference  A2  and  refers  to  the  reaction  number 
built  Into  that  computer  program. 

•The  third  part,  always  suffixed  by  the  word  ESTIMATE,  is  from  reactions 
and  rates  estimated  by  the  authors  of  this  document. 
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